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Abstract 


Nonlinear  optical  processes  provide  a  physical  mechanism  for  converting  the  frequency 
of  light.  This  allows  the  generation  of  tunable  light  sources  at  wavelengths  inaccessible 
with  lasers,  leading  to  a  diverse  set  of  applications  in  helds  such  as  spectroscopy, 
sensing,  and  metrology.  To  make  these  processes  efficient  has  conventionally  required 
relatively  exotic  materials  that  are  incompatible  with  state  of  the  art  nanofabrication, 
resulting  in  large-area  devices  that  operate  at  high  optical  powers  and  cannot  be 
integrated  with  on-chip  optical  and  electronic  circuits. 

This  dissertation  shows  how  optical  nanocavities,  by  localizing  light  into  sub-cubic 
optical  wavelength  volumes  with  long  photon  storage  times,  can  greatly  enhance  the 
efficiency  of  nonlinear  frequency  conversion  processes  in  III-V  semiconductors,  while 
simultaneously  shrinking  the  device  footprint,  reducing  the  operating  power,  and 
providing  a  scalable  on-chip  platform.  This  approach  also  enables  on-chip  quantum 
frequency  conversion  interfaces,  which  are  crucial  for  the  construction  of  quantum 
networks. 

First,  photonic  crystal  nanocavities  in  gallium  phosphide  are  shown  to  generate 
second  harmonic  radiation  with  only  nanowatts  of  coupled  optical  powers,  and  effi¬ 
ciency  many  orders  of  magnitude  greater  than  in  previous  nanoscale  devices.  This 
approach  is  then  extended  to  demonstrate  sum-frequency  generation  in  GaP  pho¬ 
tonic  crystal  cavities  with  multiple  cavity  modes,  as  well  as  broadband  upconversion 
employing  photonic  crystal  waveguides.  The  nanocavity-enhanced  second  harmonic 
generation  is  then  integrated  with  a  single  quantum  dot  to  create  a  single  photon 
source  triggered  at  300  MHz  by  a  telecommunication  wavelength  laser  coupled  with 
an  external  electro-optic  modulator,  a  simpler  and  faster  conhguration  than  standard 


approaches. 

The  efficiency  of  all  the  aforementioned  processes  can  be  further  improved  through 
resonant  photonic  crystal  nanocavities  allowing  large  frequency  separation,  which  are 
described  in  this  thesis.  Finally,  this  dissertation  presents  spectroscopic  measurements 
of  quantum  systems  that  emit  photons  at  visible  wavelengths,  which  are  promising 
candidates  for  a  number  of  quantum  and  classical  applications  and  well-suited  for 
integration  with  on-chip  frequency  conversion. 
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are  plotted  at  the  X  point).  Between  the  frequencies  of  the  dielectric 
band  and  air  band,  there  are  no  guided  modes  in  the  slab.  Boundary 
between  gray  and  white  regions  indicates  the  light  line.  Image  courtesy 
of  Jelena  Vuckovic . 


3.4  (a)  SEM  image  of  an  L3  cavity,  (b)  Local  electric  field  intensity  inside 

the  cavity  at  the  resonant  frequency.  .  21 

3.5  (a)  SEM  image  of  an  nanobeam  cavity,  (b)  Ey  inside  the  beam  cavity.  21 


4.1  (a)  Scanning  electron  microscope  image  of  a  fabricated  GaP  photonic 
crystal  membrane  after  undercut  of  sacrificial  layer  with  lattice  con¬ 
stant  a=246  nm,  hole  radius,  r/a=0.32,  slab  thickness  (i/a=0.57.  (b) 

Finite  difference  time  domain  simulation  of  electric  held  intensity  inside 

the  cavity  for  the  high  Q  mode .  25 

4.2  (a)  Scanning  electron  microscope  image  of  a  1  /rm-wide  trench  etched 

in  GaP.  (b)  SEM  of  200  nm-wide  trenches,  (c)  SEM  of  100  nm-wide 
trenches .  26 

4.3  SEM  images  of  chips  undercut  a  sufficient  amount  of  time  to  also  etch 

GaP  membrane  photonic  crystal  holes,  transforming  circular  photonic 
crystal  holes  into  squares .  26 

4.4  (a)  Experimental  setup.  The  V-polarized  cavity  is  probed  in  a  cross- 
polarized  setup  using  a  PBS  and  a  HWP.  The  cavity  signal  observed 
at  the  output  follows  a  sin (46*)  dependence  where  theta  is  the  HWP 
setting  (6*=0  corresponds  to  V  polarization),  (b)  Spectrum  of  a  cavity 
measured  with  conhguration  of  Fig.  4.4(a)  with  HWP  at  6*  =  67°. 

The  cavity  resonance  is  at  722  nm  at  room  temperature.  A  £t  to  a 
Lorentzian  (solid  line)  gives  Q=1100.  (d)  (Gircles)  Normalized  counts 
in  cavity  peak  as  a  function  of  HWP  angle  (background  subtracted). 

(Line)  Fit  with  sin(46*) .  27 

4.5  (a)  Shift  of  photonic  crystal  resonance  at  room  temperature  as  a  and 

r/a  are  changed.  Solid  lines  are  fits  to  Lorentzians.  (b)  Gavity  res¬ 
onance  measured  at  room  temperature  (blue)  and  low  temperature 
(red).  Solid  lines  indicate  fits  to  Lorentzians .  28 
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5.1 


(a)  Microscope-based  setup  for  second  harmonic  generation.  HWP: 
half  wave  plate,  NPBS:  nonpolarizing  beamsplitter,  OL:  objective  lens, 

PBS:  polarizing  beamsplitter,  SPF:  short  pass  hlter,  PD:  photodiode. 

The  incident  light  traces  the  red  line  into  the  cavity  sample.  The  second 
harmonic  light  follows  the  blue  line  into  the  spectrometer,  photodiode, 
or  camera.  The  polarization  of  the  incident  light  is  controlled  by  the 
polarizer  and  HWP;  the  polarization  of  the  second  harmonic  radiation 
is  measured  using  HWP  and  PBS.  (b)  SEM  image  of  a  fabricated 
structure.  Scale  bar  indicates  1  /rm.  (c)  Spectrum  of  generated  second 
harmonic  light  with  8  nW  power  at  1497.4  nm  coupled  to  the  cavity 
(160  nW  incident) .  33 

5.2  (a)  Finite  difference  time  domain  (FDTD)  simulation  of  electric  held 

inside  the  cavity  for  the  fundamental  TE-like  cavity  resonance  in  the 
center  of  the  slab.  Cavity  held  axes  are  E'^  and  Ey  (b)  Illustration  of 
orientation  of  cavity  relative  to  crystal  axes.  Cavities  axes  Ej,,  E'y  are 
rotated  from  crystal  axes  Ex,  Ey  by  an  angle  9.  Fields  along  crystal 
axes  are  determined  by  projection .  34 

5.3  (a)  FDTD  simulation  of  TM-like  photonic  bands  for  the  same  triangu¬ 

lar  lattice  photonic  crystal.  Red  indicates  band  positions.  White  solid 
lines  indicate  light  line;  black  solid  lines  indicate  numerical  aperture 
of  lens.  White  box  indicates  mode  at  second  harmonic  frequency,  a: 
lattice  constant  of  photonic  crystal,  (b)  E^  held  patterns  of  degenerate 
TM-like  mode  at  second  harmonic  frequency  at  the  F  point .  35 

5.4  (a)  Spectrum  of  fundamental  resonance  probed  in  cross-polarized  re- 

hectivity  with  a  broadband  source.  Lorentzian  ht  gives  a  quality  factor 
of  5600.  (b)  Spectrum  at  second  harmonic  as  exciting  laser  frequency 
is  tuned  across  the  cavity  resonance.  Solid  line  shows  ht  to  Lorentzian 
squared  with  cavity  quality  factor  of  6000 .  36 
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5.5  (a)  Second  harmonic  power  as  a  function  of  fundamental  wavelength 
power  coupled  into  the  cavity.  We  estimate  coupling  efficiency  into  the 
cavity  to  be  5%.  Solid  line  indicates  £t,  with  slope  2.02,  indicating 
quadratic  power  dependence.  Output  power  measurements  are  cor¬ 
rected  for  measured  losses  from  optics,  but  do  not  include  corrections 

for  collection  efficiency  into  the  objective  lens .  38 

5.6  (a)  Dependence  of  second  harmonic  power  on  incident  light  polariza¬ 

tion.  The  horizontal  axis  corresponds  to  the  angle  between  the  input 
polarization  and  the  polarization  of  the  cavity  mode.  Solid  line  shows 
fit  to  cos^(0).  (b)  and  (c)  Second  harmonic  radiation  imaged  on  a  cam¬ 
era  with  polarizer  oriented  in  x'  (b)  and  y'  (c)  direction  [the  orientation 
of  cavity  relative  to  the  axes  is  shown  in  Fig.  5.2(b)].  The  gray  box 
indicates  the  approximate  location  of  the  photonic  crystal  structure, 
with  scale  bar  indicating  approximately  5  pm.  Camera  integration 
time  for  the  same  for  both  images .  39 

5.7  (a)  Scanning  electron  microscope  image  of  a  fabricated  photonic  crys¬ 
tal  membrane  after  undercut  of  sacrihcial  layer.  Scale  bar  indicates 
1  pm.  (b)  White  light  reflectivity  from  cavity,  with  FDTD  simula¬ 
tions  showing  electric  held  of  modes  of  interest.  The  two  modes,  both 
transverse  electric-like  (TE-like)  are  polarized  orthogonally:  the  fun¬ 
damental  mode  is  y-polarized;  the  higher  order  mode  is  x-polarized. 

(c)  Fit  to  Lorentzian  of  fundamental  mode  (peak  indicated  by  black 
box)  with  Q=3800 .  43 
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5.8  (a)  Schematic  of  sum  frequency  generation  process.  Incident  photons 

with  frequencies  Ui  and  U2  are  converted  to  frequency  (b)  Exper¬ 
imental  setup  for  SFG.  HWP:  half  wave  plate,  NPBS:  nonpolarizing 
beamsplitter,  OL;  objective  lens.  The  incident  light  traces  the  red  line 
into  the  cavity  sample;  the  sum  frequency  light  follows  the  blue  line 
into  the  spectrometer.  The  polarization  of  the  incident  light  is  con¬ 
trolled  by  a  polarizer  and  HWP.  (c)  Simulated  electric  held  pattern  of 
degenerate  TM-like  Bloch  modes  of  the  crystal  at  wavelength  of  sum 
frequency .  44 

5.9  (a)  Sum-frequency  generation  and  second  harmonic  generation  when 
two  lasers  are  used  to  couple  into  modes  shown  in  Fig.  5.7.  Middle 
peak  (green  box)  is  sum  frequency  generation;  right  (black  dotted  and 
red  boxes)  and  left  peaks  are  second  harmonic  generation  from  each 
mode.  Black  box  indicates  second  harmonic  from  fundamental  mode 
over  which  laser  is  scanned  to  measure  power  and  wavelength  depen¬ 
dence  in  Fig.  5.9b  and  5.9c.  (b)  Power  dependence  of  second  harmonic 
and  sum  frequency  generation  as  a  function  of  incident  laser  power 
at  wavelength  of  fundamental  cavity  mode.  A  polynomial  ht  (plotted 
on  log-log  scale)  gives  a  coefficient  of  0.9  for  sum-frequency  generation 
(green  squares)  and  1.8  for  second  harmonic  (red  circles).  The  incident 
power  of  the  second  laser  at  higher  order  cavity  mode  wavelength  is  6 
mW.  (c)  Wavelength  dependence  of  sum  frequency  (green  squares)  and 
second  harmonic  counts  (red  circles)  as  a  function  of  scanning  laser 
wavelength  at  fundamental  cavity  mode.  Fit  to  Lorentzian  of  SFG 
counts  gives  Q  of  4300;  ht  to  Lorentzian  squared  of  SHG  counts  gives 
Q  of  4100,  both  matching  Q  of  the  fundamental  mode  in  Fig.  5.7,  as 
expected.  The  other  laser  is  at  1504  nm.  (d)  Wavelength  dependence 
of  sum  frequency  counts  as  a  function  of  scanning  laser  wavelength  at 
higher  order  cavity  mode.  Fit  to  Lorentzian  of  SFG  counts  gives  Q  of 
272,  matching  Q  of  the  higher  order  mode.  Second  laser  is  at  1567.4 

nm .  45 
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5.10  (a)  Dispersion  diagram  for  TE-like  modes  (E- field  at  the  center  of  the 
slab  in  plane)  for  W1  waveguide.  Parameters  are:  hole  radius  r/a=0.25 
and  thickness  d/ a=0.25  where  a  is  periodicity  in  x  direction.  The  white 
area  between  the  gray  shaded  regions  indicates  the  photonic  band  gap 
of  the  triangular  lattice  photonic  crystal.  The  green  line  indicates 
the  waveguide  mode  with  even  symmetry  of  the  held  component 
relative  to  the  x-z  plane  (including  waveguide  axis).  The  blue  line 
indicates  waveguide  mode  with  odd  symmetry.  Circles  indicate  FDTD- 
calculated  solutions,  and  solid  lines  indicate  interpolated  bands.  Red 
line  indicates  light  line.  The  inset  shows  the  FDTD  simulation  of  the 
Bz  held  component  at  the  center  of  the  slab  for  the  band  plotted  in 
green  at  the  kx  =  n/ a  point  (circled)  with  frequency  a/A=0.36.  (b) 

SEM  image  of  30-periodic  PC  photonic  crystal  waveguide  fabricated 
in  160  nm  thick  GaP  membrane  with  a=560  nm.  Scale  bar  indicates  2 
pm.  Circular  grating  at  left  or  modihed  holes  forming  coupler  (right) 

can  be  used  to  couple  into  waveguide  modes  from  free  space .  51 

5.11  (a)  Transmission  measurement  through  photonic  crystal  waveguide  us¬ 
ing  white  light  source,  (b)  Group  index  calculated  from  FDTD  sim¬ 
ulation  (Fig.  5.10a)  and  experimental  data  (Fig.  5.11a).  Maximum 
group  index  measurement  experimentally  is  25.  (c)  Second  harmonic 
intensity  measured  from  a  diherent  structure  as  a  function  of  incident 
laser  wavelength.  Inset:  second  harmonic  counts  at  the  FP  peaks  and 
calculated  as  a  function  of  wavelength,  (d)  Second  harmonic  counts 
measured  as  a  function  of  incident  laser  power.  Red  line  indicates  lin¬ 
ear  ht  of  log-log  data  with  slope  1.9.  The  laser  wavelength  is  1561.5 

nm;  the  calculated  group  index  at  this  wavelength  is  30 .  52 
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5.12  (a)  FDTD  simulation  of  Ez  field  component  of  the  monopole  TM-like 
mode  near  frequency  of  second  harmonic  in  triangular  photonic  crystal 
lattice.  (a/A=0.67)  (b)  Calculated  far-field  radiation  pattern  of  mode 
in  (a).  White  circle  indicates  numerical  aperture  of  lens,  (c)  FDTD 
simulation  of  Ez  in  photonic  crystal  waveguide  at  second  harmonic 
frequency  (a/A=0.71).  (d)  Calculated  far-field  radiation  pattern  of 
\Ex\'^  for  mode  in  (c).  (e)  Calculated  far-field  radiation  pattern  oi\Ey\'^ 
for  mode  in  (c).  (f)  Calculated  far-field  radiation  pattern  of  \E\^  for 
mode  in  (c).  (g)  Measured  radiation  pattern  for  \Ex\‘^.  (h)  Measured 
radiation  pattern  for  \Ey\^.  (i)  Measured  radiation  pattern  \E\^.  ...  53 

6.1  Photonic  crystal  cavity  design  and  fabrication,  (a)  SEM  image  of  fab¬ 
ricated  suspended  membrane  photonic  crystal  cavity.  Scale  bar  indi¬ 
cates  1  yum.  (b)  Reflectivity  measurement  of  fundamental  cavity  mode. 
Lorentzian  fit  gives  Q  of  7000.  (c)  (i)  and  (ii)  Simulated  electric  field 
components  for  fundamental  mode  of  L3  photonic  crystal  cavity,  used 
for  resonantly  enhanced  upconversion.  (d)  (i)  and  (ii)  Electric  field 
patterns  for  TE  mode  closest  to  emission  frequency  of  quantum  dot. 

The  low  Q  (<100)  mode  is  formed  through  weak  confinement  of  an  air 
band  mode .  56 

6.2  Characterization  of  single  quantum  dot  excited  with  on-chip- up  converted 

1500  nm  laser,  (a)  Spectrum  measured  from  CW  second  harmonic  ex¬ 
citation  of  quantum  dot.  Inset:  Streak  camera  measurement  of  quan¬ 
tum  dot  lifetime.  Counts  were  summed  over  spectral  window  indi¬ 
cated  by  box.  (b)  Photon  correlation  measurement  of  quantum  dot 
emission  under  frequency  doubled  CW  1550nm  excitation.  Fit  gives 
g(2)(o)=0.43±0.04 .  60 
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6.3  Experimental  setup.  A  tunable  C-band  laser  (red  excitation  path)  is 

modulated  by  an  electro-optic  modulator  (EOM)  driven  by  a  pulse 
pattern  generator  (PPG),  which  receives  an  external  clock  signal  from 
a  synthesizer.  The  laser  polarization  is  adjusted  to  match  that  of  the 
cavity  by  a  half  wave  plate  (HWP).  Pulses  incident  on  the  cavity  via  a 
non-polarizing  beamsplitter  (BS)  excite  quantum  dots  by  above  band 
absorption  of  laser  pulses  upconverted  through  second  harmonic  gen¬ 
eration  in  the  surrounding  GaAs  matrix.  Emission  from  the  quantum 
dot  (green  path)  is  spectrally  hltered,  then  analyzed  either  by  an  HBT 
setup  for  photon  correlation  measurements  or  a  spectrometer.  For  di¬ 
rect  lifetime  measurements  of  the  quantum  dot,  an  ultrashort  pulse  (3 
ps)  Ti;Sapphire  laser  at  80  MHz  repetition  rate  excites  the  dot  (blue 
excitation  path);  a  streak  camera  is  used  for  detection.  Inset  shows 
scanning  electron  microscope  image  of  photonic  crystal  cavity .  61 

6.4  Schematic  illustrating  principle  of  fast  up  conversion-based  telecom  wavelength- 
triggered  quantum  dot  single  photon  source.  Pulses  at  1500  nm  are 
coupled  into  the  cavity  resonance,  leading  to  strongly  enhanced  second 
harmonic  generation.  The  light  at  the  second  harmonic  frequency  can 
create  carriers  in  the  GaAs  membrane  that  relax  into  the  quantum 

dot,  which  are  then  emitted  as  single  photon  pulses.  The  rate  of  single 
photon  pulses  is  determined  by  the  rate  at  which  the  laser  is  modulated.  62 

6.5  Photon  correlation  measurement  of  quantum  dot  emission  when  trig¬ 

gered  with  frequency  doubled  telecom  wavelength  laser  pulses  from  an 
externally  modulated  telecom  wavelength  laser,  (a)  Second  order  au¬ 
tocorrelation  function  measurement  for  100  MHz  repetition  rate  with 
duty  cycle  20%.  g*'^^(0)  =  0.48  <  0.5  indicates  emission  from  a  sin¬ 
gle  quantum  dot.  (b)  Second  order  autocorrelation  measurement  for 
repetition  rate  300  MHz  with  duty  cycle  50%  and  g^^^(0)=0.40 .  63 
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7.1  (a)  Schematic  illustration  of  multiply  resonant  orthogonal  nanobeam 

cavity,  (b)  Illustration  of  photonic  nanobeam.  Red  box  shows  unit 
cell,  which  is  tiled  periodically  in  the  x  direction.  Parameters  are  pe¬ 
riodicity  a,  width  w,  hole  size  hx  and  hy  and  thickness  out-of-plane 
t  (not  shown),  (c)  Wavelengths  of  dielectric  and  air  bands  of  GaAs 
nanobeam  as  lattice  constant  is  varied.  Parameters  are  w/a=1.65, 
hi/w=0.6,  ^2/0=0. 5.  Solid  lines  are  plotted  for  t=160  nm  (hxed  ab¬ 
solute  slab  thickness)  and  wavelength-dependent  index  of  refraction 
n;  dotted  lines  are  plotted  for  t/a  =  0.35  (hxed  relative  slab  thick¬ 
ness)  and  n=3.37.  Field  patterns  of  Ey  at  kx  =  n/a  for  dielectric 
(top)  and  air  band  (bottom)  are  indicated  by  black  circles  and  arrows, 
(d)  Normalized  freqnencies  of  dielectric  and  air  bands  (left  axis)  of 
GaAs  nanobeam  as  beam  width  is  varied.  Parameters  are  hi/w=0.6, 
^2/0=0. 5,  t/a=0.35.  Right  axis  shows  change  in  size  of  photonic  band 
gap  with  beam  width . 
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7.2  (a)  Illustration  of  cavity  design,  showing  intersecting  orthogonal  nanobeams 
with  taper  and  mirror  regions,  as  well  as  central  cavity,  (b)  Detail  of 
white  box  in  (a).  Parameters  used  to  form  resonance,  shown  for  cavity 

in  horizontal  (subscript  h)  beam:  Ih  indicates  cavity  length;  dhx,N  and 
dhy,N  indicate  hole  sizes  in  mirror  region;  dhx,i  and  dhy,i  indicate  hole 
sizes  in  hrst  taper  period;  ah^N  indicates  periodicity  in  mirror  region; 
a/l  l  indicates  periodicity  in  hrst  taper  period,  Wh  indicates  beam  width. 

The  corresponding  parameters  are  similarly  introduced  for  the  vertical 
beam  (with  subscript  v).  The  thickness  of  both  beams  (in  the  z  direc¬ 
tion)  is  t.  Parameters  are  changed  linearly  inside  the  taper,  (c)  Field 
pattern  of  Ey  for  cavity  mode  localized  by  horizontal  beam.  Param¬ 
eters  are:  a/i,Ar=453  nm,  a„,Ar=272  nm,  dhx,i/dhx,N  =  dhy,i/dhy,N=0.5, 
^h^l/^h^N  ^v,l/^v,N  0.7,  Ifi/^h^N  ^v/^v,N  0.83,  Wfij 1.65, 

Wxj/Cty^N  1.8,  djiy^JSl  /  W}i  dyx^N  /  0.7,  djix^N  /  dyy  f\J  !  0.5, 

refractive  index  n  =  3.37,  with  slab  thickness  t/ah,N=^-^5,  N=5,  and 
6  mirror  periods  for  both  beams.  Resonant  wavelength  is  1.55  /rm  with 
Q=19,000  and  V=0.35(A/n)^.  (d)  Field  pattern  of  Ex  for  cavity  local¬ 
ized  by  vertical  beam,  n  =  3.46  and  other  parameters  same  as  in  (c). 
Resonant  wavelength  is  1103  nm  with  Q=1900  and  V=0.47(A/n)^.  .  70 

7.3  (a)  Change  in  cavity  resonant  wavelength  and  Q  in  crossed  beam  struc¬ 

ture  for  cavity  mode  localized  in  horizontal  beam  as  cavity  length  is 
varied.  Maximum  Q  occurs  for  1^/ ah=1.2b.  (b)  Change  in  horizontal 
(long  wavelength)  cavity  resonant  wavelength  and  Q  in  crossed  beam 
structure  as  Wh  is  varied,  (c)  Change  in  vertical  beam  quality  fac¬ 
tor  and  frequency  conversion  hgure  of  merit  EOM  =  Q1Q2/ \/ViV2  as 
a  function  of  Wh-  (d)  Change  in  in-plane  quality  factor  in  x  direction 
(i.e.  radiated  power  collected  at  x  and  -x  edges  of  simulation  space)  for 
resonance  localized  by  vertical  beam  as  a  function  of  Wh  for  structure 
shown  in  Fig.  7.2  (c)/(d) .  71 
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7.4  (a)  Dielectric  constant  and  (b)  FDTD-simulated  electric  field  Ey  for 

nanobeam  cavity  crossed  with  orthogonal  waveguide  for  =  0.75tc/i. 
(c)/(d)  and  (e)/(f)  show  the  same  for  increasing  width  of  the  orthog¬ 
onal  waveguide,  illustrating  cavity  confinement  weakens  as  orthogo¬ 


nal  waveguide  width  is  increased.  For  (c)/(d),  =  tCh,  for  (e)/(f) 

=  1.67tc/j .  72 


7.5  (a)  Crossed  nanobeam  cavity  design,  showing  intersecting  orthogonal 

nanobeams  with  taper  and  mirror  regions,  as  well  as  central  cavity. 

(b)  Detail  of  white  box  in  (a).  Parameters  used  to  form  resonance, 
shown  for  cavity  in  horizontal  (subscript  h)  beam:  indicates  cavity 

length;  dhx,N  and  dhy,N  indicate  hole  sizes  in  mirror  region;  dhx,i  and 
dhy,i  indicate  hole  sizes  in  first  taper  period;  ah,N  indicates  periodic¬ 
ity  in  mirror  region;  indicates  periodicity  in  first  taper  period,  Wh 
indicates  beam  width.  The  corresponding  parameters  are  similarly  in¬ 
troduced  for  the  vertical  beam  (with  subscript  n).  The  thickness  of 
both  beams  (in  the  z  direction)  is  t.  Parameters  are  changed  linearly 
inside  the  taper,  (c)  3D  FDTD  simulation  of  field  pattern  of  Ey  for 
cavity  localized  in  horizontal  beam  by  tapering  hole  dimensions  and  lat¬ 
tice  constant  in  central  region.  Parameters  are:  ah,N  =  at;,Af=449  nm, 
dhx,l/dhx,N  dhy  i/dhy^jsi  0.5,  Clh,l/(^h,N  0.7,  lh/(^h,N 

lvl^v,N  1*^5  ^h/^h,N  "^v/^VyN  1.65,  dfiy^jsi j Wfi  dyx,N /^v  0.7, 

dhx,N/cih  =  dyy^p^/ax,N=^-5,  refractive  index  n  =  3.37,  with  slab  thick¬ 
ness  t/ah,N=^-^5,  N=8,  and  6  mirror  periods  for  both  beams.  Res¬ 
onant  wavelength  is  1.55  /xm  with  Q=12,000  and  V=0.44(A/n)^.  (d) 

Field  pattern  of  Ex  for  cavity  localized  in  vertical  beam.  Parameters 
are  same  as  in  (c) .  73 
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7.6  (a)  Scanning  electron  microscope  image  of  crossbeam  structures  with 
identical  parameters  in  both  beams.  Structures  are  fabricated  by  e- 
beam  lithography,  dry  etching,  and  wet  etching,  (b)  Experimental 
setup  for  cross-polarized  reflectivity  measurements  to  characterize  cav¬ 
ity  resonances.  PBS  indicates  polarizing  beamsplitter.  Cavity  polar¬ 
ization  is  oriented  45  degrees  {\H  +  V))  from  orthogonal  input  (|E)) 
and  measurement  {\H))  polarizations,  (c)  Cross-polarized  reflectivity 
measurement  of  structure  in  (a),  showing  two  resonances  at  1571.2  nm 
(Q=4700)  and  1573.9  nm  (Q=7200).  Solid  line  indicates  fit  to  sum  of 

two  Fano  lineshapes .  75 

7.7  (a)  Second  harmonic  characterization  of  structure  with  two  resonances 

nearly  degenerate  in  frequency  as  a  function  of  laser  wavelength  for  3 
polarizations.  Two  modes  with  orthogonal  polarization  are  visible,  (b) 
Second  harmonic  intensity  as  a  function  of  incident  laser  polarization. 
Vertical  axis  indicates  wavelength  of  laser;  horizontal  axis  indicates 
angle  of  polarization.  Color  indicates  second  harmonic  intensity.  Dot¬ 
ted  horizontal  lines  indicate  traces  in  (c).  (c)  Line  plots  of  second 

harmonic  generation  measured  at  different  polarizations  for  three  laser 
wavelengths  shown  in  (b).  Red  lines  indicate  hts  for  two  cavity  modes 
with  polarizations  separated  by  exactly  90  degrees .  76 

7.8  (a)  Schematic  of  sum  frequency  generation.  Light  from  two  CW  lasers 
is  coupled  into  two  cavity  resonances  at  1552.8  nm  and  1558.9  nm. 
Nonlinear  frequency  conversion  produces  light  at  the  second  harmonic 
frequencies  of  each  laser,  as  well  as  at  the  sum  frequency,  777.9  nm.  (b) 
Sum-frequency  generation  from  structure  with  resonances  as  indicated 

in  (a) .  77 
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7.9  (a)  SEM  image  of  cross-beam  structures  fabricated  in  164  nm  thick 

GaAs  membrane  with  a/i_Ar=470  nm  and  a„^Ar=320  nm  fabricated  in 
164  nm  membrane,  (b)  Thresholded  binary  image  of  SEM  used  for 
simulating  fabricated  structure,  (c)  FDTD  simulation  of  Ey  for  cavity 
resonance  at  1477  nm.  (d)  Simulated  for  cavity  resonance  at  1043 
nm.  (e)  Reflectivity  measurement  of  cavity  mode  at  1482.7  nm  with 
Q=6600.  (f)  Reflectivity  measurement  of  cavity  resonance  at  1101  nm 


with  Q=1000 .  80 

7.10  Cross-polarized  reflectivity  of  crossbeam  structure  with  resonances  sep¬ 
arated  by  523  nm.  (a)  Resonance  at  1546.6  nm  (Q=1600).  (b)  Reso¬ 
nance  at  1023  nm  (Q=500) .  81 


7.11  Concept  illustration  of  doubly  resonant  second  harmonic  generation  in 
photonic  crystal  crossbeam  nanocavity.  Incident  light  (red)  is  coupled 
into  the  structure  via  a  grating,  transmitted  to  the  cavity,  frequency 
converted,  and  outcoupled  through  a  separate  grating  (green) .  81 

8.1  (a)  SEM  image  of  a  fabricated  photonic  crystal  cavity  in  CaP.  Scale 

bar  indicates  200  nm.  (b)  FDTD  simulation  of  electric  held  intensity 
of  the  fundamental  cavity  mode.  The  mode  is  primarily  y-polarized  (c) 
Schematic  illustrating  coupling  of  molecule  to  cavity,  (i)  DNQDI/PMMA 
is  deposited  over  the  entire  structure,  (ii)  DNQDI/PMMA  is  litho¬ 
graphically  dehned  on  cavity  region,  (d)  Bulk  photoluminescence  spec¬ 
trum  of  DNQDI  when  excited  with  a  633  nm  HeNe  laser.  The  molecule 
has  a  peak  in  its  absorption  at  this  excitation  wavelength,  (e)  Chemical 
structure  of  DNQDI  molecule .  84 
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8.2  (a)  Cross-polarized  reflectivity  measurement  of  a  cavity.  Box  indicates 
fundamental  cavity  mode,  (b)  Reflectivity  spectrum  of  high  quality 
factor  fundamental  cavity  mode  (box  in  (a)).  Spectrum  shows  addi¬ 
tional  peaks  at  shorter  wavelengths  from  higher  order  cavity  modes. 

Solid  line  shows  Lorentzian  fit  with  quality  factor  10,000.  (c)  Photo¬ 
luminescence  collected  from  the  same  photonic  crystal  cavity  in  (a-b) 
after  molecules  are  deposited  on  cavity,  x-polarized  emission  is  shown 
in  blue;  y-polarized  emission  is  shown  in  red.  Inset:  PL  measurement 
of  fundamental  cavity  mode  (black  box).  Line  indicates  Lorentzian  fit 
with  Q=10,000.  (d)  Quality  factors  measured  from  reflectivity  before 
molecule  deposition  and  photoluminescence  after  molecule  deposition 
from  the  high  Q  cavity  mode  for  structures  with  lattice  constant  a  and 
hole  radius  r/a  tuned  so  that  the  fundamental  cavity  resonance  shifts 
across  the  photoluminescence  spectrum  of  the  molecule.  Blue  open  cir¬ 
cles  indicate  reflectivity  measurements  for  the  cavities  that  were  also 
measured  in  PL  (blue  closed  circles) .  87 

8.3  (a)  Scanning  confocal  image  of  photoluminescence  from  DNQDI  doped 
PMMA  float-coated  onto  a  photonic  crystal  membrane.  Pixel  size  is 
200  nm  and  scale  bar  indicates  2  pm.  (b)  Scanning  confocal  image 
of  DNQDI  PL  after  electron  beam  lithography  is  used  to  remove  all 
molecules,  except  for  the  ones  coating  the  cavity  region  at  the  center. 

The  same  imaging  laser  power  as  in  (a)  was  used.  Pixel  size  is  80  nm 
and  scale  bar  indicates  2  pm.  (c)  PL  spectrum  from  the  fundamental 
mode  of  photonic  crystal  cavity  after  selective  removal  of  molecules 
by  e-beam  lithography,  (d)  Atomic  force  microscopy  image  showing 
localization  of  DNQDI-doped  PMMA  to  the  cavity  region.  PMMA 
thickness  is  12  nm.  Scale  bar  indicates  500  nm .  88 
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8.4  (a)  DNQDI  photoluminescence  as  a  function  of  temperature.  For  tem¬ 

peratures  below  50K,  narrow  lines  appear,  (b)  Individual  spectral  fea¬ 
tures  qualitatively  match  the  expected  emission  from  a  zero  phonon 
line  and  phonon  wing.  As  much  as  30%  of  molecule  emission  is  mea¬ 
sured  to  come  from  the  narrow  spectral  line .  89 

8.5  Spin-grown  thin-hlm  crystals  of  DNQDI  molecules  on  different  sub¬ 
strates  formed  by  using  different  spin  speeds,  as  indicated.  Scale  bars 

all  indicate  10  pm .  91 

8.6  (a)  Optical  microscope  image  of  a  crystal  formed  by  spinning  at  3000 

rpm  on  a  glass  substrate,  (b)  AFM  image  of  crystal  spun  at  2500  rpm 
on  glass,  indicating  crystal  height  of  few  hundred  nm  and  lateral  size 
of  ~1  pm.  (c)  Photoluminescence  from  a  spun  crystal  excited  by  a 
633  nm  HeNe  laser,  (d)  Confocal  photoluminescence  scan  of  region  of 
crystal  in  (c)  measured  by  a  Si  single  photon  counter,  c-axis  indicates 
photoluminescence  counts .  92 

8.7  (a)  QD  PL  from  an  unprocessed  region  of  sample  as  temperature  is 

varied  between  25K  and  300K.  The  pump  is  a  405  nm  CW  laser  diode 
at  700  pW.  (b)  Photoluminescence  as  a  function  of  pump  power  at  lOK 
showing  broadening  at  higher  energies  with  increasing  power,  indicative 
of  the  presence  of  excited  states,  (c)  Integrated  counts  of  low  energy 
half  of  PL  spectrum  (to  minimize  contribution  of  excited  states,  which 
varies  with  temperature).  Intensity  is  decreased  by  a  factor  of  4  at 
300K.  (d)  Semilog  plot  of  integrated  counts  of  low  energy  half  of  PL 
spectrum  versus  inverse  temperature.  Fit  to  exponential  (red  line) 
gives  activation  energy  77^=161  meV .  96 
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8.8  (a)  Experimental  setup  for  time-resolved  measurements.  3  ps  pulses  at 
800  nm  from  Ti; Sapphire  laser  are  frequency  doubled  in  a  BiBO  crystal. 

A  grating  is  used  to  filter  the  400  nm  light,  which  is  passed  through 
a  dichroic  mirror  onto  the  sample.  The  photoluminescence  emitted  by 
the  quantum  dots  is  transmitted  through  the  dichroic,  passed  through 
a  long-pass  hlter  (LPF)  to  remove  any  residual  pump,  and  sent  to  a 
spectrometer  or  streak  camera  (for  time-resolved  measurements),  (b) 
Time-resolved  streak  camera  measurements  showing  lifetime  integrated 
across  low  energy  half  of  PL  spectrum  for  different  temperatures.  Inset: 
excitation  pulse  at  400  nm.  Red  line  indicates  ht  with  decay  time  12  ps, 
limited  by  instrument  resolution,  (c)  Exponential  hts  of  time-resolved 
data  at  250K,  showing  initial  and  hnal  decay  times  r*  and  r/.  (d) 

Measured  lifetimes  as  a  function  of  sample  temperature .  97 

8.9  (a)  SEM  image  of  photonic  crystal  nanocavity,  (b)  PL  measured  at  12K 

with  405  nm  CW  pump  from  unprocessed  region  of  thinner  sample 
used  for  photonic  crystal  measurements.  QD  wavelength  is  slightly 
blueshifted  from  Fig.  1.  (c)  PL  measurement  indicating  emission  of 
quantum  dots  coupled  into  cavity  modes.  Fundamental  cavity  mode 
is  indicated  by  black  circle,  (d)  Lorenztian  £t  of  fundamental  cavity 
mode  at  681.7  nm  with  Q=2800.  (e)  FDTD-simulated  electric  held 
intensity  for  fundamental  cavity  mode .  98 

8.10  (a)  PL  from  a  different  photonic  crystal  cavity  structure  with  higher 
order  cavity  mode  aligned  to  wavelength  of  strongest  QD  PL  measured 
with  405  nm  CW  pump.  Black  circle  indicates  mode  of  interest.  Inset: 
FDTD-simulated  electric  held  intenstiy  for  higher  order  cavity  mode 
indicated  by  circle  in  (a),  (b)  High  resolution  spectra  of  QD  PL  from 
cavity  in  (a).  Measurements  are  taken  at  lOK.  The  two  single  lines  at 
low  power  are  indicative  of  single  quantum  dots,  (c)  Change  in  wave¬ 
length  of  quantum  dots  (measured  at  500  nW)  and  cavity  (measured 
at  100  /iW)  as  a  function  of  temperature.  Solid  lines  indicate  quadratic 

hts.  (d)  Change  in  QD  linewidth  as  a  function  of  temperature .  100 
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Chapter  1 
Introduction 


1.1  Motivation 

Nonlinear  optical  frequency  conversion  is  the  basis  of  numerous  photonic  devices  from 
green  laser  pointers  to  optical  parametric  oscillators.  Since  the  hrst  demonstration 
of  second  harmonic  generation  from  a  ruby  crystal  in  1961  [3],  on  the  heels  of  the 
development  of  the  laser [4],  the  conventional  view  of  nonlinear  frequency  conversion 
has  involved  high  optical  pump  powers  (often  Watts),  exotic  materials  such  KDP  and 
BBO  that  are  disconnected  from  semiconductor  electronics,  and  large  crystals  with 
dimensions  of  the  order  of  millimeters  to  centimeters. 

As  the  held  of  nanophotonics  has  emerged  over  the  past  decade,  the  boundary 
between  semiconductor  electronics  and  photonics  has  vastly  diminished.  In  nonlinear 
optics,  however,  the  divide  remains.  This  thesis  concerns  the  questions  of  how  to 
perform  nonlinear  frequency  conversion  processes  with  materials,  optical  power  lev¬ 
els,  and  device  sizes  amenable  to  on-chip  semiconductor  integration.  Such  integrated 
nonlinear  optical  frequency  conversion  devices  could  lead  to  miniaturized  devices  for 
applications  such  as  on-chip  spectroscopy  and  sensing  of  ultrasmall  volumes  of  ana¬ 
lyte  from  visible  wavelengths  to  mid-infrared.  Furthermore,  emerging  on-chip  optical 
routing  applications  requiring  nonlinear  functionality  are  fundamentally  incompatible 
with  conventional  nonlinear  optics.  Many  of  the  optical  elements  such  as  beamsplit¬ 
ters  and  detectors  that  can  be  used  for  quantum  optical  computing  have  now  been 
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demonstrated  on-chip  [5],  replacing  an  optical  table’s  worth  of  carefully  aligned  com¬ 
ponents.  One  element  that  remains  off-chip,  however,  has  been  a  source  of  entangled 
photons  coming  from  a  nonlinear  optical  process.  Classical  optical  information  pro¬ 
cessing  and  routing  could  also  beneht  from  on-chip  nonlinear  elements.  Four-wave 
mixing [6]  has  been  proposed  as  a  way  to  enable  frequency  multiplexed  on-chip  sources; 
additionally,  many  computational  and  signal  processing  applications  require  nonlinear 
elements. 

Dispersion  engineered  nanophotonic  structures  can  simultaneously  bring  nonlinear 
frequency  conversion  on-chip  and  also  robustly  satisfy  phase  matching  to  achieve 
efficient  devices.  This  thesis  concerns  the  design,  fabrication,  characterization,  and 
analysis  of  nanophotonic  structures  and  semiconductor  materials  that  can  meet  these 
requirements. 

1.2  Contributions 

The  contributions  in  this  dissertation  include: 

•  The  hrst  photonic  crystal  cavities  in  gallium  phosphide,  which  at  the  time  of 
demonstration  were  the  shortest  wavelength  photonic  crystal  cavities  fabricated 
in  a  high  (>3.0)  refractive  index  material 

•  The  demonstration  of  430%/W  (10“^  for  11  pW  coupled  power)  continuous  wave 
second  harmonic  generation  in  gallium  phosphide  photonic  crystal  cavities,  an 
8-orders  of  magnitude  improvement  on  previous  work  in  InP  [7] 

•  The  hrst  demonstration  of  continuous  wave  sum-frequency  generation  in  pho¬ 
tonic  crystal  cavities,  fabricated  in  GaP 

•  The  hrst  demonstration  of  second  harmonic  generation  in  photonic  crystal 
waveguides,  fabricated  in  GaP 

•  The  demonstration  of  the  fastest  optically  pumped  single  photon  source,  based 
on  a  quantum  dot  pumped  by  intracavity  second  harmonic  generation  from  a 
telecom-wavelength  laser 
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•  The  design,  fabrication,  and  characterization  of  multiply  resonant  photonic  crys¬ 
tal  cavities  with  record  frequency  separation  between  resonances  (more  than  600 
nm  resonant  wavelength  separation  measured  experimentally) 

•  The  hrst  demonstration  of  lithographic  techniques  to  align  thin  dims  of  molecules 
to  high-Q  photonic  crystal  cavities 

•  Spectral  and  temporal  characterization  of  InGaAs/GaP  quantum  dot  lumines¬ 
cence,  including  the  first  observation  of  narrow  spectral  lines  indicative  of  single 
quantum  dots  via  the  enhanced  outcoupling  of  a  photonic  crystal  cavity 


1.3  Dissertation  outline 

Ghapters  2  and  3  contain  introductory  material,  with  Ghapter  2  providing  background 
on  second  order  nonlinear  optical  frequency  conversion  and  Ghapter  3  discussing  thin 
membrane  two  dimensional  photonic  crystals. 

Ghapter  4  describes  the  design,  fabrication,  and  characterization  of  the  first  GaP 
photonic  crystal  cavities,  a  crucial  ingredient  for  the  results  in  much  of  the  later 
chapters  of  this  dissertation.  These  results  were  published  in  Ref.  [8],  with  Andrei 
Faraon  as  a  collaborator. 

Ghapter  5  presents  three  experiments  demonstrating  three  wave  mixing  in  trian¬ 
gular  lattice  photonic  crystal  structures  including  second  harmonic  generation  and 
sum  frequency  generation  in  linear  three  hole  defect  cavities,  as  well  as  second  har¬ 
monic  generation  in  W1  photonic  crystal  waveguides.  These  results  were  published  in 
Refs.  [9,  10,  11].  Ziliang  Lin  collaborated  on  cavity  experiments  and  Sonia  Buckley 
on  waveguide  experiments;  the  material  for  all  three  experiments  was  grown  by  Fariba 
Hatami  with  support  from  W.  Ted  Masselink. 

Ghapter  6  uses  the  frequency  conversion  techniques  demonstrated  in  Ghap.  5  to 
realize  a  fast  (few  hundred  MHz)  quantum  dot  single  photon  source  triggered  using 
intracavity  second  harmonic  generation  from  a  telecom-wavelength  source.  This  work 
was  a  collaboration  with  Sonia  Buckley  and  Arka  Majumdar,  with  material  provided 
by  Hyochul  Kim  and  Pierre  Petroff  and  was  published  in  Ref.  [12]. 
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Chapter  7  describes  the  design,  fabrication  and  characterization  of  multiply  reso¬ 
nant  orthogonally  crossed  beam  photonic  crystal  structures  This  work  was  a  collab¬ 
oration  with  Sonia  Buckley  and  was  published  in  Refs.  [13,  14], 

Chapter  8  presents  specotroscopic  measurements  of  two  systems  of  visible  quan¬ 
tum  emitters:  DNQDI  molecules  embedded  in  thin  PMMA  dims  and  InGaAs/GaP 
quantum  dots,  coupled  to  GaP  photonic  crystal  cavities.  The  work  on  molecules  was 
performed  jointly  with  Anika  Kinkhabwala  in  collaboration  with  W.E.  Moerner,  with 
molecules  provided  by  Yuri  Avlasevich  and  Klaus  Muellen  and  GaP  material  grown 
by  Fariba  Hatami  with  support  from  W.  Ted  Masselink;  the  principal  results  were 
published  in  Ref.  [15].  Material  for  the  InGaAs/GaP  quantum  dot  spectroscopy, 
published  in  Ref.  [16],  was  grown  by  Yuncheng  Song  with  support  from  Minjoo  Larry 
Lee. 

Finally,  Ghapter  9  concludes  the  thesis  and  provides  an  outlook  for  future  work 
in  this  held. 

All  work  presented  in  this  dissertation  was  performed  under  the  supervision  of 
Jelena  Vuckovic. 


Chapter  2 

Nonlinear  optical  frequency 
conversion 


The  goal  of  frequency  conversion  is  to  change  the  frequency  of  the  signal;  at  optical 
frequencies,  this  corresponds  to  changing  the  color  of  light.  Such  processes  hnd  appli¬ 
cations  in  a  diverse  range  of  helds,  including  sensing,  spectroscopy,  metrology,  commu¬ 
nications,  and  quantum  information.  More  specihcally,  the  uses  of  optical  frequency 
conversion  include  generation  of  light  sources  at  frequencies  otherwise  inaccessible 
with  available  lasers,  upconversion  and  downconversion  of  a  hxed  frequency  signal 
(e.g.  from  an  atom  or  semiconductor  light  emitter)  to  a  more  desirable  frequency, 
and  generation  of  entangled  photon  pairs  via  spontaneous  parametric  downconver¬ 
sion.  This  chapter  introduces  the  basic  physics  for  these  processes;  further  details  can 
be  fonnd  in  many  references  including  [17,  18,  19],  which  this  chapter  follows  in  part. 

2.1  Classical  description  of  atomic  polarization 

Because  the  emission  of  photons  at  optical  freqnencies  is  based  on  electrons  in  atoms, 
which  can  be  approximated  classically  as  linear  simple  harmonic  oscillators  with 
qnadratic  potentials,  optical  freqnency  conversion  cannot  occur  if  the  incident  light  is 
monochromatic  and  continnous  wave  (time-independent).  The  motion  of  an  atomic 
electron  of  mass  m  and  charge  q  with  phenomenological  damping  constant  7  and 
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spring  constant  K  driven  by  an  electromagnetic  field  E  =  Eocos(a;t)  can  be  modeled 
by  the  simple  linear  differential  equation: 

^ +  +a;^r  =  — Eocos(a;t)  (2.1) 

dP  dt  m 

where  (5  =  '^ /m  and  ujq  =  Kjm.  This  equation  has  steady-state  solution 

r(t)  =  —  =Eocos(a;t  -|-  S)  (2.2) 

y/(a;2  -  a;^)^  + 

2  2 

where  tan  (5)  =  — producing  a  linear  polarization  p  =  gr  =  aE  where  a  is 
the  polarizability.  The  total  polarization  will  then  be  P  =  NaE  where  N  is  the 
atom  density  (number  of  atoms  per  unit  volume).  From  this,  it  is  clear  that  a  linear 
combination  of  sinusoidal  input  signals  (into  which  any  other  signal  can  be  decom¬ 
posed,  via  Fourier  analysis)  produces  a  linear  superposition  of  sinusoidal  outputs  at 
the  same  frequencies.  Accordingly,  to  perform  frequency  conversion  requires  either  a 
time- dependent  perturbation  to  the  atom,  or  that  we  change  our  model  to  include  a 
nonlinearity.  In  fact,  a  more  accurate  model  does  include  such  a  nonlinearity  due  to 
the  anharmonicity  of  the  potential  for  the  atomic  electron  at  large  values  of  electric 
held.  The  held  at  which  this  anharmonicity  becomes  significant  can  be  estimated  by 
assuming  such  a  held  would  produce  a  displacement  of  the  electron  in  the  atom  on 
the  scale  of  the  Bohr  radius  Oq  =  dTreoh^/me^,  where  Cq  is  the  permittivity  of  free 
space,  h  is  the  reduced  Planck’s  constant,  m  is  the  mass  of  the  electron,  and  — e  is 
the  charge  on  the  electron.  This  results  in  an  atomic  electric  held 
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dvreoao 


(2.3) 


giving  E  =  5.14  x  10^^  V/m  . 

This  anharmonicity  adds  quadratic  and  higher  terms  in  r  to  Eqn.  2.2  (or  cubic  and 
higher  terms  in  the  potential),  producing  terms  with  new  frequencies  in  the  solution, 
starting  with  a  frequency  of  2u  for  the  lowest  order  nonlinearity.  This  results  in  the 
addition  of  nonlinear  terms  to  the  polarization. 


CHAPTER  2.  NONLINEAR  OPTICAL  PREQUENCY  CONVERSION 


7 


p  =  aE  +  q;2E^  +  ttsE^... 


(2.4) 


For  a  crystal  with  a  centrosymmetric  structure,  however,  only  higher  order  terms 
with  even  powers  of  r  appear  in  the  potential  (or  odd  powers  of  r  in  the  restoring 
force),  that  is: 


Pcentrosymmetric  oE  -|-  Q;3E  ... 


(2.5) 


2.2  Maxwell’s  equations-based  description  of  non¬ 
linear  polarization 

The  nonlinear  polarizaton  can  also  be  described  macroscopically  by  Maxwell’s  Equa¬ 
tions.  For  propagation  of  electromagnetic  waves  in  a  medium  with  no  free  charges  or 
currents,  the  key  equations  are: 


(9B 

V  X  E  = - —  (2.6a) 

dt  ^  ^ 

(9D 

V  X  H  =  —  (2.6b) 

dt 

where  B  =  poH  for  a  nonmagnetic  material  where  po  is  the  permeability  of  free 
space,  D  =  eoE  +  P  =  e^eE  -|-  P^^,  eo  is  the  permittivity  of  free  space  and  P^^  is  the 
nonlinear  polarization. 

This  leads  to  a  driven  wave  equation 


V^E 


e  d‘^E 

C  dC 


^2pNL 


(2.7) 


meaning  that  a  nonlinear  polarization  will  act  as  a  source  for  a  wave  with  the  fre¬ 
quency  contained  in  this  nonlinear  polarization.  Returning  to  Eqn.  2.4,  we  can  write 
the  macroscopic  nonlinear  polarization 


P  =  e„x'"’E  +  e„xP>E"  +  ... 


(2.8) 


CHAPTER  2.  NONLINEAR  OPTICAL  PREQUENCY  CONVERSION 


where  1  +  =  e.  For  a  monochromatic  scalar  input  field  E{t)  =  EQCOs{ut),  this 

leads  to 

P^^{t)  =  (2-9) 

Using  the  trigonometric  identity  cos(2a;)  =  this  then  results  in  a  term  oscil¬ 

lating  with  frequency  2u  (as  well  as  a  constant  term,  which  converts  oscillating  fields 
to  DC  and  is  known  as  optical  rectification),  which  can  act  as  a  source  in  the  wave 
equation  2.7.  As  this  produces  a  term  oscillating  at  twice  the  frequency  of  the  original 
light,  the  process  is  known  as  second  harmonic  generation.  Returning  to  Eqn.  2.3, 
this  indicates  displacements  of  the  scale  of  the  Bohr  radius  in  the  electron  leading  to 
optical  nonlinearities  could  be  produced  by  a  nonlinear  optical  coefficient  ~  2 
pm/V. 

Similarly,  for  two  incident  scalar  fields  with  different  frequencies  ui  and  UJ2-,  that 
is.  Exit)  =  EiflCos{uit)  and  7^2 (t)  =  E2fiCos{u2t),  this  produces  additional  terms  in 
the  nonlinear  polarization  with  frequencies  equal  to  the  sum  and  difference  of  the  two 
incident  frequencies,  Ui  +  U2  and  cui  —  U2'. 

^(2) 

2  +  El, 0^2*06'^"'”"^^*]  + 

2  [Ei,oEi,g  -|-  E2, 0^2,0]  }  -|-  c.c.  (2.10) 

where  c.c.  indicates  the  complex  conjugate,  the  terms  on  the  first  line  indicate  second 
harmonic,  terms  on  the  second  line  indicate  sum  and  difference  frequency,  and  terms 
on  the  third  line  indicate  optical  rectification. 

More  generally,  the  second  order  suspectibility  Xijk  takes  the  form  of  a  tensor, 
in  this  case  third-rank,  resulting  in  a  27-element  matrix  that  depends  on  the  polar¬ 
izations  of  all  three  waves  relative  to  the  crystal  structure.  The  number  of  elements 
in  this  tensor  can  be  reduced  to  18  by  intrinsic  permutation  symmetry  (arbitrary 
swapping  of  j  and  k),  allowing  the  use  of  a  contracted  notation,  resulting  in; 
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Pl{ui  +  LJ2) 
P2{u)i  +  U2) 
-^3(1^1  +  1^2) 


dll 

di2 

di3 

di4 

di5 

die 

4eo 

d2i 

d22 

d23 

d24 

d25 

d26 

dsi 

d32 

<^33 

dzi 

dzz 

<^36 

Ex{uJi)Ex{uJ2) 

Ey{u)i)Ey{uJ2) 

Ez{uJi)Ez{uJ2) 

Ey{uJi)Ez{uJ2)  +  Ey{oJ2)Ez{oJi) 
Ex{oJi)Ez{oJ2)  +  Ex{uJ2)Ez{uJi) 
Ex{oJi)Ey{oJ2)  +  Ex{i^2)EY{(^l) 


For  second  harmonic  generation,  The  number  of  independent  el¬ 

ements  in  the  tensor  can  then  be  reduced  to  10  unique  elements  with  all  3  indices 
permuted  (full  permutation  symmetry)  when  the  process  occurs  far  away  from  mate¬ 
rial  resonances.  The  crystal  structure  and  symmetries  determines  which  of  these  10 
elements  are  non-zero. 

Under  certain  approximations,  such  as  the  slowly  varying  envelope  approximation 
(which  allows  the  order  of  the  differential  equation  to  be  reduced),  this  equation  can 
be  converted  into  coupled  wave  equations  for  each  frequency  for  plane  waves,  focused 
beams,  and  several  other  common  geometries.  In  certain  regimes  of  operation,  usually 
where  the  wave  at  one  frequency  is  undepleted,  simple  solutions  to  these  coupled  wave 
equations  exist.  For  example,  for  sum  frequency  generation  with  =  ui  +  U2,  for 
waves  propagating  a  distance  L: 


h  =  (AkLI2)  (2.12) 

nin2n^eQC 

where  the  deff  is  the  nonlinearity  determined  by  Eqn.  2.11,  n  is  the  refractive  index. 
This  relationship  is  illustrated  in  Figure  2.1. 
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A  k  L/(27i) 


Figure  2.1:  Reduction  in  sum-frequency  signal  over  a  distance  L  due  to  wavevector 
mismatch  (as  described  in  Eqn.  2.12). 

2.3  Phase  matching  and  materials  for  efficient  fre¬ 
quency  conversion 

In  order  to  generate  new  frequencies  efficiently,  two  important  constraints  must  be 
satisfied:  (1)  energy  conservation  and  (2)  momentum  conservation,  better  known  in 
this  context  as  phase  matching.  Energy  conservation  stipulates  that  the  frequency  at 
which  the  nonlinear  polarization  is  generated  is  the  frequency  radiated.  Momentum 
conservation,  or  phase  matching,  means  that  to  continue  to  add  up  coherently  over 
long  distances,  the  different  frequencies  must  propagate  at  the  same  phase  velocity; 
otherwise,  the  efficiency  will  be  reduced.  The  phase  velocity  mismatch  Ak  can  be 
written  as  Ak  =  —  k2  —  ki  for  sum  frequency  or  difference  frequency  generation 

and  Ak  =  k2  —  2ki  for  second  harmonic  generation.  For  second  harmonic  generation 
with  monochromatic  plane  waves,  Ak  =  ^{n2uj  —  meaning  that  the  indices  of 
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refraction  at  the  two  frequencies  should  be  identical  to  minimize  phase  mismatch. 
In  most  materials,  however,  normal  dispersion  leads  to  a  difference  in  the  refractive 
indices  at  different  frequencies. 

In  bulk  materials,  birefringence  in  materials  such  as  LiNbOs  and  KDP  can  be 
used  to  satisfy  phase  matching,  as  beams  at  different  frequencies  can  propagate  along 
different  crystal  directions  with  different  refractive  indices.  In  practice,  this  requires 
careful  control  over  the  angle  and  frequently  the  temperature  as  well,  as  refractive 
indices  are  temperature  dependent.  For  a  birefringent  material,  the  linear  dielectric 
constant  epsilon  can  be  represented  in  a  tensor  form: 


e  = 


exx  0  0 

0  eyy  0 
0  0  ezz 


(2,13) 


where  exx,  eyy,  and  ezz  can  have  different  values,  depending  on  the  crystal  structure, 
allowing  phase  matching  to  be  satished. 

When  considering  which  materials  to  use  for  nonlinear  optical  frequency  conver¬ 
sion,  conventionally,  large  nonlinear  optical  coefficient,  transparency  at  the  frequen¬ 
cies  of  interest,  and  birefringence  are  the  key  parameters,  with  other  properties  such 
as  high  refractive  index  and  large  thermal  conductivity  also  desirable.  This  typi¬ 
cally  results  in  a  relatively  exotic  choice  of  materials  for  nonlinear  optical  frequency 
conversion  that  optimize  these  quantities  of  interest,  but  are  difficult  to  integrate 
with  semiconductor  processing  and  electronics.  These  relevant  properties  for  semi¬ 
conductor  materials,  both  Si  and  two  III-V  semiconductors,  GaAs  and  GaP,  and  their 
comparison  to  the  frequently  used  nonlinear  optical  material  LiNbOs  are  shown  below 
in  Table  2.3. 

We  can  observe  that  Si  does  not  have  a  second  order  optical  nonlinearity,  due  to 
the  centrosymmetric  nature  of  the  crystal,  making  it  a  poor  choice  for  second-order 
nonlinear  optical  processes.  GaAs  and  GaP  have  large  optical  nonlinearities,  broad 
transparency  ranges  (with  GaP  extending  into  the  visible  and  GaAs  far  into  the 
infrared),  large  thermal  conductivity  for  reducing  heat  buildup,  and  large  refractive 
index  useful  for  conhning  light  by  total  internal  reflection.  However,  the  zincblende 
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Si 

GaAs 

GaP 

LiNbOs 

dij  (at  2  yum) 

~0 

107 

45 

23 

Transparency  (yum) 

1.1-5 

0.9-17 

0.5-11 

0.4-4.5 

Refractive  index  (at  1.5  yum) 

3.5 

3.4 

3.1 

2.1-2.2 

Thermal  conductivity  (W/m  K) 

150 

50 

110 

5 

Birefringent 

No 

No 

No 

Yes 

Table  2.1:  Relevant  parameters  for  nonlinear  optical  frequency  conversion  for  semi¬ 
conductors  Si,  GaAs,  and  GaP  and  traditional  nonlinear  optical  material  lithium 
niobate.  d  =  di4  for  GaP,  GaAs,  d  =  ^33  for  LiNbOs.  Low  frequency  side  of  absorp¬ 
tion  region  determined  by  onset  of  two  phonon  absorption  in  the  lattice.  Values  from 
Refs.[l,  2]. 


lattice  of  GaAs  and  GaP  leads  to  an  isotropic  refractive  index,  meaning  that  phase 
matching  cannot  be  accomplished  via  birefringence.  Instead,  phase  matching  can  be 
achieved  by  structuring  the  material. 

2.4  Phase  matching  in  III-V  semiconductors 

The  hrst  methods  introduced  to  support  phase  matching  in  III-V  structures  involved 
structuring  the  material  on  a  macroscopic  scale  to  periodically  reverse  the  sign  of  the 
nonlinear  optical  coefficient  every  coherence  length.  This  could  be  accomplished  by 
gluing  chunks  of  the  material  together  (stack  of  plates  approach),  or  more  recently, 
by  etching  and  regrowth  of  the  material  (known  as  orientation  patterning[20,  21,  22]). 
Such  methods  are  known  as  “quasi  phase  matching”  and  allow  a  continued  increase 
in  nonlinearly  generated  signal  as  the  light  propagates  along  a  crystal. 

An  alternative  approach  is  to  structure  the  material  on  a  smaller  scale  to  change 
the  dispersion  relation,  and  accordingly  the  Ak.  For  example,  for  a  second  harmonic 
generation  process,  perfect  phase  matching  could  be  accomplished  with  an  arbitrary 
dispersion  relation  containing  the  two  {k,u)  points  {ki,ui)  and  (2A;i,  2c<;i). 

Several  possible  methods  exist  to  engineer  dispersion.  A  grating  with  periodicity 
a  allows  momentum  to  be  converted  in  integer  multiples  of  the  reciprocal  grating 
vectors  27i/a-,  accordingly  the  maximum  possible  phase  mismatch  in  such  a  structure 
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would  be  set  by  the  size  of  the  first  Brillouin  zone.  Nanoscale  waveguides  have  also 
been  designed  to  have  anomalous  dispersion  [6],  which  can  be  engineered  to  cancel 
normal  dispersion  at  a  particular  wavelength  for  a  particular  waveguide  geometry. 

Finally,  it  is  also  helpful  to  think  of  phase  matching  more  generally  as  an  overlap 
integral  between  the  electric  helds  at  different  frequencies:  in  the  case  of  planes 
waves,  this  leads  to  a  product  of  recovering  the  familiar  Ak  relation.  For  a 

nanoscale  cavity  structure,  however,  that  exists  in  a  region  smaller  than  the  coherence 
length  and  is  not  dehned  by  a  single  k  vector;  this  reduces  to  a  spatial  overlap  of  the 
electric  helds,  e.g.  assuming  all  helds  real,  the  overlap  integral  (3  for  the  three  helds 
can  be  dehned 


p  ^ _ Inl  dVEiE2E3 _ 

^ J  dVe\Ei\^^ J  J  dVe\Es\^ 

so  the  electric  helds  are  normalized  by  the  energy  and  the  numerator  integral  is  taken 
only  in  the  nonlinear  material. 

The  radiated  power  at  the  sum  or  diherence  frequency  generation  will  then  be 
proportional  to  A  nonlinear  coupled  mode  theory  approach  including  this 

overlap  can  be  then  used  to  model  the  frequency  conversion  process  if  the  spatial 
form  of  all  modes  is  known [23,  24,  25]. 

Note  that  such  a  held  overlap  also  falls  out  naturally  from  a  quantum  mechanical 
description  of  three  wave  mixing,  as  the  interaction  Hamiltonian  term  comes  from 
pNL.^;  procoss  is  thon  proportional  to  the  square  of  the  relevant  matrix 

elements,  again  resulting  in  a  dependence.  In  this  case,  the  total  Hamiltonian  can 
be  written: 


(2.14) 


H  =  hwittiai  +  huj2a2a2  +  hg  ^020]*“^  +  a^a^j  (2-15) 

where  the  second  term  describes  the  nonlinear  interaction  and 

a  =  j  'iyxfA)El(r)Ei(r)El(r)  (2.16) 

where  E  is  normalized  to  one  photon  per  mode.  In  the  limit  g  >  k  =  u/{2Q),  the 
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second  harmonic  and  first  harmonic  would  experience  Rabi  oscillation  [26]. 

In  order  to  a  define  a  quantity  normalized  to  1,  it  can  also  be  useful  to  define  a 
two-held  overlap  (this  can  be  normalized  as  powers  of  electric  held  are  the  same  in 
all  integrals  in  numerator  and  denominator;  this  type  of  expression  can  be  used  to 
describe  spatial  overlap  of  two  nanoscale  modes,  where  the  third  frequency  will  have 
a  classical  form).  In  this  case,  we  can  introduce  a  modihed  form  of  (3  as  7 [27],  which 
will  be  used  later  in  this  thesis. 

/nl  dVEiE2 

7  =  /  -  /  =  (2.17) 

Jj  dVe\E^\\  j  dVe\E2\^ 


Chapter  3 

Photonic  crystal  cavities 


Optical  components  such  as  high  reflectivity  mirrors  and  ultrahigh-Q  cavities  have 
formed  the  cornerstone  of  many  optical  devices  including  lasers  and  other  types  of 
oscillators.  Historically,  these  components  have  been  machined  on  the  macroscale. 
In  the  past  decade,  herculean  efforts  have  been  made  to  borrow  the  technology  of 
semiconductor  micro-  and  nanofabrication  to  dramatically  reduce  the  size  of  opti¬ 
cal  components  and  facilitate  their  integration  into  on-chip  devices.  Planar  pho¬ 
tonic  crystal  cavities,  the  focus  of  this  chapter,  are  an  example  of  this  technology, 
which  provides  a  means  to  form  high  reflectivity  micron-scale  mirrors  and  nanoscale 
high-Q  cavities.  Other  miniaturized  technologies  with  different  tradeoffs  in  quality 
of  confinement  compared  to  size  of  resonator  also  exist,  including  microdisks [28], 
microrings [29,  30],  microspheres [31],  etc.  Many  of  these  technologies  confine  light  in 
all  three  dimensions  by  total  internal  reflection  by  forming  the  nanophotonic  struc¬ 
ture  inside  a  higher  refractive  index  than  the  surrounding  medium.  Here,  we  discuss 
photonic  crystal  technologies,  which  use  distributed  Bragg  reflection  for  confinement 
in  at  least  one  dimension. 


3.1  Distributed  Bragg  reflection 

If  we  consider  the  interface  between  two  materials  with  different  refractive  indices,  we 
can  calculate  the  reflection  and  transmission  for  light  impinging  on  the  interface.  For 
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a  single  interface  at  normal  incidence,  this  leads  to  an  intensity  reflection  coefficient 
R  =  I  |2  fi  and  no  are  the  refractive  indices  of  the  two  materials.  If 

we  then  stack  a  set  of  these  interfaces,  we  must  consider  how  the  phases  of  the 
reflected  light  at  successive  surfaces  add  up.  In  general,  by  carefully  choosing  the 
thicknesses  and  refractive  indices  of  the  two  materials,  the  phases  of  the  reflected 
light  can  interfere  constructively  leading  to  greater  reflection,  or  destructivity  leading 
to  reduced  reflection.  This  principle,  known  as  distributed  Bragg  reflection  (DBR), 
provides  a  straightforward  method  for  creating  very  high  reflectivity  or  low  reflectivity 
mirrors  with  a  geometry  that  is  simple  to  understand.  A  one-dimensional  DBR,  also 
known  as  a  Bragg  stack,  is  shown  in  Fig.  3.1.  In  one  dimension,  at  normal  incidence, 
very  high  reflectivity  can  be  achieved  by  choosing  the  thicknesses  of  the  materials  to 
be  equal  to  a  quarter  of  the  optical  wavelength  in  the  material,  The  reflectivity  of 
the  such  a  structure  will  increase  as  the  number  of  periods  and  can  have  bandwidth  of 
a  few  hundred  nanometers.  The  dielectric  constant  for  the  one  dimensionally  periodic 
structure  can  be  simply  described  if  the  dielectric  constant  for  a  single  period  is  known: 
e{y)  =  e{y  +  a).  Such  a  structure  can  also  be  described  not  only  by  the  periodicity  in 
real  space,  but  also  by  the  periodicity  in  the  Fourier  domain,  or  the  /c-space  vectors 
that  can  be  used  to  create  the  perfectly  periodic  lattice.  The  ID  Bragg  stack  in 
Fig.  3.1  could  be  described  by  the  reciprocal  lattice  vector  ky  =  27i/a  where  a  is  the 
periodicity. 

3.2  Planar  photonic  crystal  cavities 

This  periodic  tiling  of  refractive  index  and  principle  of  distributed  Bragg  reflection 
can  also  be  implemented  in  multiple  dimensions.  To  facilitate  fabrication,  it  is  helpul 
to  use  planar  structures,  resulting  in  distributed  Bragg  reflection  in  at  most  two 
dimensions.  Many  2D  periodic  tilings  in  space  exist;  for  photonic  crystals,  the  most 
commonly  used  have  been  the  square  lattice  and  triangular  lattice.  A  triangular 
lattice  2D  photonic  crystal  is  shown  in  Fig.  3.2.  In  this  structure,  for  proper  choice  of 
hole  size  and  periodicity,  the  structure  can  realize  distributed  Bragg  reflection  for  all 
in-plane  wavevectors  over  a  broad  range  of  frequencies;  total  internal  reflection  can 
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Figure  3.1:  One- dimensional  photonic  crystal  with  periodicity  a  and  two  planar  mate¬ 
rials  per  period  (pink  and  purple)  with  different  refractive  indices.  For  correct  choice 
of  thickness  and  refractive  index  in  each  period,  reflections  summed  over  all  interfaces 
add  coherently  to  produce  very  high  reflection. 

confine  light  to  the  slab  in  the  out  of  plane  direction. 

The  exact  dispersion  relation  between  wavevector  and  frequency  can  be  found  by 
solving  Maxwell’s  equations  for  the  structure.  In  practice,  this  is  frequently  done 
by  discretizing  Maxwell’s  equations  and  solving  within  a  single  unit  cell,  while  ap¬ 
plying  periodic  boundary  conditions  to  represent  the  inhnitely  tiled  structure.  One 
popular  approach  is  a  hnite  difference  frequency  domain  simulation,  where  Maxwell’s 
equations  are  cast  into  the  Hermitian  eigenvalue  form [32] 

1 

Vx-VxH=— H  (3.1) 

e  C 

A  second  approach  is  known  as  finite  difference  time  domain  simulation,  where 
an  initial  current  source,  dipole  source,  or  field  pattern  is  allowed  to  evolve  in  time 
according  to  Maxwell’s  equations;  resonant  modes  are  found  by  taking  the  Fourier 
transform  of  the  time-evolved  electromagnetic  helds.  The  result  of  such  a  simulation 
is  shown  in  Fig.  3.3.  The  white  region  indicates  the  area  confined  to  the  slab  by 
total  internal  reflection;  the  gray  region  indicates  modes  in  this  area  can  couple  to  free 
space.  The  lowest  frequency  band  of  solutions  is  known  as  the  “dielectric  band”,  as 
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Figure  3.2:  (a)  Two  dimensional  distributed  Bragg  reflection  in  a  triangular  photonic 
crystal  lattice  of  air  holes  in  a  high  refractive  index  slab.  Light  is  strongly  reflected 
by  the  photonic  crystal  for  all  directions  in  the  plane  of  the  slab,  (b)  Light  is  confined 
inside  the  slab  by  total  internal  reflection  between  the  high  refractive  index  slab  and 
surrounding  low-index  air. 


electric  held  is  concentrated  in  the  higher  dielectric  constant  material  for  these  modes. 
At  higher  frequencies,  there  are  additional  bands  of  solutions,  the  first  of  which  has 
electric  held  primarily  concentrated  in  air,  and  is  known  as  “air  band” ;  higher  bands 
have  solutions  orthogonal  to  the  dielectric  and  air  bands.  In  the  frequencies  between 
the  dielectric  and  air  band,  there  are  no  solutions  to  Maxwell’s  equations;  this  region 
serves  as  a  ’’photonic  band  gap”  preventing  light  from  propagating  inside  the  crystal. 

3.3  Photonic  crystal  cavities 

By  making  a  small  perturbation  to  the  perfectly  periodic  photonic  crystal,  it  is  pos¬ 
sible  to  create  defects  that  can  locally  change  the  frequency  of  a  mode  and  pull  it 
into  the  band  gap  (similar  to  doping  impurities  into  the  electronic  band  structure  of 
semiconductors) .  For  a  small  increase  in  the  size  of  a  single  hole  ( “acceptor”  cavity) , 


CHAPTER  3.  PHOTONIC  CRYSTAL  CAVITIES 


19 


the  frequency  of  the  dielectric  band  is  locally  increased,  forming  a  defect  mode  inside 
the  band  gap  that  cannot  propagate  through  the  crystal.  On  the  other  hand,  de¬ 
creasing  the  size  of  or  eliminating  a  single  hole  locally  decreases  the  frequency  of  the 
electromagnetic  mode,  pulling  a  mode  down  from  the  air  band  into  the  gap  (“donor” 
cavity).  In  this  way,  it  is  possible  to  form  localized  cavities  that  can  have  high  quality 
factor  while  approaching  the  small-volume  limit  of  a  dielectric  cavity  ((A/(2n))^). 

In  the  hrst  chapters  of  this  thesis,  we  will  use  a  donor-type  cavity  design,  where 
three  consecutive  holes  in  a  single  row  are  removed  (also  known  as  the  linear  three 
hole  defect  or  L3  cavity,  developed  by  Susumu  Noda’s  group[33]),  and  the  adjacent 
holes  are  slightly  shifted  to  increase  the  Q  (which  can  be  thought  of  as  adjusting 
the  boundaries  of  the  cavity  to  match  the  phase  [34]).  The  dielectric  constant  for  the 
cavity  and  local  electromagnetic  held  intensity  localized  by  the  mode  are  shown  in 
Fig.  3.4. 

It  is  also  possible  to  form  planar  photonic  crystals  and  cavities  using  distributed 
Bragg  rehection  in  one  dimensional  and  total  internal  rehection  in  two  dimensions, 
for  example,  by  forming  a  narrow  beam  with  holes[35];  the  central  region  of  the  beam 
can  then  be  locally  perturbed  to  form  a  cavity,  as  shown  in  Fig.  3.5. 

3.4  Figures  of  merit  for  photonic  crystal  cavities 

3.4.1  Quality  factor 

The  measure  of  the  quality  of  the  optical  cavity  resonance,  Q  =  is  proportional 
to  the  ratio  of  energy  stored  in  the  cavity  W  to  the  power  radiated  by  the  cavity 
P.  The  Q  can  also  be  written  in  terms  of  the  photon  lifetime  r^/i,  which  can  be 
thought  of  as  the  amount  of  time  the  photon  spends  in  the  cavity  before  it  leaks  out, 
where  Q  =  UTph-  Experimentally,  it  is  often  more  practical  to  measure  the  spectrum 
to  hnd  the  Q,  rather  than  the  lifetime;  this  can  neatly  be  expressed  by  rewriting 
Q  =  ^  =  -k-  Therefore, 


CHAPTER  3.  PHOTONIC  CRYSTAL  CAVITIES 


20 


Band  Diagram  -  Triangular  Lattice 
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Figure  3.3:  Solutions  to  Maxwell’s  Equations  (black  circles)  for  light  polarized  in  the 
plane  of  the  slab  in  the  inhnite  triangular  photonic  crystal  lattice  (scanning  electron 
microscope  (SEM)  of  fabricated  structure  shown  in  top  left).  Lowest  frequency  set 
of  solutions  forms  the  dielectric  band;  next  higher  frequency  set  of  solutions  forms 
the  air  band  (electric  helds  are  plotted  at  the  X  point).  Between  the  frequencies  of 
the  dielectric  band  and  air  band,  there  are  no  guided  modes  in  the  slab.  Boundary 
between  gray  and  white  regions  indicates  the  light  line.  Image  courtesy  of  Jelena 
Vuckovic. 
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Figure  3.4:  (a)  SEM  image  of  an  L3  cavity,  (b)  Local  electric  field  intensity  inside 
the  cavity  at  the  resonant  frequency. 
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Figure  3.5:  (a)  SEM  image  of  an  nanobeam  cavity,  (b)  Ey  inside  the  beam  cavity. 
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The  quality  factor  of  the  L3  cavity  for  index  contrast  of  3.5:1  (e.g.  semiconduc- 
tor:air)  is  approximately  50,000.  Higher  quality  factor  designs  using  smaller  critical 
dimensions  can  also  be  engineered [36,  37].  The  nanobeam  cavity  can  also  be  designed 
to  yield  Q  of  more  than  1  million  [38,  39]. 


3.4.2  Mode  volume 

The  mode  volume  can  be  dehned  as  the  physical  space  over  which  the  mode  would 
extend  if  the  energy  was  uniformly  distributed,  or  alternatively: 


hmode 


/deiElr)!' 


(3,3) 


max  [e(r)|E(r)|^] 

This  L3  cavity  and  nanobeam  cavity  both  have  mode  volume  on  the  order  of  (A/n)^. 


3.4.3  Nonlinear  frequency  conversion 

In  general,  higher  quality  factor  cavities  are  desirable  for  maximizing  the  photon 
storage  time  and  enhancing  the  interaction  between  light  and  the  material  forming 
the  cavity.  Higher  Q  however,  does  ultimately  limit  the  speed  of  any  cavity-based 
device,  as  well  as  the  bandwidth  of  pulses  that  can  be  coupled  into  the  cavity,  so  an 
ultrahigh-Q  is  not  always  desirable.  The  strongest  interaction  will  also  be  obtained 
for  smallest  mode  volumes,  as  the  concentration  of  light  is  enhanced. 

For  each  held  in  nonlinear  three-wave  mixing,  there  will  be  an  enhancement  of 
the  nonlinear  optical  interaction  proportional  to  the  quality  factor  of  the  mode  at 
that  frequency.  That  is,  for  second  harmonic  generation  at  U2  =  2cni,  there  would  be 
a  Q1Q2  enhancement;  for  sum-frequency  generation  the  enhancement  would  be  pro¬ 
portional  to  QiQ2Q3-  a  smaller  mode  volume  will  also  enhance  the  nonlinear  optical 
interaction.  As  the  electromagnetic  held  strength  is  proportional  to  l/i/I4iode,  the 
electromagnetic  energy  concentration  will  be  enhanced  by  a  factor  related  to  l/hAode- 
However,  this  enhancement  will  only  be  realized  if  the  spatial  mode  prohles  of  the 
electromagnetic  held  are  identical  (as  described  by  7  and  (3  in  Chap.  2);  alternatively, 
an  ehective  mode  volume  can  be  dehned[40]. 


Chapter  4 

GaP  photonic  crystal  cavities 

4.1  Motivation 

As  described  in  Chap.  3,  photonic  crystal  nanocavities  conhne  light  in  ultrasmall  vol¬ 
umes,  making  them  ideal  for  low-power,  on-chip  optoelectronic  devices[41,  42,  43]  as 
well  as  for  exploring  fundamental  light-matter  interactions[44,  45].  Most  experimen¬ 
tal  work  in  photonic  crystal  cavities  has  been  done  in  telecom  bands  [45].  This  has 
been  facilitated  by  well  established  fabrication  procedures  for  semiconductors  such  as 
silicon  and  gallium  arsenide,  which  have  band  gaps  in  the  infrared.  Photonic  crystal 
devices  operating  (via  linear  or  nonlinear  processes)  in  the  visible  part  of  the  spec¬ 
trum,  however,  could  serve  as  light  sources  and  spectroscopic  devices  operating  below 
750  nm.  Additionally,  for  devices  operating  at  higher  optical  powers,  such  as  nonlin¬ 
ear  optical  processes  where  the  desired  output  is  superlinearly  dependent  on  incident 
power,  a  material  with  an  electronic  band  gap  in  the  visible  minimizes  two  photon 
absorption  of  telecom  wavelength  inputs.  Furthermore,  photonic  crystal  cavities  with 
resonances  in  the  visible  could  also  be  coupled  to  novel  light  emitters  such  as  nitro¬ 
gen  vacancy  centers[46,  47],  fluorescent  molecules [48],  and  visible  colloidal  quantum 
dots  [49]. 

To  date,  most  photonic  crystal  devices  operating  in  the  visible  have  been  fabricated 
in  one  of  three  material  systems:  gallium  nitride,  silicon  nitride,  and  AlGaInP.  While 
devices  with  wavelengths  as  short  as  blue  can  be  fabricated  in  gallium  nitride  [50, 
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51],  fabrication  processes  in  this  system  are  difficult,  and  the  low  refractive  index  of 
gallium  nitride  (n=2.4)  limits  the  device  quality  factors[52].  Low  refractive  index  is 
also  a  drawback  of  the  silicon  nitride  system [53,  54],  as  is  the  lack  of  a  second-order 
optical  nonlinear  coefficient.  Two  groups  have  fabricated  AlGalnP-based  photonic 
crystals  in  the  red [55,  56]  with  InGaP  quantum  wells.  This  system,  however,  requires 
precise  control  over  the  fraction  of  each  material  in  the  quaternary  system  to  obtain 
the  proper  band  gap  and  lattice  constant,  and  is  lattice  matched  to  a  GaAs  substrate, 
which  is  absorbing  in  the  visible. 

In  this  chapter,  we  show  that  photonic  crystal  devices  at  visible  wavelengths,  which 
could  be  used  as  linear  and  nonlinear  optical  light  sources,  can  be  fabricated  in  gallium 
phosphide.  Gallium  phosphide  is  a  high  refractive  index  (n=3.25  at  700  nm,  n=3.44 
at  555  nm  at  room  temperature  [5 7])  111-V  semiconductor  with  an  indirect  band  gap 
at  555  nm  at  room  temperature.  Although  gallium  phosphide  is  not  typically  used  for 
devices  requiring  very  high  brightness  because  of  its  indirect  band  gap,  incorporating 
InP  quantum  dots  or  quantum  wells  [58],  or  colloidal  quantum  dots  should  greatly 
increase  the  quantum  efficiency.  Additionally,  gallium  phosphide  has  a  relatively  high 
thermal  conductivity  (see  Ghap.  2)  and  a  low  surface  recombination  rate,  estimated  to 
be  roughly  an  order  of  magnitude  less  than  that  of  gallium  arsenide [59]  because  of  its 
high  phosphorus  content,  which  has  a  passivating  effect [60].  Surface  recombination  is 
particularly  important  for  photonic  crystal  devices  because  of  the  additional  exposed 
surface  area  from  etched  air  holes. 


4.2  Cavity  design  and  fabrication 

Linear  three  hole  defect  (L3)[33]  cavities  were  hrst  simulated  and  then  fabricated  in 
a  140  nm  thick  GaP  membrane  (Fig.  4.1).  Simulation  results  indicate  that  quality 
factors  above  10,000  can  be  obtained  for  this  type  of  cavity  (three  holes  removed, 
holes  next  to  cavity  displaced  outward  by  0.15a)  using  a  triangular  lattice  of  lattice 
constant  a  with  slab  thickness  d/a=0.55  and  hole  radius  r/a=0.3  with  no  other  design 
modihcations.  The  electric  held  prohle  for  the  fundamental  L3  cavity  mode  is  shown 
in  Fig.  4.1b. 
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Figure  4.1:  (a)  Scanning  electron  microscope  image  of  a  fabricated  GaP  photonic 
crystal  membrane  after  undercut  of  sacrificial  layer  with  lattice  constant  a=246  nm, 
hole  radius,  r/a=0.32,  slab  thickness  d/a=0.57.  (b)  Finite  difference  time  domain 
simulation  of  electric  held  intensity  inside  the  cavity  for  the  high  Q  mode. 


To  fabricate  these  structures,  the  140  nm  thick  GaP  membrane  was  grown  on 
the  top  of  a  1  pm  thick  sacrihcial  AlxGai_xP  layer.  Patterns  were  hrst  dehned  in 
ZEP520A  resist  by  electron-beam  lithography  and  then  transferred  into  the  GaP 
membrane  by  a  chlorine-based  reactive  ion  etch.  The  recipe  used  was  similar  to  the 
one  listed  in  Appendix  A.  To  calibrate  the  dry  etch  rate  of  small  features  in  GaP, 
trenches  of  various  widths  were  etched  and  then  cleaved  and  measured  in  a  scanning 
electron  microscope  (SEM),  as  shown  in  Fig.  4.2.  As  evidenced  in  the  SEM  images, 
smaller  features  have  less  vertical  etching  prohles.  Excess  photoresist  was  removed 
with  oxygen  plasma,  and  the  sacrihcial  layer  was  undercut  with  hydrohuoric  acid  to 
yield  suspended  membrane  structures  with  high  index  contrast.  For  this  value  of 
aluminum  content  (x=0.6)  in  AlGaP,  the  selectivity  of  the  wet  etch  was  insufficient 
to  remove  the  sacrihcial  AlGaP  layer  without  some  ehect  on  the  GaP  membrane  as 
well.  Scanning  electron  microscope  images  of  one  of  the  fabricated  membranes  can 
be  seen  in  Fig.  4.1;  in  this  case,  wet  etch  time  is  minimized  so  the  ehect  on  GaP 
photonic  crystal  is  not  readily  visible.  For  other  samples  fabricated  with  slightly 
longer  wet  etch  times,  the  circular  photonic  crystal  holes  etch  to  squares  along  the 
crystallographic  directions  of  the  GaP,  as  shown  in  Fig.  4.2.  High  aluminum  content 
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Figure  4.2:  (a)  Scanning  electron  microscope  image  of  a  1  /rm-wide  trench  etched  in 
GaP.  (b)  SEM  of  200  nm-wide  trenches,  (c)  SEM  of  100  nm-wide  trenches. 


Figure  4.3:  SEM  images  of  chips  undercut  a  sufficient  amount  of  time  to  also  etch 
GaP  membrane  photonic  crystal  holes,  transforming  circular  photonic  crystal  holes 
into  squares. 


(x>0.8)  was  subsequently  found  to  eliminate  this  problem  by  increasing  the  selectivity 
of  the  undercut. 


4.3  Reflectivity  measurement  of  cavity  resonances 

The  fabricated  photonic  crystal  resonators  were  probed  using  a  confocal  cross-polarized 
reflectivity  measurement  [61,  62]  technique  as  depicted  in  Fig.  4.4a.  This  technique 
uses  polarization  control  to  achieve  a  high  signal-to-noise  ratio  and  allows  probing 
of  cavities  with  no  internal  light  source.  A  vertically  polarized  (V)  probe  is  directed 
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Figure  4.4:  (a)  Experimental  setup.  The  V-polarized  cavity  is  probed  in  a  cross- 
polarized  setup  using  a  PBS  and  a  HWP.  The  cavity  signal  observed  at  the  output 
follows  a  sin(46')  dependence  where  theta  is  the  HWP  setting  (6*=0  corresponds  to  V 
polarization),  (b)  Spectrum  of  a  cavity  measured  with  conhguration  of  Fig.  4.4(a) 
with  HWP  at  0  =  67°.  The  cavity  resonance  is  at  722  nm  at  room  temperature.  A  £t 
to  a  Lorentzian  (solid  line)  gives  Q=1100.  (d)  (Circles)  Normalized  counts  in  cavity 
peak  as  a  function  of  HWP  angle  (background  subtracted).  (Line)  Fit  with  sin(46'). 
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Figure  4.5:  (a)  Shift  of  photonic  crystal  resonance  at  room  temperature  as  a  and 
r/a  are  changed.  Solid  lines  are  hts  to  Lorentzians.  (b)  Cavity  resonance  measured 
at  room  temperature  (blue)  and  low  temperature  (red).  Solid  lines  indicate  hts  to 
Lorentzians. 
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through  a  polarizing  beam  splitter  (PBS)  and  a  half  wave  plate  (HWP)  onto  the 
photonic  crystal  cavity,  which  has  a  mode  that  is  also  vertically  polarized  (V).  The 
reflected  output  is  observed  through  the  PBS,  which  acts  as  a  horizontal  (H)  polarizer. 
The  primary  polarization  of  the  cavity  mode  is  vertical  (V).  When  the  angle  between 
the  HWP  fast  axis  and  the  vertical  direction  (9)  is  set  to  zero,  the  cavity  coupled 
light  is  reflected  with  V-polarization  and  does  not  transmit  through  the  beamsplitter. 
Rotation  of  the  HWP  allows  part  of  the  cavity-coupled  light  to  be  transmitted  at  the 
PBS  into  the  output  port  with  intensity  following  a  sin(46*)  dependence  (Fig.  4.4(b,c)). 
A  tungsten  halogen  lamp  source  was  used  as  a  wide  broadband  signal  at  the  input. 
The  output  held  was  measured  using  a  spectrometer  with  a  liquid  nitrogen  cooled 
CCD. 

To  fabricate  cavities  with  resonant  wavelengths  spanning  a  large  interval,  lat¬ 
tice  constant  a  (210  nm<  a  <246  nm)  and  hole  radius  (0.25<  r/a  <0.35)  were 
modihed.  These  changes  in  the  fabrication  parameters  resulted  in  cavities  with  res¬ 
onance  wavelength  spanning  from  645  nm  to  750  nm  as  shown  in  Fig.  4.5.  (Plots 
in  Fig.  4.5  are  shown  with  background  subtracted.)  The  measured  resonances  have 
shorter  wavelengths  than  predicted  from  hnite  difference  time  domain  simulations  (by 
about  0.05a/A).  We  believe  this  difference  is  a  result  of  imperfect  selectivity  between 
the  sacrihcial  layer  and  membrane  during  the  wet  etch  step,  which  reduces  the  mem¬ 
brane  thickness.  Quality  factors  measured  experimentally  are  between  500  and  1700. 
This  is  about  one  order  of  magnitude  smaller  than  predicted  by  simulations  because 
of  imperfections  introduced  in  the  fabrication  process. 

The  shortest  resonant  cavity  wavelength  observed  at  room  temperature  was  645 
nm  with  quality  factor  610.  The  resonance  wavelength  depends  on  the  temperature- 
dependent  refractive  index  of  the  membrane.  To  shift  the  cavity  resonance  to  lower 
wavelengths,  the  sample  was  cooled  to  IIK  in  a  continuous  flow  liquid  helium  cryostat. 
The  change  in  temperature  caused  the  resonance  to  shift  from  645  nm  to  641  nm 
(in  the  red  wavelenth  range)  and  the  quality  factor  to  decrease  to  490  as  shown  in 
Fig.  4.5.  With  improved  fabrication,  devices  operating  close  the  555  nm  band  gap 
of  gallium  phosphide  could  be  made;  such  devices  could  be  used  as  green  LEDs  and 
lasers  in  combination  with  embedded  light  emitters  at  these  wavelengths.  Quality 
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factors  could  also  be  improved  by  using  modified  designs  with  simulated  Q  values 
more  than  a  order  of  magnitude  higher [63].  However,  these  Qs  can  still  enable  a 
signihcant  Purcell  effect,  up  to  185,  as  mode  volume  is  only  0.7(A/n)^. 

4.4  Conclusions 

In  summary,  in  this  chapter,  we  demonstrate  a  new  GaP-based  materials  system 
for  photonic  crystals  in  the  visible.  We  probe  linear  three-hole  defect  cavities  in 
reflectivity  and  observe  resonances  from  645  nm  to  750  nm  at  room  temperature  with 
quality  factors  as  high  as  1700.  To  our  knowledge,  these  are  the  shortest  wavelength 
photonic  crystal  cavities  fabricated  in  a  high-index  material. 

We  expect  this  materials  system  to  allow  access  to  a  new  range  of  photonic  crystal 
devices  in  the  visible,  including  LEDs  and  lasers.  Additionally,  the  gallium  phosphide 
system  will  provide  opportunities  to  couple  photonic  crystal  cavities  to  emitters  in 
the  visible  such  as  colloidal  quantum  dots,  nitrogen  vacancy  centers,  biomolecules,  or 
organic  molecules. 


Chapter  5 

Three  wave  mixing  in  GaP 
photonic  crystals 


This  chapter  presents  three  experiments  in  three  wave  mixing  using  GaP  photonic 
crystals.  In  the  hrst,  we  demonstrate  record  efficiency  second  harmonic  generation 
using  photonic  crystal  cavities;  we  then  employ  a  similar  approach  to  demonstrate 
sum- frequency  generation.  Finally,  we  demonstrate  a  broadband  enhancement  of 
second  harmonic  generation  by  using  a  photonic  crystal  waveguide. 

5.1  Motivation 

III-V  semiconductors  such  as  GaAs  and  GaP  are  considered  promising  candidates 
for  nonlinear  optical  devices  [22,  21,  64]  because  of  their  large  second  order  non¬ 
linearity  [65],  transparency  over  a  wide  wavelength  range  (870  nm-17  /rm  for  GaAs 
and  550  nm-11  /rm  for  GaP),  and  ease  of  integration  with  semiconductor  processing. 
The  cubic  symmetry  of  the  zincblende  lattice  of  III-V  semiconductors,  however,  does 
not  exhibit  birefringence,  so  achieving  phase  matching  of  the  two  different  frequen¬ 
cies  typically  requires  employing  quasi-phase  matching  techniques  or  an  additional 
birefringent  material  [20,  66,  67,  68].  In  addition  to  the  challenges  posed  by  phase 
matching,  nonlinear  optical  devices  are  also  constrained  by  the  bulky,  macroscopic  res¬ 
onant  cavities  typically  used  to  enhance  conversion  efficiency  and  reduce  the  required 
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input  power  [69].  High  quality  factor  micro  cavities  have  the  potential  to  achieve  sim¬ 
ilar  conversion  efficiencies  with  a  vastly  reduced  size,  and  could  be  integrated  with 
nanophotonic  technology.  In  these  microcavities,  the  phase  matching  condition  is 
satished  by  the  spatial  overlap  between  the  fundamental  and  second  harmonic  held 
patterns  [24,  70,  71].  Experimentally,  this  resonant  enhancement  has  been  demon¬ 
strated  in  silica  microdisks  [72] ,  where  green  third  harmonic  radiation  was  observable 
with  hundreds  of  microwatts  incident  continuous  wave  IR  power.  Enhanced  second 
harmonic  generation  in  photonic  crystal  cavities  [7]  and  third  harmonic  generation 
in  photonic  crystal  waveguides  [73]  have  also  been  demonstrated.  However,  these 
experiments  suffered  from  poor  conversion  efficiency,  in  part  because  the  second  har¬ 
monic  was  above  the  band  gap  of  the  employed  semiconductor  and  therefore  strongly 
absorbed  in  it  (where  efficiency  is  dehned  as  the  ratio  of  the  output  second  harmonic 
to  the  coupled  input  power  at  1550  nm,  third  harmonic  Pout /.Pin, coupled  =  5  x  10“^° 
for  peak  Pin  =  IW  [73];  second  harmonic  Pout/Pin  =  10“^^  for  CW  P^  =  300  fTW 

[7]). 

5.2  Second  harmonic  generation  in  GaP  photonic 
crystal  cavities 

In  this  section,  we  show  that  photonic  crystal  nanocavities  resonant  with  the  pump 
wavelength  can  be  used  to  generate  second  harmonic  radiation  with  input  power 
orders  of  magnitude  smaller  than  previously  demonstrated. 

5.2.1  Experimental  setup  and  fabrication  of  cavities 

The  experiment  is  described  in  Fig.  5.1(a).  Our  resonator  is  a  modihed  linear  three 
hole  defect  photonic  crystal  cavity  [33].  A  scanning  electron  microscope  (SEM)  image 
is  shown  in  Fig.  5.1(b).  A  tunable  infrared  laser  (Agilent  81680A)  with  wavelength 
range  around  1500  nm  at  normal  incidence  is  spatially  aligned  to  the  cavity  location 
and  spectrally  aligned  to  the  cavity  resonance;  the  laser  polarization  is  also  aligned  to 
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Figure  5.1:  (a)  Microscope-based  setup  for  second  harmonic  generation.  HWP:  half 
wave  plate,  NPBS:  nonpolarizing  beamsplitter,  OL:  objective  lens,  PBS:  polarizing 
beamsplitter,  SPF:  short  pass  hlter,  PD:  photodiode.  The  incident  light  traces  the 
red  line  into  the  cavity  sample.  The  second  harmonic  light  follows  the  blue  line 
into  the  spectrometer,  photodiode,  or  camera.  The  polarization  of  the  incident  light 
is  controlled  by  the  polarizer  and  HWP;  the  polarization  of  the  second  harmonic 
radiation  is  measured  using  HWP  and  PBS.  (b)  SEM  image  of  a  fabricated  structure. 
Scale  bar  indicates  1  pm.  (c)  Spectrum  of  generated  second  harmonic  light  with  8 
nW  power  at  1497.4  nm  coupled  to  the  cavity  (160  nW  incident). 


match  that  of  the  cavity  mode.  In  this  conhguration,  enhanced  second  harmonic  ra¬ 
diation  is  generated;  this  radiation  is  either  analyzed  by  a  spectrometer,  measured  by 
a  femto-Watt  photodetector,  or  imaged  onto  a  camera.  A  second  harmonic  spectrum 
with  8  nW  power  coupled  to  the  cavity  (160  nW  incident  on  the  sample)  is  shown 
in  Fig.  5.1(c).  The  electronic  band  gap  of  GaP  minimizes  both  absorption  of  the 
harmonic  radiation  at  750  nm  and  two  photon  absorption  at  the  fundamental  wave¬ 
length  1500  nm.  The  samples  were  grown  by  gas-source  molecular  beam  epitaxy  on  a 
(lOO)-oriented  GaP  wafer.  A  160  nm  thick  GaP  membrane  was  grown  on  the  top  of 
a  1  pm  thick  sacrihcial  AlGaP  layer.  Structures  were  fabricated  with  e-beam  lithog¬ 
raphy  and  etching,  as  described  in  [8].  The  photonic  crystal  cavities  are  three  hole 
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linear  defects  [33]  with  lattice  constant  a=500-560  nm,  hole  radius  r/a^  0.2  —  0.25, 
and  slab  thickness  d/a  ^  0.3.  We  use  a  perturbation  design  for  our  photonic  crystal 
cavities  [74]  to  increase  the  coupling  efficiency  between  the  cavity  and  objective  lens. 


5.2.2  Simulation  of  second  harmonic  inside  photonic  crystal 
cavities 
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Figure  5.2:  (a)  Finite  difference  time  domain  (FDTD)  simulation  of  electric  held 
inside  the  cavity  for  the  fundamental  TE-like  cavity  resonance  in  the  center  of  the 
slab.  Cavity  held  axes  are  E'^  and  Ey  (b)  Illustration  of  orientation  of  cavity  relative 
to  crystal  axes.  Cavities  axes  E'^,  Ey  are  rotated  from  crystal  axes  E^,  Ey  by  an  angle 
6.  Fields  along  crystal  axes  are  determined  by  projection. 


To  efficiently  radiate  power  at  the  second  harmonic,  the  held  patterns  at  the  fun¬ 
damental  and  second  harmonic  frequency  must  have  constructive  overlap.  In  our 
experiment,  the  fundamental  frequency  is  set  to  that  of  the  cavity  mode,  which  dic¬ 
tates  the  held  pattern  at  this  frequency.  The  second  harmonic  frequency  is  matched 
to  a  higher  order  photonic  band  edge  mode,  which  determines  the  second  harmonic 
held  pattern.  We  simulate  the  held  patterns  inside  the  structure  at  the  fundamental 
and  second  harmonic  wavelengths  using  3D  hnite  diherence  time  domain  methods. 
The  in-plane  electric  held  prohle  in  the  center  of  the  slab  for  the  fundamental  trans¬ 
verse  electric-like  (TE-like)  cavity  mode  is  shown  in  Fig.  5.2(a).  In  the  far-held,  this 
cavity  radiates  primarily  with  y'  polarization,  and  therefore  incident  light  with  this 
polarization  can  couple  into  the  cavity.  GaP  has  a  noncentrosymmetric  cubic  crystal 
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Figure  5.3:  (a)  FDTD  simulation  of  TM-like  photonic  bands  for  the  same  triangular 
lattice  photonic  crystal.  Red  indicates  band  positions.  White  solid  lines  indicate  light 
line;  black  solid  lines  indicate  numerical  aperture  of  lens.  White  box  indicates  mode 
at  second  harmonic  frequency,  a:  lattice  constant  of  photonic  crystal,  (b)  E^  held 
patterns  of  degenerate  TM-like  mode  at  second  harmonic  frequency  at  the  F  point. 


/Q'i 

structure;  the  only  non-zero  elements  of  the  bulk  Xijk  have  i  ^  j  ^  k.  Since  the 
input  signal  in  our  experiment  couples  to  the  TE-like  photonic  crystal  cavity  mode 
[with  dominant  E'^  and  Ey  in-plane  held  components  dehned  in  Figs.  5.2(a),  5.2(b)], 
the  output  second  harmonic  signal  will  then  be  primarily  5-polarized,  i.e.,  transverse 
magnetic-like  (TM-like)  mode.  The  bulk  second  harmonic  nonlinear  polarization  gen- 
erated  takes  the  form  Pz  =  2eQdi4ExEy  where  x,y  are  the  crystal  axes,  Ej.,  Ey  are 
the  corresponding  components  of  electric  held  of  the  coupled  input  signal  (pump), 
eo  is  the  permittivity  of  free  space,  and  di4  is  the  second-order  nonlinear  coefficient 
in  contracted  notation.  The  photonic  band  diagram  for  TM-like  modes  is  shown  in 
Fig.  5.3(a).  Modes  near  the  F  point  {k^  =  ky  =  0)  at  the  second  harmonic  wavelength 
can  efficiently  radiate  into  the  numerical  aperture  of  our  objective  lens  (indicated  by 
black  solid  lines).  The  mode  with  frequency  closest  to  the  second  harmonic  is  in¬ 
dicated.  The  held  patterns  for  this  doubly  degenerate  mode  [Fig.  5.3(b)]  indicate 
that  spatial  overlap  between  the  fundamental  and  second  harmonic  helds  is  best  (i.e. 
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with  the  Ez  field  pattern  that  best  matches  Pz  described  by  the  previous  expres¬ 
sion,  maximizing  |  J  dVE^Pzl)  for  cavities  oriented  at  6^  =  0°,  in  agreement  with  our 
experimental  observations  for  cavities  with  different  6.  The  fundamental  and  second 
harmonic  field  patterns  have  imperfect  spatial  overlap;  future  optimization  of  this 
overlap,  or  introduction  of  a  TM  cavity  mode  at  the  second  harmonic  frequency  [27], 
would  improve  conversion  efficiency. 


5.2.3  Experimental  characterization  of  cavity- enhanced  sec¬ 
ond  harmonic  generation 


Figure  5.4:  (a)  Spectrum  of  fundamental  resonance  probed  in  cross-polarized  reflec¬ 
tivity  with  a  broadband  source.  Lorentzian  fit  gives  a  quality  factor  of  5600.  (b) 
Spectrum  at  second  harmonic  as  exciting  laser  frequency  is  tuned  across  the  cavity 
resonance.  Solid  line  shows  fit  to  Lorentzian  squared  with  cavity  quality  factor  of 
6000. 


We  first  characterize  the  fundamental  resonance  of  the  cavities  by  probing  them 
from  vertical  incidence  using  cross-polarized  reflectivity  with  a  broadband  lamp  [8] 
and  a  tunable  laser.  The  cross-polarized  configuration  is  used  to  obtain  sufficient 
signal-to-noise  to  observe  the  resonance  above  the  reflected  background  uncoupled  to 
the  cavity.  We  select  cavities  with  fundamental  TE-like  modes  at  1480  nm-1560  nm, 
within  the  range  of  our  tunable  laser.  We  estimate  the  coupling  efficiency  into  the 
cavity  to  be  approximately  5%  by  measuring  the  reflected  and  incident  light  power. 
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in  agreement  with  simulations  for  this  type  of  cavity  [74],  (For  higher  input  coupling 
efficiency,  the  cavity  could  be  integrated  with  photonic  crystal  waveguide  couplers 
[75].)  A  typical  reflectivity  spectrum  measured  with  a  tungsten  halogen  white  light 
source  is  shown  in  Fig.  5.2.3(a);  the  measured  quality  factor  is  5600.  Finite  difference 
time  domain  (FDTD)  simulations  indicate  that  the  Q-factor  in  this  structure  is  limited 
to  15000  because  of  the  thin  slab.  The  spectral  prohle  we  observe  is  governed  by 
the  Lorentzian  density  of  states  of  the  cavity  Pci(^)  =  - ,  ,  ^  where  Un 

and  Q  are  the  cavity  frequency  and  quality  factor  respectively.  Once  we  locate  the 
fundamental  resonance  of  the  cavity,  we  use  the  continuous  wave  tunable  infrared  laser 
to  generate  second  harmonic  radiation,  as  depicted  in  Fig.  5.1.  A  typical  spectrum 
of  the  second  harmonic  signal  as  we  scan  the  laser  through  the  cavity  resonance 
is  shown  in  Fig.  5.2.3(b).  The  cavity  enhances  the  coupled  input  power  by  a  factor 
proportional  to  on  resonance;  away  from  resonance,  this  enhancement  drops  off  like 
a  Lorentzian  squared.  On  resonance,  the  total  second  harmonic  power  Pout  emitted 
depends  on  the  square  of  the  input  power  coupled  to  the  cavity  Pin^coupied)  square 
of  the  quality  factor  Q  of  the  cavity,  and  the  overlap  between  fundamental  and  second 
harmonic  mode  held  prohles.  Pout  oc  Pj^  ,coupied<5^l  /  From  the 

second  harmonic  spectrum  [Fig.  5.2.3(c)],  we  can  extract  the  quality  factor  of  the 
cavity  mode  equal  to  6000,  in  good  agreement  with  the  Q  observed  in  the  rehectivity 
measurements  at  the  fundamental  wavelength  (1500  nm). 

To  measure  the  absolute  second  harmonic  power  radiated,  we  send  the  second 
harmonic  signal  to  a  femto-Watt  photodetector.  Figure  5.5  shows  the  measured  out¬ 
put  power  as  a  function  of  the  incident  power  coupled  into  the  cavity,  assuming  an 
input  coupling  to  the  cavity  of  5%.  The  data  show  a  good  ht  to  a  quadratic  depen¬ 
dence  on  input  power.  We  measure  a  normalized  conversion  efficiency  T’out/T’i, coupled 
of  430%/W,  or  Pout /-Pin, coupled  =  5  X  10“®  for  11  pW  coupled  input  power  (Pout/-Pin  = 
2  X  10“®  for  220  pW  power  through  the  objective  lens).  The  measured  value  is  an 
underestimate  of  the  second  harmonic  power  generated  inside  the  cavity,  as  we  only 
measure  radiation  that  is  directed  vertically  and  collected  with  our  objective  lens. 

We  also  investigate  the  dependence  of  the  second  harmonic  signal  on  the  polariza¬ 
tion  of  the  incident  light  [Fig.  5.6(a)].  To  couple  power  into  the  cavity  and  generate 
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Figure  5.5:  (a)  Second  harmonic  power  as  a  function  of  fundamental  wavelength 
power  coupled  into  the  cavity.  We  estimate  coupling  efficiency  into  the  cavity  to  be 
5%.  Solid  line  indicates  fit,  with  slope  2.02,  indicating  quadratic  power  dependence. 
Output  power  measurements  are  corrected  for  measured  losses  from  optics,  but  do 
not  include  corrections  for  collection  efficiency  into  the  objective  lens. 


resonantly  enhanced  second  harmonic  radiation,  the  incident  light  must  have  polar¬ 
ization  aligned  to  the  polarization  the  cavity  radiates.  For  a  cavity  oriented  at  45° 
relative  to  the  crystal  axes,  we  see  maximum  second  harmonic  signal  when  the  input 
polarization  is  aligned  to  the  cavity  mode  polarization,  and  a  reduction  from  this 
maximum  by  a  factor  of  cos"^  {6)  as  the  incidence  angle  is  rotated  away  from  the  cav¬ 
ity,  as  expected  for  a  quadratic  process.  We  observe  a  similar  angular  dependence 
for  cavities  oriented  at  other  angles  to  the  crystal  axes.  Figures  5.6(b)  and  5.6(c) 
show  images  of  second  harmonic  radiation  recorded  on  a  camera  (DVC  710M,  10s 
integration  time),  imaged  through  the  objective  lens,  as  indicated  in  Fig.  5.1(a).  The 
radiation  shows  similar  intensities  for  x'  and  y'  polarizations,  as  expected  for  the  TM 
band  edge  mode. 
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Figure  5.6:  (a)  Dependence  of  second  harmonic  power  on  incident  light  polarization. 
The  horizontal  axis  corresponds  to  the  angle  between  the  input  polarization  and  the 
polarization  of  the  cavity  mode.  Solid  line  shows  £t  to  cos^(6').  (b)  and  (c)  Second 
harmonic  radiation  imaged  on  a  camera  with  polarizer  oriented  in  x'  (b)  and  y'  (c) 
direction  [the  orientation  of  cavity  relative  to  the  axes  is  shown  in  Fig.  5.2(b)].  The 
gray  box  indicates  the  approximate  location  of  the  photonic  crystal  structure,  with 
scale  bar  indicating  approximately  5  /rm.  Camera  integration  time  for  the  same  for 
both  images. 
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5.3  Sum- frequency  generation  in  GaP  photonic  crys¬ 
tal  cavities 

In  this  section,  we  demonstrate  continuous  wave  (CW)  sum  frequency  generation 
(SFG)  using  two  modes  of  a  photonic  crystal  cavity. 

5.3.1  Cavity  design  and  characterization 

Our  structures  are  fabricated  in  a  160-nm  thick  GaP  membrane,  which  is  grown 
on  top  of  1  fim  Alo.ssGao.isP  on  a  (lOO)-oriented  GaP  wafer  by  gas-source  molecular 
beam  epitaxy.  Photonic  crystals  are  defined  by  e-beam  lithography  and  Gl2-based  dry 
etching[8],  followed  by  a  hydrofluoric  acid  undercut  of  the  sacrificial  AlGaP  layer.  The 
cavity  is  a  modihed  three- hole  linear  defect  [33]  with  lattice  constant  a  =560  nm,  r/a  ^ 
0.25,  and  slab  thickness  d/a  ^  0.3.  We  use  a  perturbation  design  for  our  photonic 
crystal  cavities  [74]  to  increase  the  coupling  efficiency  between  the  cavity  and  objective 
lens  to  5%  (verihed  experimentally  by  measurements  of  incident  and  reflected  power). 

A  scanning  electron  microscope  (SEM)  image  of  a  fabricated  structure  is  shown  in 
Fig.  5.7a.  The  cavity  axes  are  oriented  parallel  to  the  crystal  axes. 

We  characterize  the  cavity  resonances  by  normal-incidence  cross-polarized  white 
light  reflectivity,  using  a  tungsten  halogen  lamp  as  a  source [8].  The  cross-polarized 
conhguration  is  used  to  obtain  sufficient  signal-to-noise  ratio  to  observe  the  Lorentzian 
resonance  above  the  reflected  background  uncoupled  to  the  cavity.  The  white  light 
reflectivity  spectrum  (Fig.  5.7b)  shows  several  resonant  modes  (all  with  transverse- 
electric  polarization),  as  expected  for  this  cavity  design [76],  as  well  as  the  electric 
held  intensities  for  the  two  modes  used  in  our  experiment.  Fig.  5.7c  shows  a  ht  of 
the  fundamental  cavity  mode  to  a  Lorentzian,  giving  a  quality  factor  of  3800.  A 
Lorentzian  ht  of  the  higher  order  mode  gives  a  Q  of  220.  (From  hnite  diherence  time 
domain  (FDTD)  simulations,  the  expected  quality  factor  of  the  fundamental  mode  is 
10,000,  limited  by  the  thinness  of  the  membrane.  The  expected  quality  factor  of  the 
higher  order  mode  is  400.)  The  fundamental  mode  is  primarily  polarized  along  the 
y-axis  (Fig.  5.7a);  the  second  mode  is  polarized  primarily  along  the  x-axis. 


CHAPTER  5.  THREE  WAVE  MIXING  IN  GAP  PHOTONIC  CRYSTALS 


41 


5.3.2  Cavity-enhanced  sum  frequency  generation 

The  experimental  setup  and  principle  of  the  experiment  are  shown  in  Figs.  5.8a,  5.8b. 
We  use  two  continuous  wave  lasers,  one  matched  to  the  wavelength  and  polarization 
of  the  fundamental  mode  (1565  nm),  and  one  matched  to  the  higher  order  mode 
(1504  nm);  the  nonlinear  polarization  generates  sum  frequency  generation  at  767  nm, 
as  well  as  some  second  harmonic  generation  at  752  nm  and  782.5  nm.  Light  from 
the  two  lasers  is  combined  on  a  beamsplitter  and  passed  through  an  objective  lens 
(NA=0.5)  and  focused  onto  the  sample;  sum  frequency  and  second  harmonic  signals 
are  collected  by  the  same  objective  and  sent  to  a  spectrometer. 

Sum-frequency  generation  requires  held  components  along  x,  y,  and  2;  axes  (since 
only  Xxyz  7^  0  for  (100)  GaP).  For  input  beams,  we  use  the  two  resonant  modes  of 
the  cavity  to  provide  the  x  and  y  oriented  components.  As  in  our  previous  exper¬ 
iments  on  second  harmonic  generation  (SHG)  in  these  cavities[9],  we  use  a  trans¬ 
verse  magnetic-like  (TM-like)  Bloch  mode  of  the  photonic  crystal  (Fig.  5.8c)  with 
electric  held  primarily  along  the  z-axis  to  outcouple  SFG  light.  Phase  matching  in 
the  cavity  geometry  requires  strong  spatial  overlap  of  modes  at  all  three  frequencies 
.Pout  oc  Pin, coupled, iPin, coupled, 2Q1Q2 1  J"  Xa;j/^Py,(.JiPx,aj2pz,tu3C^G|  ,  where  Pout  Is  the  gen¬ 
erated  sum  frequency  radiation,  Pin,coupied  is  the  power  coupled  into  the  structure  at 
oil  or  022,  and  Qi  and  Q2  are  the  quality  factors  of  the  resonances  at  ui  and  022,  as¬ 
suming  no  depletion  of  the  incident  beams.  Poor  overlap  of  these  three  modes  limits 
the  measured  conversion  efficiency  (PsFG/Pcoupied,cavity  to  approximately  10“®,  where 
P coupled, cavity)  IS  the  coupled  laser  power  in  the  fundamental  mode  (coupled  power  in 
the  higher  order  mode  is  an  order  of  magnitude  higher). 

Fig.  5.9a  shows  a  spectrum  when  both  lasers  are  incident  on  the  cavity;  there 
are  three  peaks,  two  due  to  second  harmonic  of  each  laser  at  the  cavity  wavelength; 
the  central  peak  is  sum  frequency  generation.  We  characterize  the  sum  frequency 
radiation  by  studying  its  dependence  on  the  wavelength  and  power  of  a  single  laser, 
in  this  case  the  laser  aligned  to  the  fundamental  mode  of  the  cavity  (black  box  shows 
second  harmonic  from  this  mode).  Fig.  5.9b  shows  the  power  dependence  of  sum 
frequency  generation  counts  (sum  of  counts  in  green  box)  and  counts  from  second 
harmonic  at  the  fundamental  wavelength  (sum  inside  red  box).  A  polynomial  £t 
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of  second  harmonic  counts  shows  a  power  dependence  of  1.8,  close  to  the  expected 
2.0;  a  £t  to  the  sum  frequency  counts  shows  a  power  dependence  of  0.9,  close  to 
the  expected  1.0.  By  adjusting  the  power,  we  can  also  control  whether  we  observe 
more  light  from  sum  frequency  generation  or  second  harmonic  generation.  Fig.  5.9c 
shows  the  dependence  of  sum  frequency  and  second  harmonic  counts  (green  and 
red  boxes  respectively)  on  the  wavelength  of  the  laser  tuned  near  the  fundamental 
cavity  mode.  The  density  of  states  at  the  output  wavelength  is  determined  by  the 
product  of  the  Lorentzian  densities  of  states  for  each  of  the  participating  input  modes 

^  where  Ac  and  Q  are  the  cavity  wavelength  and  quality  factor 


_  Q\ 


respectively  for  each  mode).  We  observe  this  spectral  profile  by  detuning  the  laser 
from  the  fundamental  wavelength  of  the  cavity,  while  keeping  the  input  beam  at 
the  higher  order  mode  fixed.  For  second  harmonic  generation,  we  fit  the  data  to  a 
Lorentzian  squared,  which  gives  a  Q  of  4100.  For  the  sum  frequency  radiation,  we  fit 
the  data  to  a  Lorenztian,  giving  a  Q  of  4300.  These  Q’s  are  in  good  agreement  with 
the  value  measured  from  cross-polarized  reflectivity.  To  increase  the  tunable  range  of 
the  output  light,  we  can  alternatively  scan  the  laser  at  the  wavelength  of  the  higher 
order  cavity  mode,  which  has  a  lower  Q  (Fig.  5.9d).  Finally,  we  believe  the  sum 
frequency  and  second  harmonic  of  both  cavity  modes  outcouple  through  the  same 
TM  Bloch  mode  of  the  crystal  (Fig.  5.8c),  as  the  radiation  pattern  collected  through 
the  objective  and  imaged  onto  a  camera  looks  similar  at  all  three  wavelengths. 


5.4  Second  harmonic  generation  in  GaP  photonic 
crystal  waveguides 

Previous  work  in  nanophotonic  nonlinear  frequency  conversion)?,  9,  10]  has  fo¬ 
cused  on  using  nanophotonic  cavities  with  high  quality  factor  (Q)  to  generate  a  res¬ 
onant  enhancement  in  conversion  efficiency  at  the  wavelength  of  the  cavity  mode; 
however,  the  narrow  bandwidth  of  a  high-Q  cavity  would  be  unsuitable  for  frequency 
conversion  of  quantum  emitters  with  a  broadband  spectrum  (such  as  the  nitrogen 
vacancy  (NV)  center)  or  frequency  conversion  of  ultrashort  pulses.  Photonic  crystal 
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Figure  5.7:  (a)  Scanning  electron  microscope  image  of  a  fabricated  photonic  crystal 
membrane  after  undercut  of  sacrificial  layer.  Scale  bar  indicates  1  fim.  (b)  White 
light  reflectivity  from  cavity,  with  FDTD  simulations  showing  electric  held  of  modes  of 
interest.  The  two  modes,  both  transverse  electric-like  (TE-like)  are  polarized  orthogo¬ 
nally:  the  fundamental  mode  is  y-polarized;  the  higher  order  mode  is  x-polarized.  (c) 
Fit  to  Lorentzian  of  fundamental  mode  (peak  indicated  by  black  box)  with  Q=3800. 


CHAPTER  5.  THREE  WAVE  MIXING  IN  GAP  PHOTONIC  CRYSTALS 


44 


(a) 


(c) 


Ez 


032(1565  nm) 
03,  (1504  nm)] 


033(767  nm) 


(b) 


Laser  2 


CCD 


olarizer  P^rize^ 


HWP? 


Sample  OL  NPBS 


I 


Of 


VIonochromator  / 


Figure  5.8:  (a)  Schematic  of  sum  frequency  generation  process.  Incident  photons 
with  frequencies  ui  and  U2  are  converted  to  frequency  (b)  Experimental  setup  for 
SFG.  HWP:  half  wave  plate,  NPBS:  nonpolarizing  beamsplitter,  OL:  objective  lens. 
The  incident  light  traces  the  red  line  into  the  cavity  sample;  the  sum  frequency  light 
follows  the  blue  line  into  the  spectrometer.  The  polarization  of  the  incident  light  is 
controlled  by  a  polarizer  and  HWP.  (c)  Simulated  electric  field  pattern  of  degenerate 
TM-like  Bloch  modes  of  the  crystal  at  wavelength  of  sum  frequency. 
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Figure  5.9:  (a)  Sum-frequency  generation  and  second  harmonic  generation  when  two 
lasers  are  used  to  couple  into  modes  shown  in  Fig.  5.7.  Middle  peak  (green  box) 
is  sum  frequency  generation;  right  (black  dotted  and  red  boxes)  and  left  peaks  are 
second  harmonic  generation  from  each  mode.  Black  box  indicates  second  harmonic 
from  fundamental  mode  over  which  laser  is  scanned  to  measure  power  and  wavelength 
dependence  in  Fig.  5.9b  and  5.9c.  (b)  Power  dependence  of  second  harmonic  and  sum 
frequency  generation  as  a  function  of  incident  laser  power  at  wavelength  of  funda¬ 
mental  cavity  mode.  A  polynomial  £t  (plotted  on  log-log  scale)  gives  a  coefficient  of 
0.9  for  sum- frequency  generation  (green  squares)  and  1.8  for  second  harmonic  (red 
circles).  The  incident  power  of  the  second  laser  at  higher  order  cavity  mode  wave¬ 
length  is  6  mW.  (c)  Wavelength  dependence  of  sum  frequency  (green  squares)  and 
second  harmonic  counts  (red  circles)  as  a  function  of  scanning  laser  wavelength  at 
fundamental  cavity  mode.  Fit  to  Lorentzian  of  SFG  counts  gives  Q  of  4300;  £t  to 
Lorentzian  squared  of  SHG  counts  gives  Q  of  4100,  both  matching  Q  of  the  funda¬ 
mental  mode  in  Fig.  5.7,  as  expected.  The  other  laser  is  at  1504  nm.  (d)  Wavelength 
dependence  of  sum  frequency  counts  as  a  function  of  scanning  laser  wavelength  at 
higher  order  cavity  mode.  Fit  to  Lorentzian  of  SFG  counts  gives  Q  of  272,  matching 
Q  of  the  higher  order  mode.  Second  laser  is  at  1567.4  nm. 
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waveguides  offer  an  alternative  to  cavities;  the  waveguide  photonic  dispersion  can 
produce  large  slow-down  factors  greatly  increasing  the  effective  nonlinear  interaction 
length  over  a  broad  wavelength  range [73,  77].  Third  harmonic  frequency  conversion 
has  been  previously  demonstrated  in  a  silicon  photonic  crystal  waveguide  [73];  how¬ 
ever,  the  efficiency  is  limited  by  the  weak  nonlinearity  of  the  material,  linear 
absorption  of  the  third  harmonic,  and  two-photon  absorption  of  the  fundamental  fre¬ 
quency.  Here,  we  demonstrate  second  harmonic  generation  using  the  nonlinearity 
of  gallium  phosphide,  which  has  a  large  electronic  band  gap  (~550  nm)  minimizing 
absorption.  The  circulating  intensity  in  the  waveguide  is  enhanced  by  the  group  in¬ 
dex  Ug,  providing  a  factor  of  enhancement  in  conversion  efficiency  compared  to 
a  conventional  waveguide  with  similar  size.  This  material  is  compatible  with  quan¬ 
tum  emitters  including  InP[78]  and  InGaAs  quantum  dots [79],  NV  centers [80],  and 
fluorescent  molecules [15]. 

5.4.1  Design  and  fabrication  of  photonic  crystal  waveguide 

Fig.  5.10  shows  the  dispersion  diagram  of  transverse  electric  (TE)-hke  modes  (in¬ 
plane  E  held  component  only  at  the  center  of  the  slab)  of  the  W1  waveguide.  The 
W1  waveguide  is  formed  by  removing  one  row  of  holes  from  the  triangular  lattice[81]. 
The  photonic  band  gap  is  indicated  by  the  white  region  between  the  shaded  gray 
areas.  The  structure  supports  two  guided  modes  inside  the  band  gap  with  even 
(mode  of  interest)  and  odd  symmetry  of  the  held  component  relative  to  the  xz- 
plane  containing  the  waveguide  axis.  The  inset  shows  the  held  component  Bz  of  the 
photonic  crystal  waveguide  mode  calculated  from  3D  hnite  diherence  time  domain 
(FDTD)  simulation  at  the  kx='x/a  point  for  the  green  band  (even  mirror  symmetry). 
Here,  no  other  changes  to  the  photonic  crystal  have  been  made,  although  the  group 
velocity  dispersion  could  be  decreased  by  using  modihed  designs  shifting  the  hrst  two 
rows  of  holes  relative  to  the  defect  [82]. 

Fig.  5.10b  shows  a  scanning  electron  microscope  (SEM)  image  of  a  photonic 
crystal  waveguide  fabricated  in  a  164  nm  thick  gallium  phosphide  membrane.  The 
membrane  is  grown  by  gas-source  molecular  beam  epitaxy  on  top  of  a  1  pm  thick  layer 
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of  Alo.85Gao.15P  on  a  (lOO)-oriented  gallium  phosphide  wafer.  The  structures  are  fab¬ 
ricated  by  e-beam  lithography  and  dry  etching  through  the  GaP  membrane  followed 
by  a  wet  etch  to  remove  the  sacrihcial  layer  yielding  a  suspended  membrane.  The 
W1  waveguide  is  fabricated  with  lattice  constant  a=560  nm,  hole  radius  r/a  ~0.25, 
and  is  30  periods  (17  /rm)  long.  To  access  wavevectors  of  the  dispersion  lying  below 
the  light  line  (which  have  low  loss  and  large  group  index)  from  normal  incidence, 
the  structure  includes  two  grating  couplers:  a  circular  \/2n  grating  on  the  left  [83] 
and  a  coupler  on  the  right  formed  by  perturbing  the  lattice  with  periodicity  2a  using 
enlarged  holes [84]  (used  only  for  transmission  measurements).  The  structures  are 
oriented  with  the  waveguide  along  a  [Oil]  crystal  direction.  For  transmission  mea¬ 
surements,  an  objective  lens  with  numerical  aperture  of  0.5  is  positioned  above  right 
grating;  the  light  source  is  focused  slightly  off-center  onto  the  left  grating.  Fig.  5.11a 
shows  transmission  measurements  using  a  tungsten  halogen  white  light  source.  Peri¬ 
odic  peaks  in  transmission  result  from  the  hnite  length  of  the  waveguide  (30  periods) 
which  supports  Fabry  Perot  (FP)-like  reflections  with  peaks  separated  by  X‘^/{2Lng) 
where  L  is  the  length  of  the  waveguide  and  Ug  is  the  group  index[85].  The  linewidths 
of  the  peaks  become  narrower  as  the  wavelength  increases  into  the  part  of  the  dis¬ 
persion  below  the  light  line.  The  derived  group  index  is  shown  in  Fig.  5.11b.  The 
maximum  group  index  measured  is  25,  lower  than  expected  from  simulation  (Fig. 
5.11),  but  still  much  larger  than  in  a  conventional  waveguide. 

5.4.2  Experimental  characterization  of  waveguide  transmis¬ 
sion  and  enhanced  second  harmonic  generation 

To  measure  second  harmonic  generation,  we  couple  light  from  a  tunable  continuous 
wave  telecommunications- wavelength  laser  (Agilent  8 1989 A)  into  the  left  grating. 
Second  harmonic  generated  inside  the  waveguide  scatters  out  of  plane  and  is  collected 
by  the  same  objective  placed  above  the  center  of  the  structure  from  the  entire  field  of 
view  onto  a  camera  for  imaging  or  a  monochromator  and  GGD  for  spectral  analysis. 
The  collected  second  harmonic  radiation  as  a  function  of  incident  laser  wavelength 
is  shown  in  Fig.  5.11c  for  a  different  structure  with  slightly  larger  r/a  that  better 
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matched  the  range  of  our  tunable  laser.  Conversion  is  observed  over  a  range  of  incident 
wavelengths,  with  periodic  peaks  again  observed  from  Fabry  Perot  (FP)  reflections 
at  the  fundamental  wavelength.  The  inset  shows  the  close  match  between  the  second 
harmonic  counts  at  the  FP  peaks  and  the  expected  enhancement,  where  Ug  is 
calculated  from  the  experimentally  measured  spacing  between  FP  peaks  in  the  second 
harmonic  data.  Fig.  5. lid  shows  the  second  harmonic  counts  measured  as  a  function 
of  incident  laser  power  transmitted  through  the  objective.  The  laser  wavelength  is 
1561.5  nm;  the  calculated  group  index  for  this  structure  at  this  wavelength  is  30.  A 
log-log  £t  of  the  data  yields  a  slope  of  1.9,  close  to  what  is  expected  for  a  second-order 
process.  The  estimated  external  conversion  efficiency,  calculated  from  the  maximum 
second  harmonic  intensity  is  10“®  for  2  mW  transmitted  through  the  objective,  or 
5xlO“’^/W.  The  power  at  second  harmonic  wavelength  was  measured  by  separately 
calibrating  the  spectrometer  counts  with  a  laser  at  the  second  harmonic  wavelength 
using  a  power  meter.  This  conversion  efficiency  is  four  orders  of  magnitude  less  than 
that  measured  for  photonic  crystal  cavities  in  GaP[9].  As  a  comparison  to  previous 
work  on  frequency  conversion  in  III-V  semiconductors,  we  note  that  our  device,  which 
has  not  been  optimized  for  incoupling  or  outcoupling  efficiency,  has  only  70,  000 
times  lower  second  harmonic  external  conversion  efficiency  for  1.5  /im  light  than 
past  experiments  (using  GaAs/AlGaAs  waveguides  phase  matched  by  orientation 
patterning [86]  that  require  significantly  more  complicated  MBE  growth)  that  used  a 
ridge  waveguide  300  times  longer. 

5.4.3  Measurement  and  simulation  of  second  harmonic  in  the 
far  field 

Because  the  only  nonzero  elements  of  the  x^xyz  tensor  for  (lOO)-grown  GaP  have 
X  y  z  and  the  guided  mode  of  the  photonic  crystal  waveguide  at  the  fundamental 
wavelength  has  electric  held  in-plane  only,  the  second  harmonic  must  have  electric 
held  out-of-plane,  e.g.  couple  to  transverse  magnetic-like  (TM-like)  mode.  A  3D 
FDTD  simulation  of  electric  held  at  the  second  harmonic  frequency  in  a  periodic 
triangular  photonic  crystal  lattice  is  shown  in  Fig.  5.12a;  the  symmetry  of  the  electric 
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fields  indicates  a  monopole  mode  with  far-held  radiation  pattern  (calculated  from  the 
Fourier  transform  of  the  held  immediately  above  the  slab [87])  shown  in  Fig.  5.12b. 
The  removal  of  a  row  of  holes  to  form  the  W1  photonic  crystal  waveguide  perturbs 
this  mode;  the  resulting  electric  held  is  shown  in  Fig.  3c,  with  far  held  patterns 
in  Figs.  5.12d-f.  The  simulated  quality  factor  of  the  mode  is  200,  more  than  an 
order  of  magnitude  lower  than  for  the  same  mode  unperturbed  by  the  waveguide. 
To  verify  experimentally  that  we  collect  second  harmonic  through  this  mode,  we 
image  the  radiation  pattern  from  the  waveguide  onto  a  camera  (Figs.  5.12g-i)  using  a 
Gian  Thompson  polarizer  to  resolve  diherent  polarizations;  the  resulting  images  show 
very  good  agreement  with  the  simulations.  The  wavelength  of  the  incident  laser  for 
these  measurements  is  1561.5  nm  (a/A  =  0.359,  in  good  agreement  with  simulations 
(Fig.  5.10)). 

5.4.4  Conclusions 

In  summary,  we  demonstrate  efficient  second  harmonic  generation  in  the  visible  from 
photonic  crystal  cavities  with  incident  powers  well  below  1  fiW  as  a  result  of  resonant 
recirculation  of  the  pump  (1500  nm)  light  in  the  GaP  nanocavity.  By  using  a  GaP 
membrane  for  our  structures,  we  minimize  absorption  losses  at  the  output  wavelength 
of  750  nm.  Additionally,  our  device  can  be  integrated  with  a  tapered  optical  hber 
(which  creates  near-held  as  well  as  mechanical  perturbatons)  to  create  a  frequency- 
tunable  second  harmonic  source  with  bandwidth  of  >10  nm  as  described  in  Ref.  [88]. 
We  have  also  shown  continuous  wave  sum  frequency  generation  using  two  modes  of  a 
photonic  crystal  cavity,  which  provides  a  nonlinear  optics-based  on-chip  source  that 
can  easily  be  tuned  by  changing  the  frequencies  of  the  cavity  resonances.  By  matching 
a  high  Q  mode  with  a  low  Q  mode,  we  create  a  source  that  is  tunable  over  10  nm 
range.  Finally,  we  demonstrate  enhancement  of  second  harmonic  generation  in  GaP 
in  a  photonic  crystal  waveguide  with  maximum  group  index  of  30.  We  measure  an 
external  conversion  efficiency  of  5xlO“^/W  or  10“®  for  2  mW  incident  power. 
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Our  structures  could  serve  as  on-chip,  low-power  sources  compatible  with  semi¬ 
conductor  fabrication  processing.  These  results  also  indicate  the  potential  of  micro¬ 
cavities  for  signihcantly  reducing  the  required  input  powers  for  other  experiments  in 
nonlinear  optics,  such  as  sum/difference  frequency  generation  and  parametric  down- 
conversion,  as  well  as  for  on-chip  upconversion  of  weak  infrared  signals  to  visible 
wavelengths [89,  90].  Our  results  using  the  waveguide  geometry  are  promising  for 
frequency  conversion  of  on-chip  integrated  emitters  having  broad  spectra  (e.g.,  NV 
centers)  or  large  inhomogeneous  broadening  (e.g.,  InAs/GaAs  quantum  dots),  as  well 
as  for  frequency  conversion  of  ultrashort  pulses. 
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Figure  5.10:  (a)  Dispersion  diagram  for  TE-like  modes  (E-field  at  the  center  of  the 
slab  in  plane)  for  W1  waveguide.  Parameters  are:  hole  radius  r/a=0.25  and  thickness 
d/a=0.25  where  a  is  periodicity  in  x  direction.  The  white  area  between  the  gray 
shaded  regions  indicates  the  photonic  band  gap  of  the  triangular  lattice  photonic 
crystal.  The  green  line  indicates  the  waveguide  mode  with  even  symmetry  of  the 
Bz  held  component  relative  to  the  x-z  plane  (including  waveguide  axis).  The  blue 
line  indicates  waveguide  mode  with  odd  symmetry.  Circles  indicate  FDTD-calculated 
solutions,  and  solid  lines  indicate  interpolated  bands.  Red  line  indicates  light  line. 
The  inset  shows  the  FDTD  simulation  of  the  Bz  held  component  at  the  center  of  the 
slab  for  the  band  plotted  in  green  at  the  =  vr/a  point  (circled)  with  frequency 
a/A=0.36.  (b)  SEM  image  of  30-periodic  PC  photonic  crystal  waveguide  fabricated 
in  160  nm  thick  GaP  membrane  with  a=560  nm.  Scale  bar  indicates  2  pm.  Circular 
grating  at  left  or  modihed  holes  forming  coupler  (right)  can  be  used  to  couple  into 
waveguide  modes  from  free  space. 
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Figure  5.11:  (a)  Transmission  measurement  through  photonic  crystal  waveguide  using 
white  light  source,  (b)  Group  index  calculated  from  FDTD  simulation  (Fig.  5.10a) 
and  experimental  data  (Fig.  5.11a).  Maximum  group  index  measurement  experimen¬ 
tally  is  25.  (c)  Second  harmonic  intensity  measured  from  a  different  structure  as  a 
function  of  incident  laser  wavelength.  Inset:  second  harmonic  counts  at  the  FP  peaks 
and  calculated  as  a  function  of  wavelength,  (d)  Second  harmonic  counts  measured 
as  a  function  of  incident  laser  power.  Red  line  indicates  linear  £t  of  log-log  data 
with  slope  1.9.  The  laser  wavelength  is  1561.5  nm;  the  calculated  group  index  at  this 
wavelength  is  30. 
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Figure  5.12:  (a)  FDTD  simulation  of  E^  field  component  of  the  monopole  TM- 
like  mode  near  frequency  of  second  harmonic  in  triangular  photonic  crystal  lattice. 
(a/A=0.67)  (b)  Calculated  far-field  radiation  pattern  of  mode  in  (a).  White  circle 
indicates  numerical  aperture  of  lens,  (c)  FDTD  simulation  of  E^  in  photonic  crystal 
waveguide  at  second  harmonic  frequency  (a/A=0.71).  (d)  Calculated  far-field  radi¬ 
ation  pattern  of  \Ex\‘^  for  mode  in  (c).  (e)  Calculated  far-field  radiation  pattern  of 
iF'yP  for  mode  in  (c).  (f)  Calculated  far-field  radiation  pattern  of  for  mode  in 
(c).  (g)  Measured  radiation  pattern  for  \Ex\^.  (h)  Measured  radiation  pattern  for 
\Ey\‘^.  (i)  Measured  radiation  pattern  iF'p. 


Chapter  6 

Fast  quantum  dot  single  photon 
source  triggered  at 
telecommunications  wavelength 


This  chapter  presents  the  demonstration  of  the  fastest  optically  pumped  single  photon 
source  by  using  intracavity  second  harmonic  generation  enhanced  by  a  GaAs  photonic 
crystal  cavity.  This  conhguration  is  simpler  and  faster  than  conventional  optically 
pumped  quantum  dot  single  photon  sources,  which  rely  on  Ti;Sapphire  lasers  that 
operate  at  a  hxed  repetition  rate  of  ~80  MHz. 

6.1  Motivation 

Single  quantum  emitters  such  as  nitrogen  vacancy  centers  [91],  molecules  [48],  semi¬ 
conductor  quantum  dots  [92],  and  atoms  [93]  emit  antibunched  light  suitable  for 
communications  requiring  single  photons[94,  95].  Among  these  systems,  semiconduc¬ 
tor  quantum  dots  have  the  highest  emission  rates  and  can  be  most  easily  integrated 
with  semiconductor  technology,  including  microcavities  with  high  quality  factor,  small 
volume,  and  directional  emission  that  increase  the  emission  rate,  efficiency,  and  in- 
distinguishability  of  the  generated  single  photons.  [96,  92] 

However,  the  generation  rate  of  demonstrated  optically  triggered  quantum  dot 
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single  photon  sources  has  been  limited  by  excitation  (Ti:Sapphire)  lasers  to  around 
80  MHz  [92],  Electrical  excitation  [97,  98]  can  circumvent  this;  however,  resonant 
optical  excitation  [96]  improves  the  indistinguishability  of  output  photons,  and  many 
desirable  microcavity  structures  such  as  photonic  crystal  cavities  have  geometries 
that  are  challenging  to  pump  electrically  [99].  Furthermore,  while  telecommunications 
wavelengths  are  desirable  for  transporting  photons  over  long  distances,  excitation  and 
emission  in  many  quantum  dot  materials  systems,  determined  by  material  parameters, 
occurs  at  much  shorter  wavelengths. 

As  described  in  the  previous  chapter [9,  10],  we  demonstrated  that  photonic  crystal 
cavities  fabricated  in  III-V  semiconductors  with  large  nonlinearities  can  greatly 
enhance  nonlinear  frequency  conversion  efficiency,  as  a  result  of  light  recirculation  in¬ 
side  an  ultrasmall  volume.  Here,  we  apply  a  similar  approach  to  excite  a  single  InAs 
quantum  dot  (with  transitions  ~900  nm)  using  a  commercially  available  telecommu¬ 
nications  wavelength  (~1550  nm)  laser  that  can  serve  as  a  trigger  at  GHz  speeds 
when  paired  with  a  lithium  niobate  electro-optic  modulator. 

6.2  Cavity  design  and  fabrication 

A  scanning  electron  microscope  image  of  the  three-hole  linear  defect  photonic  crystal 
cavity[33]  is  shown  in  Fig.  6.1a.  The  structures  are  fabricated  in  a  164  nm  thick  GaAs 
membrane  using  e-beam  lithography,  dry  etching,  and  HF  wet  etching  of  the  sacrihcial 
layer  beneath  the  membrane,  as  described  previously  [62].  Quantum  dots  are  grown 
by  molecular  beam  epitaxy  in  the  center  of  the  membrane.  The  cavity  axis  is  oriented 
along  a  [Oil]  crystal  direction. The  quantum  dot  density  is  <10  quantum  dots/pm^, 
and  the  quantum  dot  inhomogeneous  broadening  is  around  40  nm,  preventing  multiple 
quantum  dots  from  having  resonances  at  the  same  frequency. 

6.3  Continuous  wave  excitation  of  quantum  dot 

Finite  difference  time  domain  simulations  of  the  electric  held  components  (TE-like 
resonance  with  electric  held  primarily  in  the  plane  of  the  photonic  crystal  slab)  are 
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Figure  6.1:  Photonic  crystal  cavity  design  and  fabrication,  (a)  SEM  image  of  fab¬ 
ricated  suspended  membrane  photonic  crystal  cavity.  Scale  bar  indicates  1  /im.  (b) 
Reflectivity  measurement  of  fundamental  cavity  mode.  Lorentzian  £t  gives  Q  of 
7000.  (c)  (i)  and  (ii)  Simulated  electric  held  components  for  fundamental  mode  of 
L3  photonic  crystal  cavity,  used  for  resonantly  enhanced  upconversion.  (d)  (i)  and 
(ii)  Electric  held  patterns  for  TE  mode  closest  to  emission  frequency  of  quantum  dot. 
The  low  Q  (<100)  mode  is  formed  through  weak  confinement  of  an  air  band  mode. 

shown  in  Fig.  6.1c.  A  rehectivity  spectrum  of  the  cavity  showing  a  resonance  at 
~1500  nm,  measured  with  a  broadband  source  in  the  cross-polarized  conhguration 
(as  in  our  previous  work  [8])  to  maximize  signal  to  noise,  is  shown  in  Fig.  6.1b.  A 
Lorentzian  ht  gives  a  cavity  Q  of  7,000.  GaAs  has  a  noncentrosymmetric  cubic  crystal 
lattice  with  43m  symmetry;  the  only  non-zero  elements  of  the  bulk  xljk  tensor  have 
j  Since  normally  incident  light  couples  to  the  TE-like  photonic  crystal  cavity 
mode  (Fig.  6.1c,  with  dominant  E^  and  Ey  in-plane  held  components),  the  generated 
nonlinear  polarization  (and  accordingly  second  harmonic,  as  described  in  [9])  is  i- 
polarized,  i.e.,  transverse  magnetic-like  (TM-like)  mode.  For  a  photonic  crystal,  this 


CHAPTER  6.  FAST  QUANTUM  DOT  SINGLE  PHOTON  SOURCE 


57 


corresponds  to  guided  resonances  in  the  TM-like  air-band  [9].  The  second  harmonic 
at  750  nm  is  used  to  excite  the  quantum  dot  above  the  GaAs  band  gap.  Anisotropy 
of  the  quantum  dot  causes  the  dot  to  emit  with  TE  polarization;  at  the  quantum  dot 
wavelength  (~900  nm),  the  photonic  crystal  supports  an  optical  (air-band)  TE-like 
mode  (Fig.  6. Id),  weakly  localized  by  the  perturbation  from  removing  and  shifting 
holes  near  the  cavity. 

6.4  Pulsed  excitation  triggered  single  photon  source 

A  tunable  telecommunications  wavelength  continuous  wave  laser  (Agilent  8 1989 A, 
1463  nm  -  1577  nm)  is  coupled  into  the  fundamental  high-Q  mode  of  the  cavity. 
Photons  upconverted  in  the  cavity  through  second  harmonic  generation  excite  the 
quantum  dot  above  band;  the  emitted  single  photons  from  the  quantum  dot  are  spec¬ 
trally  hltered  and  sent  to  a  spectrometer  or  Hanbury  Brown-Twiss  (HBT)  setup,  in 
which  photons  are  split  by  a  50/50  beamsplitter  to  two  single  photon  counters [92],  for 
photon  statistics  analysis.  The  telecom  wavelength  laser  can  also  be  modulated  for 
triggered  single  photon  source  operation.  For  direct  recombination  lifetime  measure¬ 
ments  of  the  quantum  dot,  we  alternatively  use  an  ultrashort  pulse  (3  ps)  Ti;Sapphire 
laser  to  excite  the  dot  and  a  streak  camera  for  detection. 

Continuous  wave  (CW)  excitation  is  shown  in  Figure  6.2.  Fig.  6.2a  shows  the 
photoluminescence  spectrum  measured  when  the  continuous  wave  tunable  laser,  on 
resonance  with  the  cavity  mode  (at  ~1500  nm),  is  normally  incident  on  the  structure. 
Spectral  hlters  are  used  to  create  an  effective  bandpass  hlter  (~5  nm  FWHM)  centered 
at  898  nm.  The  measured  counts  from  the  two  strongest  lines  (red  box.  Fig.  6.2a) 
both  show  linear  dependence  on  pump  power  at  low  optical  power;  we  believe  these 
may  be  different  charge  states  of  the  single-exciton[100]. 

To  determine  the  maximum  speed  at  which  the  system  can  be  modulated,  we 
perform  an  independent  experiment  directly  measuring  the  quantum  dot  lifetime 
with  a  streak  camera  (2  ps  timing  resolution),  using  a  Ti; Sapphire  laser  at  750  nm 
with  80  MHz  repetition  rate  and  3  ps  pulses  to  excite  the  dot  (Fig.  6.2a,  inset).  To 
accumulate  sufficient  counts,  data  from  peaks  shown  in  the  red  box  in  Fig.  6.2a  were 
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summed.  The  data  were  fit  to  a  monoexponential  decay  with  time  constant  2.4±0.1 
ns  (where  given  error  is  one  standard  deviation). 

To  verify  that  the  upconverted  light  from  the  CW  tunable  laser  could  be  used  to 
efficiently  excite  a  single  quantum  dot,  we  perform  a  photon  correlation  measurement 
of  the  signal  collected  through  our  spectral  hlter  (Fig.  6.1a)  at  ~898  nm  (QD  emission 
wavelength)  using  the  HBT  setup.  Figure  6.2b  shows  the  histogram  of  coincidence 
counts  as  a  function  of  time  delay  r  between  the  two  detectors.  The  data  were 
£t  to  the  unnormalized  correlation  function  =<  I {t  +  t) I (t)  >=  ^[l  —  (1  — 

[100],  where  A,  tq  and  5'*'^^(0)  are  fitting  parameters,  and  I/tq  =  F  +  rp 
where  F  is  the  spontaneous  emission  lifetime  and  is  the  pnmp  rate.  The  symmetric 
antibnnching  in  coincidence  connts  at  zero  time  delay  g*'^^(0)=0.43±0.04,  normalized 
snch  that  g*'^^(cxD)=l,  indicates  emission  primarily  from  the  same  initial  state  of  a 
single  qnantum  dot  (since  g*'^^(0)  <  0.5)  [100].  Backgronnd  subtraction  for  a  signal  to 
noise  ratio  of  S/N  ^  10  gives  g^^^(0)=0.31±0.05.  An  exponential  fit  to  the  central 
antibnnching  dip  gives  a  decay  rate  of  2.3±0.2  ns;  by  comparison  with  independent 
lifetime  measnrements,  this  indicates  that  the  qnantum  dot  is  pumped  in  the  low 
power  regime,  far  below  satnration. 

To  realize  a  triggered  single  photon  sonrce,  the  telecom  tnnable  laser  is  modulated 
by  a  gigahertz  electro-optic  modulator  (JDSU),  which  is  driven  by  a  pulse  pattern 
generator  (Anritsu  1800A)  that  receives  an  external  clock  signal  from  a  synthesizer 
sweeper  (Agilent)  (Fig.  6.3).  Photon  correlation  measurements  were  performed  with 
the  same  HBT  setnp  as  in  the  CW  experiment  (Fig.  6.2).  The  repetition  rate  of  the 
sonrce  was  varied  between  100  and  300  MHz,  with  dnty  cycle  20-50%  (minimum  duty 
cycle  for  which  the  experiment  conld  be  performed  was  limited  by  backgronnd  from 
the  modulator  during  the  off  portion  of  the  cycle;  the  finite  excitation  dnration  likely 
increases  g*^^^(0)),  as  a  resnlt  of  the  qnantnm  dot  re-excitation.  A  schematic  of  the 
pulsed  excitation  single  photon  sonrce  measnrement  is  shown  in  Fig  6.4. 

Fig.  6.5a  shows  the  measnred  photon  statistics  for  a  repetition  rate  of  100  MHz. 
The  peaks  were  fit  to  exponentials  with  a  single  decay  rate,  and  the  ratio  of  the  areas 
of  the  fitted  peaks  was  used  to  determine  g‘^^^(0)=0.49±0.07  (g‘^^^(0)=0.38±0.08  with 
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background  subtraction,  indicating  multi-photon  probability  suppression  to  38%  rel¬ 
ative  to  an  attenuated  laser  of  the  same  power).  We  observe  a  reduction  in  g^^^  for  the 
immediately  adjacent  peaks  (up  to  two  repetition  periods)  to  the  r  =  0  peak  to  a  value 
intermediate  between  the  value  for  the  central  peak  and  the  saturated  value.  This  long 
timescale  negative  correlation  may  be  caused  by  memory  effects  resulting  from  long- 
lived  charged  states  in  the  quantum  dot  [101].  Fig.  6.5b  shows  photon  correlation 
measurements  for  repetition  rate  of  300  MHz;  g(2)(o)=0.40±0.09  (g(2)(0)=0.3±0.1 
with  background  subtraction).  At  300  MHz,  coincidences  from  consecutive  pulses  be¬ 
gin  to  overlap,  as  a  result  of  the  excitation  repetition  period  approaching  the  quantum 
dot  lifetime. 


6.5  Conclusions 

In  conclusion,  we  have  demonstrated  an  optically  triggered  single  photon  source  op¬ 
erating  at  100  MHz  (with  suppression  of  g*-^^(0)  visible  up  to  300  MHz),  excited 
with  a  telecommunications  wavelength  (~1500  nm)  laser  followed  by  an  electro-optic 
modulator.  In  our  system,  the  rate  of  single  photons  is  limited  primarily  by  the  spon¬ 
taneous  emission  rate  of  the  quantum  dot.  A  doubly  resonant  cavity  with  resonances 
at  both  telecommunications  frequency  and  the  frequency  of  the  dot  would  increase 
the  spontaneous  emission  rate  of  the  dot  via  the  Purcell  effect  [102],  and  enable  the 
realization  of  a  significantly  faster  (e.g.  >1  GHz)  triggered  single  photon  source. 
Photons  emitted  from  the  quantum  dot  could  also  be  frequency  converted  back  to 
telecommunications  wavelengths  via  difference  frequency  generation  using  the  non¬ 
linearity  of  the  surrounding  GaAs  material  [103]  providing  a  full  interface  between 
the  ~900  nm  quantum  node  and  the  optical  fiber  network. 


CHAPTER  6.  FAST  QUANTUM  DOT  SINGLE  PHOTON  SOURCE 


60 


•  M  •  3#  *  *  • 

k' j  ” 

40 


Figure  6.2:  Characterization  of  single  quantum  dot  excited  with  on-chip- up  converted 
1500  nm  laser,  (a)  Spectrum  measured  from  CW  second  harmonic  excitation  of  quan¬ 
tum  dot.  Inset:  Streak  camera  measurement  of  quantum  dot  lifetime.  Counts  were 
summed  over  spectral  window  indicated  by  box.  (b)  Photon  correlation  measurement 
of  quantum  dot  emission  under  frequency  doubled  CW  1550nm  excitation.  Fit  gives 
g(2)(o)=0.43±0.04. 
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Figure  6.3:  Experimental  setup.  A  tunable  C-band  laser  (red  excitation  path)  is 
modulated  by  an  electro-optic  modulator  (EOM)  driven  by  a  pulse  pattern  generator 
(PPG),  which  receives  an  external  clock  signal  from  a  synthesizer.  The  laser  polar¬ 
ization  is  adjusted  to  match  that  of  the  cavity  by  a  half  wave  plate  (HWP).  Pulses 
incident  on  the  cavity  via  a  non-polarizing  beamsplitter  (BS)  excite  quantum  dots  by 
above  band  absorption  of  laser  pulses  upconverted  through  second  harmonic  genera¬ 
tion  in  the  surrounding  GaAs  matrix.  Emission  from  the  quantum  dot  (green  path) 
is  spectrally  hltered,  then  analyzed  either  by  an  HBT  setup  for  photon  correlation 
measurements  or  a  spectrometer.  For  direct  lifetime  measurements  of  the  quantum 
dot,  an  ultrashort  pulse  (3  ps)  Ti:Sapphire  laser  at  80  MHz  repetition  rate  excites 
the  dot  (blue  excitation  path);  a  streak  camera  is  used  for  detection.  Inset  shows 
scanning  electron  microscope  image  of  photonic  crystal  cavity. 
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Figure  6.4:  Schematic  illustrating  principle  of  fast  upconversion-based  telecom 
wavelength-triggered  quantum  dot  single  photon  source.  Pulses  at  1500  nm  are  cou¬ 
pled  into  the  cavity  resonance,  leading  to  strongly  enhanced  second  harmonic  gener¬ 
ation.  The  light  at  the  second  harmonic  frequency  can  create  carriers  in  the  GaAs 
membrane  that  relax  into  the  quantum  dot,  which  are  then  emitted  as  single  photon 
pulses.  The  rate  of  single  photon  pulses  is  determined  by  the  rate  at  which  the  laser 
is  modulated. 
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Figure  6.5:  Photon  correlation  measurement  of  quantum  dot  emission  when  triggered 
with  frequency  doubled  telecom  wavelength  laser  pulses  from  an  externally  modulated 
telecom  wavelength  laser,  (a)  Second  order  autocorrelation  function  measurement 
for  100  MHz  repetition  rate  with  duty  cycle  20%.  g*^^^(0)  =  0.48  <  0.5  indicates 
emission  from  a  single  quantum  dot.  (b)  Second  order  autocorrelation  measurement 
for  repetition  rate  300  MHz  with  duty  cycle  50%  and  g‘^^^(0)=0.40. 


Chapter  7 

Multiply  resonant  high  quality 
photonic  crystal  nanocavities 


This  chapter  presents  the  design,  fabrication,  and  linear  and  nonlinear  characteriza¬ 
tion  of  photonic  crystal  cavities  with  multiple  resonances  separated  by  large  frequency, 
which  are  intended  for  use  in  frequency  conversion  devices. 

7.1  Motivation 

State  of  the  art  photonic  crystal  nanocavity  designs[36,  38,  39,  37]  can  be  optimized 
to  generate  quality  factors  exceeding  one  million  for  a  single  cavity  resonance.  For 
nonlinear  optical  interactions  such  as  frequency  conversion  or  stimulated  Raman  scat¬ 
tering,  however,  it  is  desirable  to  have  multiply  resonant  nanostructures [24]  with  ar¬ 
bitrary  frequency  separation  and  with  good  spatial  held  overlap.  To  enable  nonlinear 
coupling  of  cavity  resonances,  structures  can  be  fabricated  in  a  III-V  semiconductor 
such  as  GaP  or  GaAs  (for  second  or  third  order  optical  nonlinearity)  or  a  group  IV 
semiconductor  such  as  silicon  (for  third  order  nonlinearities,  such  as  Raman  scattering 
or  four  wave  mixing),  as  described  previously  in  this  thesis.  Past  experimental  work  in 
three  and  four-wave  mixing  in  semiconductor  photonic  crystal  nanocavities  has  been 
implemented  using  designs  featuring  a  single  photonic  band  gap  [9,  10,  7,  104,  105]. 
Within  a  single  photonic  band  gap,  however,  it  is  difficult  to  independently  control 
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their  frequencies;  moreover,  field  patterns  of  different  resonances  typically  have  min¬ 
imal  spatial  overlap,  and  the  absolute  frequency  separation  between  resonances  is 
limited  by  the  size  of  the  photonic  bandgap.  Photonic  crystals  and  quasicrystals [26] 
can  have  multiple  photonic  band  gaps,  but  the  size  of  the  higher  order  band  gaps 
greatly  diminishes  for  hnite  thickness  structures;  in  addition,  in  such  planar  struc¬ 
tures,  higher  order  band  gaps  are  located  above  the  light  line,  implying  that  the 
resonances  at  those  frequencies  would  also  have  low  quality  factors  (Q),  resulting 
from  the  lack  of  total-internal  reflection  conhnement.  Band  gaps  for  different  polar¬ 
izations  (e.g.  transverse  magnetic  (TM)  and  transverse  electric  (TE))[23,  27,  103,  106] 
can  generate  additional  resonant  modes;  however,  it  is  difficult  to  independently  tune 
their  frequencies,  and  TM  resonances  require  relatively  thick  membranes  that  are 
more  difficult  to  fabricate. 

Recently [13],  we  proposed  a  crossed  beam  photonic  crystal  cavity  suitable  for  non¬ 
linear  frequency  conversion  that  allows  at  least  two  individually  tunable  resonances 
with  a  frequency  separation  larger  than  the  size  of  the  photonic  bandgap  in  a  single 
nanobeam.  Here, 

To  summarize,  for  optimal  performance  in  frequency  conversion-based  devices,  a 
design  must  allow  for: 

•  Independent  tuning  of  individual  cavity  resonances  with  respect  to  each  other 
through  varying  geometric  parameters  in  the  photonic  crystal 

•  Large  maximum  separation  between  cavity  frequencies 

•  High  spatial  overlap  of  electric  helds  in  the  nonlinear  material 

•  High  quality  and  small  mode  volume  cavities 

•  Separate  channels  for  waveguide  coupling  at  each  frequency 

•  Non- zero  effective  nonlinear  tensor  for  the  polarizations,  material,  and  material 
crystal  orientation  of  choice 

In  this  chapter,  we  propose  and  experimentally  demonstrate  a  crossed  nanobeam 
photonic  crystal  cavity  design  that  can  satisfy  these  requirements  allows  at  least  two 
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individually  tunable  resonances  with  a  frequency  separation  larger  than  the  size  of  the 
photonic  bandgap  in  a  single  nanobeam.  We  also  describe  the  experimental  linear  and 
nonlinear  characterization  of  structures  with  resonances  that  are  nearly  degenerate 
with  orthogonal  polarization,  as  well  as  the  design  and  linear  characterization  of 
structures  with  resonances  separated  by  more  than  500  nm. 


7.2  Design  of  multiply  resonant  photonic  crystals 


(a) 


(c)  (d) 


Figure  7.1:  (a)  Schematic  illustration  of  multiply  resonant  orthogonal  nanobeam 
cavity,  (b)  Illustration  of  photonic  nanobeam.  Red  box  shows  unit  cell,  which  is 
tiled  periodically  in  the  x  direction.  Parameters  are  periodicity  a,  width  w,  hole 
size  hx  and  hy  and  thickness  out-of-plane  t  (not  shown),  (c)  Wavelengths  of  dielec¬ 
tric  and  air  bands  of  GaAs  nanobeam  as  lattice  constant  is  varied.  Parameters  are 
w/a=1.65,  hi/w=0.6,  /i2/a=0.5.  Solid  lines  are  plotted  for  t=160  nm  (fixed  absolute 
slab  thickness)  and  wavelength-dependent  index  of  refraction  n;  dotted  lines  are  plot¬ 
ted  for  t/a  =  0.35  (hxed  relative  slab  thickness)  and  n=3.37.  Field  patterns  of  Ey 
at  kx  =  n/a  for  dielectric  (top)  and  air  band  (bottom)  are  indicated  by  black  circles 
and  arrows,  (d)  Normalized  frequencies  of  dielectric  and  air  bands  (left  axis)  of  GaAs 
nanobeam  as  beam  width  is  varied.  Parameters  are  hi/w=0.6,  h2/a=0.5,  t/a=0.35. 
Right  axis  shows  change  in  size  of  photonic  band  gap  with  beam  width. 
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An  illustration  of  the  design  in  shown  in  Fig.  7.1a.  The  basic  element  of  our  design 
is  the  nanobeam[38,  39,  35,  107,  34,  108],  a  ID  periodic  photonic  crystal  waveguide 
clad  in  the  other  two  directions  by  air,  shown  in  Fig.  7.1b.  A  single  unit  cell  of  a 
single  nanobeam,  comprising  a  rectangular  airhole  with  dimensions  hx  and  hy  in  a 
high  index  dielectric  slab  of  width  w,  is  shown  inside  the  red  box  (thickness  t  of  the 
air-bridged  membrane  is  out  of  the  plane).  The  structure  has  periodicity  a  in  the  x 
direction;  this  periodicity  produces  a  ID  TE  photonic  band  gap  with  size  determined 
by  the  difference  in  frequency  between  the  first  two  bands  (dielectric  band  and  air 
band)  at  the  wavevector  kx  =  n /a.  We  first  investigate  the  periodicities  required  in  a 
nanobeam  structure  to  achieve  different  wavelength  resonances  using  3D  simulations 
performed  using  plane  wave  expansion  with  supercell  approach  (MIT  Photonic  Bands 
(MPB))[109].  Fig.  7.1c  plots  the  wavelengths  of  the  dielectric  and  air  bands  obtained 
from  simulations  as  a  function  of  lattice  constant  a.  Solid  lines  show  wavelengths 
for  hxed  absolute  thickness  of  the  membrane  and  are  calculated  using  the  refractive 
index  n  of  GaAs  at  each  wavelength;  dotted  lines  are  calculated  by  scaling  the  result 
for  wavelength  1550  nm  (i.e.  fixed  t/a,  fixed  n).  Because  a  real  structure  has  a  hxed 
membrane  thickness,  decreasing  the  periodicity  in  one  beam  relative  to  the  other  leads 
to  larger  relative  thickness  t/a,  redshifting  the  wavelengths  of  the  bands;  therefore 
achieving  resonances  with  relative  frequency  /2//1  requires  superlinear  scaling  of  the 
feature  size,  e.g.  02/01  >  /2//1,  as  shown  by  the  divergence  between  solid  and  dotted 
lines.  (There  is  a  second  less  important  contribution  because  refractive  index  increases 
for  higher  frequencies,  from  3.37  to  3.53  for  range  plotted  in  Fig.  7.1c).  Fig.  7. Id  shows 
the  frequencies  of  dielectric  band  and  air  band,  as  well  the  relative  size  of  the  band 
gap,  as  the  width  of  the  beam  is  changed. 

A  cavity  is  formed  in  each  beam  by  introducing  a  central  region  with  no  holes 
(cavity  length  /)  and  tapering  the  lattice  constant  and  hole  size  near  the  cavity 
region[34,  39].  In  each  beam,  confinement  along  the  periodic  direction  is  provided 
by  distributed  Bragg  reflection;  conhnement  out  of  plane  is  provided  by  total  internal 
reflection.  Conhnement  in  the  in-plane  direction  orthogonal  to  the  beam  axis  is  pro¬ 
vided  by  total  internal  rehection,  and  in  the  case  of  beams  with  overlapping  photonic 
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band  gaps,  also  by  distributed  Bragg  reflection  (as  in  structures  designed  for  minimiz¬ 
ing  crosstalk  in  waveguide  intersections  [110]).  This  structure  allows  nearly  indepen¬ 
dent  tuning  of  each  resonant  frequency  by  tuning  the  parameters  (e.g.  width,  lattice 
constant,  cavity  length)  of  each  beam.  Additionally,  the  structure  has  natural  chan¬ 
nels  for  coupling  through  each  beam  to  an  access  waveguide  at  each  wavelength  [111]. 
To  optimize  the  design,  parameters  of  cavity  length,  lattice  constant  (a),  number  of 
taper  periods  N,  distance  between  holes  in  the  taper  region  (oj,  i  =  1,2,  ...N),  hole 
size  in  the  taper  region  (di,  i  =  1,2, ...,  N),  and  beam  width  (w)  were  varied  in  each 
beam.  Fig.  7.2c  shows  the  3D  finite  difference  time  domain  (FDTD)-simulated  Ey 
held  pattern  for  a  resonance  at  1.55  pm,  with  mode  volume  0.35(A/n)^  where  n  is 
the  refractive  index,  and  Q=19000,  limited  by  loss  in  the  vertical  direction.  Fig.  7.2d 
shows  the  held  pattern  of  E^  for  a  resonance  primarily  localized  by  the  vertical  beam 
with  wavelength  1.1  pm  (also  limited  by  vertical  loss)  and  mode  volume  0.47(A/n)^. 
The  shorter  wavelength  mode  likely  has  a  lower  quality  factor  in  part  due  to  its  nar¬ 
rower  width  than  the  horizontal  beam  and  partly  due  to  FDTD  discretization  error 
(12  points  per  period).  Decreasing  the  vertical  lattice  constant  further  (10  points  per 
lattice  period)  leads  to  resonances  with  separation  of  574  nm  (1550  nm  and  976  nm). 

Fig.  7.3a  shows  tuning  of  the  quality  factor  of  an  individual  resonance  calculated 
by  FDTD  by  adjusting  the  cavity  length  (a  minimum  cavity  length  is  required  to  avoid 
overlap  between  the  holes  forming  cavities  in  each  beam).  The  resonant  wavelength 
as  well  as  quality  factor  of  the  beam  can  also  be  tuned  by  varying  other  parameters; 
Fig.  7.3b  shows  the  change  in  resonant  frequency  and  quality  factor  as  width  of  the 
horizontal  beam  is  varied.  To  study  the  overall  effect  f  varying  a  single  parameter,  we 
plot  a  figure  of  merit  (FOM)  for  nonlinear  frequency  conversion  FOM  =  Q1Q2/ \/VA^ 
in  Fig.  7.3c  (right  axis),  which  is  seen  to  depend  primarily  on  Q2-  As  mentioned 
previously,  for  large  horizontal  beamwidths,  the  quality  factor  of  the  higher  frequency 
mode  is  limited  by  diffraction  into  the  orthogonal  beam;  this  principle  is  illustrated 
in  Fig.  7.3(d)  and  Fig.  7.4. 

For  nonlinear  frequency  conversion  applications,  it  is  important  for  cavity  held 
patterns  to  have  large  spatial  overlap.  Defining  the  nonlinear  overlap,  normalized  to 
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1,  as[27] 


Xvl 

^  j  dVe\E^\^^  j  dVe\E2\^ 


(7.1) 


where  NL  indicates  nonlinear  material  only,  we  calculate  7=0.02  for  the  structures 
shown  in  Fig.  7.2.  This  number  could  be  increased  to  0.07  by  decreasing  the  number  of 
taper  periods  from  5  to  3  to  further  localize  the  held  to  the  central  region;  however, 
the  quality  factors  are  reduced  to  1440  and  1077  respectively,  although  this  could 
likely  be  increased  by  reoptimization  of  other  parameters. 


7.3  Degenerate  frequency  orthogonal  polarization 
photonic  crystal  cavities 

A  cavity  with  orthogonally  polarized  resonances  degenerate  in  frequency  can  be 
formed  by  using  the  same  parameters  for  each  beam.  This  could  be  used  for  ap¬ 
plications  such  as  coupling  to  spin  states  of  embedded  quantum  emitters [112]  or  for 
building  polarizaton  entangled  photon  sources  based  on  a  bi-excitonic  cascade  from 
a  single  quantum  dot[113,  114].  Figs.  7.5c,  7.5d  show  the  held  patterns  from  3D 
hnite  diherence  time  domain  (FDTD)  simulations  for  a  doubly  degenerate  structure. 
The  cavities  also  have  several  additional  higher  order  modes  the  number  of  which 
is  determined  primarily  by  cavity  length  (two  per  beam  for  the  structures  shown  in 
Fig.  7.5)  with  additional  held  pattern  nodes.  The  design  of  a  degenerate  frequency 
cavity  is  shown  in  Fig.  7.5. 

The  value  of  7  for  the  degenerate  cavity  is  0.02.  The  value  of  7  for  a  single  reso¬ 
nance  of  the  degenerate  structure  overlapped  with  itself  is  0.89  (i.e.  7  is  determined 
by  the  fraction  of  held  in  the  nonlinear  material),  indicating  the  size  of  7  is  limited  by 
the  amount  of  held  concentrated  in  the  central  region.  This  overlap  is  smaller  than 
that  for  TEoo  and  TMqo  modes  in  a  single  nanobeam[27];  however,  our  design  uses 
a  thinner  membrane  which  is  easier  to  fabricate  and  can  support  larger  frequency 
separations. 


CHAPTER  7.  MULTIPLY  RESONANT  PHOTONIC  CRYSTAL  CAVITIES  70 


4.x 

1 

<o 

CD 

CZ? 

CD 

O 

o 

2 

o 

o 

<z> 

CD 

<o 

CD 

CD 

O 

o 

CD 

o 

0=19000 

=0.35(A/n)^ 

DDCWDto 

•ttDIODD 

CD 

c=. 

CD 

Q=1900 

CD 

O 

'=0.47(A/n)^ 

OOOOOQo® 

f 

oooQOOQQ 

Q 

CD 

C3 

O 

CD 

CD 

CD 

{=5 

Figure  7.2:  (a)  Illustration  of  cavity  design,  showing  intersecting  orthogonal 

nanobeams  with  taper  and  mirror  regions,  as  well  as  central  cavity,  (b)  Detail  of 
white  box  in  (a).  Parameters  used  to  form  resonance,  shown  for  cavity  in  horizontal 
(subscript  h)  beam:  R  indicates  cavity  length;  dhx,N  and  dhy,N  indicate  hole  sizes  in 
mirror  region;  dhx,i  and  dhy,i  indicate  hole  sizes  in  hrst  taper  period;  ah,N  indicates 
periodicity  in  mirror  region;  ah^i  indicates  periodicity  in  hrst  taper  period,  Wh  indi¬ 
cates  beam  width.  The  corresponding  parameters  are  similarly  introduced  for  the 
vertical  beam  (with  subscript  v).  The  thickness  of  both  beams  (in  the  z  direction)  is 
t.  Parameters  are  changed  linearly  inside  the  taper,  (c)  Field  pattern  of  Ey  for  cav¬ 
ity  mode  localized  by  horizontal  beam.  Parameters  are:  ah,N=4:53  nm,  a.y,Ar=272 
nm,  dhx,l/ dhx,N  dfiy^l/dfiy^N  0.5,  0‘h,l/(^h,N  0.7,  lh/cbh,N  f'2) 

lv/^v,N  0.83,  Wfij 1.65,  1.8,  dhy^j^fw^  dyx,Nl'^v  0.7,  d^x^N / 

dvy,N/av,N=0.5,  refractive  index  n  =  3.37,  with  slab  thickness  t/ah,N=0.35,  N=5,  and 
6  mirror  periods  for  both  beams.  Resonant  wavelength  is  1.55  pm  with  Q=19,000 
and  V=0.35(A/n)^.  (d)  Field  pattern  of  Ex  for  cavity  localized  by  vertical  beam. 
n  =  3.46  and  other  parameters  same  as  in  (c).  Resonant  wavelength  is  1103  nm  with 
Q=1900  and  V=0.47(A/n)^. 
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Figure  7.3:  (a)  Change  in  cavity  resonant  wavelength  and  Q  in  crossed  beam  structure 
for  cavity  mode  localized  in  horizontal  beam  as  cavity  length  is  varied.  Maximum  Q 
occurs  for  //j/a/i=1.25.  (b)  Change  in  horizontal  (long  wavelength)  cavity  resonant 
wavelength  and  Q  in  crossed  beam  structure  as  Wh  is  varied,  (c)  Change  in  vertical 
beam  quality  factor  and  frequency  conversion  hgure  of  merit  FOM  =  Q1Q2/ 
as  a  function  of  Wh-  (d)  Change  in  in-plane  quality  factor  in  x  direction  (i.e.  radiated 
power  collected  at  x  and  -x  edges  of  simulation  space)  for  resonance  localized  by 
vertical  beam  as  a  function  of  Wh  for  structure  shown  in  Fig.  7.2  (c)/(d). 
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Figure  7.4:  (a)  Dielectric  constant  and  (b)  FDTD-simulated  electric  field  Ey  for 
nanobeam  cavity  crossed  with  orthogonal  waveguide  for  Wy  =  0.75wh-  (c)/(d)  and 
(e)/(f)  show  the  same  for  increasing  width  of  the  orthogonal  waveguide,  illustrating 
cavity  confinement  weakens  as  orthogonal  waveguide  width  is  increased.  For  (c)/(d), 

Wy  =  Wh,  for  (e)/(f)  Wy  =  1.67wh. 
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Figure  7.5:  (a)  Crossed  nanobeam  cavity  design,  showing  intersecting  orthogonal 
nanobeams  with  taper  and  mirror  regions,  as  well  as  central  cavity,  (b)  Detail  of 
white  box  in  (a).  Parameters  used  to  form  resonance,  shown  for  cavity  in  horizontal 
(subscript  h)  beam:  R  indicates  cavity  length;  dhx,N  and  dhy,N  indicate  hole  sizes  in 
mirror  region;  dhx,i  and  dhy,i  indicate  hole  sizes  in  hrst  taper  period;  ah,N  indicates 
periodicity  in  mirror  region;  ah,i  indicates  periodicity  in  hrst  taper  period,  Wh  indi¬ 
cates  beam  width.  The  corresponding  parameters  are  similarly  introduced  for  the 
vertical  beam  (with  subscript  v).  The  thickness  of  both  beams  (in  the  z  direction) 
is  t.  Parameters  are  changed  linearly  inside  the  taper,  (c)  3D  FDTD  simulation  of 
held  pattern  of  Ey  for  cavity  localized  in  horizontal  beam  by  tapering  hole  dimen¬ 
sions  and  lattice  constant  in  central  region.  Parameters  are:  ah,N  =  C!,i;,Ar=449  nm, 
dhx,l/dhx,N  dfiy^x/dfiy^N  0.5,  Clh,l/(^h,N  (^v,l/(^v,N  0.7,  lh/Oih,N  R/(^v,N  f'4) 

Wh/ Ojh,N  Wx/ojv,n  1.65,  dhy^n/wh  diix,N /wy  0.7,  dhx,N / Oih  dyy^]\f/(iy^]\f  0.5,  re¬ 
fractive  index  n  =  3.37,  with  slab  thickness  t/a/i,Ar=0.35,  iV=8,  and  6  mirror  periods 
for  both  beams.  Resonant  wavelength  is  1.55  pm  with  Q=12,000  and  V=0.44(A/n)^. 
(d)  Field  pattern  of  E^  for  cavity  localized  in  vertical  beam.  Parameters  are  same  as 
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7.3.1  Linear  characterization  of  cavity  resonances 

We  now  present  characterization  of  cavities  with  nearly  degenerate  resonant  frequen¬ 
cies.  A  scanning  electron  microscope  image  of  structure  fabricated  in  GaAs  is  shown 
in  Fig.  7.6a.  The  structures  are  dehned  by  e-beam  lithography  and  dry  etching,  as 
well  as  wet  etching  of  a  sacrihcial  AlGaAs  layer  underneath  the  164  nm  thick  GaAs 
membrane.  To  characterize  the  resonant  frequencies  of  the  structures,  we  measure 
reflectivity,  using  vertical  incidence  through  an  objective  lens  and  free  space  coupling, 
in  the  cross-polarized  configuration  using  a  tungsten  halogen  lamp  as  a  broadband 
light  source[8].  The  principle  of  the  measurement  is  illustrated  in  Fig.  7.6b.  The 
input  and  measurement  polarizations  are  orthogonal;  the  cavity  polarization  is  ori¬ 
ented  45  degrees  to  both  output  and  input.  The  cavity  acts  as  a  frequency-selective 
polarization  rotator,  hltering  light  at  the  frequencies  of  the  cavity  resonances,  result¬ 
ing  in  a  high-signal  to  noise  measurement  of  the  cavity  reflectivity  from  which  the 
quality  factor  can  be  extracted.  The  result  of  this  measurement  for  the  structure 
in  Fig.  7.6a  is  shown  in  Fig.  7.6c.  There  are  two  resonances  visible  at  1571.2  nm 
(Q=4700)  and  1573.9  nm  (Q=7200).  The  difference  between  the  frequencies  of  the 
resonances  is  due  to  fabrication  imperfections;  improving  fabrication  or  employing 
local  tuning[114,  115,  116]  could  assist  in  yielding  perfectly  degenerate  resonances. 
Quality  factors  are  extracted  by  fitting  the  measured  data  to  the  sum  of  two  Fano 
hneshapes[117],  plus  linear  background  (to  account  for  the  measured  background  from 
the  white  light  source). 

7.3.2  Nonlinear  characterization  of  nearly  frequency  degen¬ 
erate  cavity  resonances 

This  cross-polarized  conhguration,  however,  cannot  easily  distinguish  between  cav¬ 
ities  with  vertical  and  horizontal  polarizations,  and  also  cannot  resolve  resonances 
with  frequencies  separated  by  less  than  a  linewidth.  To  more  sensitively  measure 
cavity  resonances,  we  use  intracavity  second  harmonic  generation[9,  104].  This  allows 
us  to  absolutely  determine  the  far-held  polarization  of  the  cavity,  since  second  har¬ 
monic  generation  is  maximized  when  incident  light  is  aligned  to  the  polarization  of 
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Figure  7.6:  (a)  Scanning  electron  microscope  image  of  crossbeam  structures  with  iden¬ 
tical  parameters  in  both  beams.  Structures  are  fabricated  by  e-beam  lithography,  dry 
etching,  and  wet  etching,  (b)  Experimental  setup  for  cross-polarized  reflectivity  mea¬ 
surements  to  characterize  cavity  resonances.  PBS  indicates  polarizing  beamsplitter. 
Cavity  polarization  is  oriented  45  degrees  {\H  +  V))  from  orthogonal  input  (|E)) 
and  measurement  {\H))  polarizations,  (c)  Cross-polarized  reflectivity  measurement 
of  structure  in  (a),  showing  two  resonances  at  1571.2  nm  (Q=4700)  and  1573.9  nm 
(Q=7200).  Solid  line  indicates  £t  to  sum  of  two  Fano  lineshapes. 


the  cavity.  Additionally,  because  the  measured  lineshape  from  the  second  harmonic 
is  proportional  to  the  square  of  the  lineshape  from  the  cavity,  we  can  more  easily 
distinguish  cavities  separated  by  small  differences  in  frequency.  A  second  harmonic 
measurement  of  a  different  structure  with  nearly  degenerate  resonances  in  shown  in 
Fig.  7.7.  Fig.  7.7a  shows  the  second  harmonic  counts  measured  for  three  different 
incident  laser  polarizations  as  the  laser  wavelength  is  varied.  The  angle  of  alignment 
between  input  polarization  and  the  polarization  of  each  cavity  mode  determines  the 
fraction  of  power  coupled  into  the  cavity,  and  accordingly  affects  the  magnitude  of  sec¬ 
ond  harmonic  generated.  The  intensity  of  second  harmonic  generation  for  each  mode 
also  depends  on  the  free  space  coupling  efficiency,  quality  factor  of  the  mode,  spatial 
overlap  of  fundamental  and  second  harmonic  helds,  and  detuning  of  laser  wavelength 
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from  the  cavity  resonance [9].  Two  modes  with  orthogonal  polarization  (with  resonant 
wavelengths  1541.75  nm  and  1542.1  nm),  as  expected,  are  clearly  visible.  To  further 
illustrate  this,  we  scan  additional  incident  polarizations  and  wavelengths,  as  shown  in 
Figs.  7.7b,  7.7c.  Fig.  7.7c  shows  2D  slices  at  three  incident  laser  wavelengths.  The  red 
lines  show  hts  for  two  resonances  with  polarization  separated  by  exactly  90  degrees, 
in  good  agreement  with  the  measured  data. 


(a)  (b)  (c) 


Angle  [deg] 


Figure  7.7:  (a)  Second  harmonic  characterization  of  structure  with  two  resonances 
nearly  degenerate  in  frequency  as  a  function  of  laser  wavelength  for  3  polarizations. 
Two  modes  with  orthogonal  polarization  are  visible,  (b)  Second  harmonic  intensity 
as  a  function  of  incident  laser  polarization.  Vertical  axis  indicates  wavelength  of 
laser;  horizontal  axis  indicates  angle  of  polarization.  Color  indicates  second  harmonic 
intensity.  Dotted  horizontal  lines  indicate  traces  in  (c).  (c)  Line  plots  of  second 
harmonic  generation  measured  at  different  polarizations  for  three  laser  wavelengths 
shown  in  (b).  Red  lines  indicate  hts  for  two  cavity  modes  with  polarizations  separated 
by  exactly  90  degrees. 


We  also  perform  nonlinear  mixing  with  two  resonant  modes  detuned  in  frequency. 
Fig.  7.8  shows  sum-frequency  generation  performed  in  a  structure  with  resonances  at 
1552.8  nm  and  1558.9  nm  (selected  to  have  two  resonances  overlapping  with  the  high 
power  range  of  our  tunable  lasers  and  large  enough  frequency  separation  to  observe 
the  SFG  peak).  The  central  sum- frequency  peak  is  smaller  due  to  the  smaller  spatial 
overlap  between  the  two  modes  than  with  a  single  mode.  There  is  no  cavity  resonance 
at  either  the  second  harmonic  or  sum  frequency. 


CHAPTER  7.  MULTIPLY  RESONANT  PHOTONIC  CRYSTAL  CAVITIES  77 


(a)  (b) 


Figure  7.8:  (a)  Schematic  of  sum  frequency  generation.  Light  from  two  CW  lasers  is 
coupled  into  two  cavity  resonances  at  1552.8  nm  and  1558.9  nm.  Nonlinear  frequency 
conversion  produces  light  at  the  second  harmonic  frequencies  of  each  laser,  as  well  as 
at  the  sum  frequency,  777.9  nm.  (b)  Sum-frequency  generation  from  structure  with 
resonances  as  indicated  in  (a). 


7.4  Fabrication  and  characterization  of  multiply 
resonant  photonic  crystals  with  large  frequency 
separation 

A  scanning  electron  microscope  (SEM)  image  of  a  fabricated  structure  is  shown  in 
Fig.  7.9a.  The  structures  are  defined  by  e-beam  lithography  and  dry  etching,  as  well 
as  wet  etching  of  a  sacrificial  AlGaAs  layer  underneath  the  164  nm  thick  GaAs  mem¬ 
brane.  To  compare  the  experimental  structure  with  the  proposed  design,  we  simulate 
the  fabricated  structure  by  converting  the  SEM  image  to  the  binary  refractive  indices 
of  air  and  GaAs[118].  The  thresholded  SEM  image  used  for  simulation  is  shown  in 
Fig.  7.9b.  Figs.  7.9c,  7.9d  show  the  simulated  electric  field  for  the  two  resonances  of 
the  fabricated  structure,  with  expected  resonances  at  1477  nm  (Q=4100)  and  1043 
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nm  (Q=540).  Finally,  we  experimentally  characterize  our  design  by  performing  a  re¬ 
flectivity  measurement  in  the  cross-polarized  con£guration[8]  (Figs.  7.9e,  7.9f)  using 
light  from  a  tungsten  halogen  lamp  which  is  linearly  polarized  using  a  Gian  Thompson 
polarizer  and  polarizing  beamsplitter.  The  cavity  is  oriented  at  45  degrees  to  both  the 
polarization  of  the  incident  light  and  orthogonal  polarization  used  for  measurement. 
We  measure  quality  factors  of  6600  at  1482.7  nm  and  1000  at  1101  nm.  Differences 
between  simulation  of  fabricated  structures  and  experimentally  measured  resonances 
are  most  likely  due  to  computational  error  in  binary  thresholding  of  the  refractive 
index.  Experimentally,  the  largest  separation  between  resonances  we  have  measured 
is  506  nm. 

Fig.  7.10  shows  cross-polarized  reflectivity  measurements  using  a  broadband  source 
of  structures  with  resonances  at  1546.6  nm  (Q=1600)  and  1023  nm  (Q=500),  a  fre¬ 
quency  separation  of  523.6  nm.  We  believe  the  experimental  Q  factors  are  limited  by 
fabrication  inaccuracies  caused  by  the  small  feature  size  of  the  central  holes  in  the 
vertical  beam,  which  are  located  in  regions  of  high  held  for  both  resonances.  We  could 
not  perform  nonlinear  characterization  of  these  structures  due  to  lack  of  availability 
of  a  laser  at  1023  nm. 


7.5  Improving  the  crossbeam  frequency  conversion 
platform 

For  highly  efficient  three  wave  mixing,  it  would  be  desirable  to  have  all  three  frequen¬ 
cies  confined  by  cavity  modes.  This  could  be  achieved  by  additionally  using  a  higher 
order  mode  in  one  beam.  However,  all  cavity  resonances  in  our  design  have  electric 
held  primarily  in  plane  (TE),  while  nonlinear  frequency  conversion  in  III-V  semicon¬ 
ductors  grown  on  (lOO)-oriented  wafers  requires  at  least  one  frequency  to  have  out  of 
plane  (TM)  polarization  (because  the  only  non-zero  elements  of  the  bulk  Xifk  have 
i  ^  j  ^  k[9]).  This  can  be  circumvented  by  using  a  wafer  with  different  crystal  orien¬ 
tation  such  as  (111) [119,  26]  so  that  the  effective  tensor  is  rotated  to  have  more 
non-zero  elements.  Fig.  7.11  shows  a  concept  illustration  of  highly  efficient  doubly 
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resonant  second  harmonic  generation  using  the  crossbeam  structure. 

Our  design  is  also  interesting  for  integration  with  single  semiconductor  emit¬ 
ters,  which  is  important  for  creating  quantum  photonic  interfaces  between  different 
quantum  emitters  and  telecommunications  wavelength.  As  described  in  the  previous 
chapter [12],  we  demonstrated  second  harmonic  excitation  of  a  InAs  single  quantum 
dot  emitting  at  900  nm,  creating  a  fast  single  photon  source  triggered  at  telecommuni¬ 
cation  wavelengths  (triggered  single  photons  at  100  MHz,  with  nonclassical  statistics 
visible  at  300  MHz)  .  The  speed  of  that  source  could  be  improved  by  using  a  second 
cavity  resonance,  as  in  our  crossed  beam  design,  to  enhance  the  spontaneous  emission 
rate  of  the  quantum  dot  through  the  Purcell  effect  [120],  with  possible  speeds  exceed¬ 
ing  1  GHz.  Finally,  our  design  is  also  promising  for  intracavity  frequency  conversion 
of  single  photons  from  an  integrated  quantum  emitter [103]. 

7.6  Conclusions 

In  conclusion,  we  have  proposed  and  demonstrated  the  design  of  a  photonic  crystal 
cavity  with  multiple  nearly  independent  resonances,  with  experimentally  measured 
Q>1000  for  both  resonances.  We  demonstrate  cavities  with  nearly-degenerate  fre¬ 
quencies,  as  well  as  frequencies  separated  by  more  than  500  nm,  which  we  characterize 
by  linear  and  nonlinear  spectroscopy.  Finally,  we  propose  extensions  of  this  work  to 
allow  high  efficiency  frequency  conversion,  which  might  be  used  in  conjunction  with 
integrated  single  emitters.  The  design  also  features  independent  waveguide  channels 
for  two  resonances,  which  facilitates  their  spatial  separation,  as  is  desirable  in  many 
applications.  Our  design  is  promising  for  resonantly  enhancing  on-chip  nonlinear  fre¬ 
quency  interactions,  such  as  stimulated  Raman  scattering  (which  could  be  achieved 
in  semiconductors  such  as  GaAs  or  Si  by  a  slight  detuning  between  the  parameters 
of  the  two  beams)  and  sum/difference  frequency  generation. 
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Figure  7.9:  (a)  SEM  image  of  cross-beam  structures  fabricated  in  164  nm  thick  GaAs 
membrane  with  ah,N=4:70  nm  and  a^,7v=320  nm  fabricated  in  164  nm  membrane, 
(b)  Thresholded  binary  image  of  SEM  used  for  simulating  fabricated  structure,  (c) 
FDTD  simulation  of  Ey  for  cavity  resonance  at  1477  nm.  (d)  Simulated  E^  for  cavity 
resonance  at  1043  nm.  (e)  Reflectivity  measurement  of  cavity  mode  at  1482.7  nm 
with  Q=6600.  (f)  Reflectivity  measurement  of  cavity  resonance  at  1101  nm  with 
Q=1000. 
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Figure  7.10:  Cross-polarized  reflectivity  of  crossbeam  structure  with  resonances  sep¬ 
arated  by  523  nm.  (a)  Resonance  at  1546.6  nm  (Q=1600).  (b)  Resonance  at  1023 
nm  (Q=500). 


Figure  7.11:  Concept  illustration  of  doubly  resonant  second  harmonic  generation  in 
photonic  crystal  crossbeam  nanocavity.  Incident  light  (red)  is  coupled  into  the  struc¬ 
ture  via  a  grating,  transmitted  to  the  cavity,  frequency  converted,  and  outcoupled 
through  a  separate  grating  (green). 


Chapter  8 

Visible  quantum  emitters  and  GaP 
photonic  crystal  cavities 


This  chapter  describes  two  such  systems  and  their  properties  when  coupled  to  GaP 
photonic  crystal  cavities,  described  in  Chap.  4.  The  high  refractive  index  of  the 
material  enables  a  large  photonic  band  gap  and  cavities  with  high  quality  factor, 
while  the  large  electronic  band  gap  prevents  absorption  in  the  near-IR  and  part 
of  the  visible.  The  hrst  system  consists  of  dinaphthoquaterrylene  diimide  (DNQDI) 
molecules  embedded  in  a  PMMA  thin  him.  This  him  can  then  be  deposited  on  a  GaP 
photonic  crystal  to  couple  molecules  to  cavity,  and  e-beam  lithography  techniques  can 
be  used  to  localize  the  molecule  to  the  cavity  spatial  position.  The  second  system 
is  InGaAs/GaP  quantum  dots  grown  by  molecular  beam  epitaxy;  when  coupled  to  a 
photonic  crystal,  the  outcoupling  is  sufficiently  enhanced  to  allow  the  observation  of 
narrow  spectral  lines  consisting  of  single  quantum  dot  emission. 

Beyond  this  thesis,  the  author  has  also  collaborated  on  relevant  work  with  other 
visible  quantum  emitters  and  GaP  photonic  crystals  including  nitrogen  vacancy  cen¬ 
ters  in  nanocrystals  [80],  Gs  atoms [121],  and  InP/InGaP  quantum  dots. 
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8.1  Motivation 

Coupled  photonic  crystal  cavity-emitter  systems  studied  so  far  are  primarily  based  on 
gallium  arsenide  and  silicon  materials,  which  absorb  strongly  at  wavelengths  shorter 
than  the  electronic  band  gap  of  the  material.  This  precludes  the  use  of  emitters 
atoms,  molecules,  defect  centers  in  crystals  such  as  the  nitrogen  vacancy  (NV“)  cen¬ 
ter,  as  well  as  many  quantum  dot  systems  having  room  temperature  conhnement, 
which  typically  have  resonances  at  visible  wavelengths.  Some  of  these  emitters  have 
particular  advantages  for  devices  resulting  from  long  electronic  and  spin  coherence 
times  and  room  temperature  coherence.  Second,  emitters  at  these  wavelengths  are 
easily  integrated  with  the  frequency  conversion  techniques  describing  earlier,  as  the 
frequency  of  these  emitters’  transitions  can  be  accessed  by  the  sum  frequency  of 
telecom-wavelength  lasers.  Additionally,  hybrid  systems  in  which  the  emitter  (such 
as  free-space  atoms  or  molecules  embedded  in  thin  hms)  and  cavity  exist  in  separate 
materials  allow  deterministic  positing  of  emitter  with  respect  to  cavity,  which  is  not 
possible  for  self-assembled  quantum  dots  grown  by  molecular  beam  epitaxy. 

8.2  Lithographic  positioning  of  molecules  on  high- 
Q  photonic  crystal  cavities 

As  described  previously  in  this  thesis,  photonic  crystal  nanocavities  have  been  used  to 
demonstrate  nanoscale  on-chip  devices  and  to  probe  fundamental  quantum  interac¬ 
tions  between  light  and  matter [45,  62,  122,  41].  Experiments  in  this  regime,  however, 
are  limited  by  the  precision  with  which  cavity  and  emitters  can  be  spatially  aligned 
and  by  the  spectral  range  of  emitters  that  can  be  coupled  to  cavity.  Emitters  are 
most  often  distributed  randomly  in  the  photonic  crystal  slab,  and  spatial  alignment 
to  the  photonic  crystal  cavity  occurs  by  chance.  Recently,  several  techniques  have 
been  developed  to  position  emitters  with  respect  to  cavities:  these  techniques  rely 
primarily  on  either  a  mechanical  transfer  process  to  bring  an  emitter  to  the  surface 
of  the  cavity[123,  124]  or  the  fabrication  of  a  cavity  at  the  location  of  a  previously 
detected  emitter [125,  126].  Neither  method  is  easily  scalable  to  arrays  of  cavities  and 
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emitters,  nor  achievable  with  conventional  semiconductor  fabrication  processes.  In 
this  section,  we  demonstrate  cavities  with  quality  factors  above  10,000  at  wavelengths 
compatible  with  near-IR  fluorophores  and  show  that  we  can  selectively  position  these 
molecules  on  top  of  a  nanocavity  using  conventional  lithography  techniques. 


8.2.1  Fabrication  of  cavity  and  thin  film  of  molecules 


(a)  (b)  (e) 


Wavelength  [nm] 


Figure  8.1:  (a)  SEM  image  of  a  fabricated  photonic  crystal  cavity  in  GaP.  Scale  bar 
indicates  200  nm.  (b)  FDTD  simulation  of  electric  held  intensity  of  the  fundamental 
cavity  mode.  The  mode  is  primarily  y-polarized  (c)  Schematic  illustrating  coupling 
of  molecule  to  cavity,  (i)  DNQDI/PMMA  is  deposited  over  the  entire  structure,  (ii) 
DNQDl/PMMA  is  lithographically  dehned  on  cavity  region,  (d)  Bulk  photolumines¬ 
cence  spectrum  of  DNQDl  when  excited  with  a  633  nm  HeNe  laser.  The  molecule 
has  a  peak  in  its  absorption  at  this  excitation  wavelength,  (e)  Chemical  structure  of 
DNQDl  molecule. 

The  cavity  is  a  linear  three-hole  defect  (L3)[33]  fabricated  in  a  125-nm  gallium 
phosphide  membrane  grown  by  gas-source  molecular  beam  epitaxy.  A  scanning  elec¬ 
tron  microscope  (SEM)  image  of  a  fabricated  cavity  and  the  simulated  electric  held 
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intensity  of  the  fundamental  high-Q  cavity  mode  are  shown  in  Fig.  8.1(a)  and  8.1(b). 
Cavities  are  fabricated  as  described  in  Chap.  4.  The  molecule  we  use  is  dinaphtho- 
quaterrylene  diimide  (DNQDI)  [127],  which  was  chosen  for  its  broadband  emission  over 
the  desired  wavelength  range  (700-850  nm),  good  photostability,  and  high  photolumi¬ 
nescence  (PL)  quantum  yield  (40%).  The  structure  of  the  molecule  and  its  emission 
spectrum  are  shown  in  Fig.  8.1.  To  couple  DNQDI  to  photonic  crystal  cavities  (Fig. 
1(c)),  the  molecule  was  dissolved  into  a  solution  of  1%  poly(methylmethacrylate) 
(PMMA)  in  distilled  toluene.  In  standard  lithographic  processing,  this  solution  is 
then  spun  onto  a  surface,  leaving  behind  a  smooth,  thin  him  of  dye-doped  resist. 
However,  spinning  onto  an  uneven  surface,  such  as  a  photonic  crystal  membrane, 
causes  unwanted  aggregation  of  the  dye-doped  PMMA.  Instead,  the  solution  was 
hoat-coated  [128],  whereby  the  photonic  crystal  sample  is  submerged  into  a  water 
bath  and  a  single  drop  of  the  dye-doped  PMMA  in  toluene  solution  is  dropped  onto 
the  surface  of  the  water  bath.  The  drop  quickly  disperses  across  the  surface  leaving 
a  locally  uniform  layer  of  hydrophobic  dye-doped  resist  hoating  on  top  of  the  water 
bath.  The  water  is  then  pipetted  away,  allowing  the  PMMA  layer  to  fall  on  top  of 
the  photonic  crystal  sample.  The  sample  is  baked  at  90°C  for  30  minutes  to  ensure 
that  all  the  water  is  fully  evaporated.  The  concentration  of  DNQDI  in  the  PMMA 
layer  is  approximately  5  molecules/ 100  nm^. 

8.2.2  Characterization  of  cavity  resonances 

We  hrst  characterize  cavities  passively  prior  to  depositing  molecules.  We  probe  cav¬ 
ity  resonances  using  cross-polarized  normal-incidence  reflectivity,  described  in  Ch.  4 
with  a  tungsten  halogen  white  light  source  [8].  The  cross-polarization  conhguration  is 
used  to  obtain  a  sufficient  signal-to-noise  ratio  to  observe  the  cavity  resonance  above 
the  reffected  background  uncoupled  to  the  cavity.  A  typical  reflectivity  spectrum  is 
shown  in  Fig.  8.2(a),  showing  the  multiple  resonances  of  the  L3  cavity;  the  fundamen¬ 
tal  mode  is  denoted  with  a  black  box.  The  spectrum  of  the  fundamental  mode  (Fig. 
8.2(b))  is  £t  to  a  Lorentzian,  giving  a  quality  factor  of  10,000.  (The  improvement  in 
quality  factor  from  [8]  is  due  to  better  fabrication.)  After  depositing  the  molecules 
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over  the  entire  structure,  we  measure  photoluminescence  of  the  molecule  (Fig.  8.2(c)) 
using  a  633-nm  helium-neon  excitation  laser  in  a  confocal  microscope  setup.  Above 
the  broad  emission  from  molecules  not  coupled  to  the  cavity,  we  observe  sharp  po¬ 
larized  resonances  identical  to  those  in  our  reflectivity  measurements,  demonstrating 
the  molecules  are  coupled  to  the  cavity  modes.  The  quality  factor  of  the  fundamen¬ 
tal  mode  is  measured  to  be  10,000,  indicating  that  deposition  of  molecules  onto  the 
membrane  does  not  degrade  the  properties  of  the  cavity,  in  agreement  with  hnite 
difference  time  domain  simulations  for  a  thin  (<40  nm-thick)  layer  of  PMMA.  Af¬ 
ter  deposition  of  molecules,  we  observe  a  small  (several  nm)  redshift  in  the  cavity 
resonance,  as  expected  from  simulations.  With  no  DNQDI/PMMA  present,  only 
background  counts  are  detectable  over  the  entire  spectral  range.  We  vary  the  spatial 
periodicity  of  the  photonic  crystal  holes  and  hole  radius  to  tune  the  fundamental  cav¬ 
ity  resonance  through  the  photoluminescence  spectrum  of  the  molecule.  We  measure 
high  cavity  quality  factors  up  to  12,000  via  photoluminescence  (Fig.  8.2(d))  across  a 
range  of  more  than  100  nm,  from  735  nm-860  nm.  The  cavity  Q  is  higher  at  longer 
wavelengths,  where  we  fabricate  most  of  our  cavities,  as  fabrication  imperfections 
are  reduced  because  the  feature  size  is  larger.  Small  differences  in  cavity  Q  between 
reflectivity  and  photoluminescence  measurements  (Fig.  8.2(d))  are  primarily  due  to 
fit  error. 

8.2.3  Lithographic  localization  of  molecules  to  cavity  region 

Since  the  molecules  are  doped  into  PMMA,  an  electron-beam  lithography  resist,  it  is 
straightforward  to  selectivity  expose  and  develop  the  resist  using  e-beam  hthography[129] 
so  molecules  and  PMMA  remain  only  at  the  location  of  the  photonic  crystal  cavity 
(Fig.  8.1(c)).  The  size  of  the  unexposed  region  at  the  center  of  the  photonic  crys¬ 
tal  cavity  is  approximately  700  nm  by  400  nm.  While  float  coating  deposits  resist 
uniformly  over  a  small  region,  PMMA  thickness  variations  were  observed  from  one 
coating  to  the  next,  so  electron  beam  doses  were  varied  for  different  cavities  on  one 
sample.  Fig.  8.3(a)  shows  a  scanning  confocal  image  of  photoluminescence  from  a 
photonic  crystal  cavity  coated  with  DNQDI-doped  PMMA.  The  photoluminescence 
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Figure  8.2:  (a)  Cross-polarized  reflectivity  measurement  of  a  cavity.  Box  indicates 
fundamental  cavity  mode,  (b)  Reflectivity  spectrum  of  high  quality  factor  fundamen¬ 
tal  cavity  mode  (box  in  (a)).  Spectrum  shows  additional  peaks  at  shorter  wavelengths 
from  higher  order  cavity  modes.  Solid  line  shows  Lorentzian  fit  with  quality  factor 
10,000.  (c)  Photoluminescence  collected  from  the  same  photonic  crystal  cavity  in 
(a-b)  after  molecules  are  deposited  on  cavity,  x-polarized  emission  is  shown  in  blue; 
y-polarized  emission  is  shown  in  red.  Inset:  PL  measurement  of  fundamental  cavity 
mode  (black  box).  Line  indicates  Lorentzian  fit  with  Q=10,000.  (d)  Quality  factors 
measured  from  reflectivity  before  molecule  deposition  and  photoluminescence  after 
molecule  deposition  from  the  high  Q  cavity  mode  for  structures  with  lattice  constant 
a  and  hole  radius  r/a  tuned  so  that  the  fundamental  cavity  resonance  shifts  across  the 
photoluminescence  spectrum  of  the  molecule.  Blue  open  circles  indicate  reflectivity 
measurements  for  the  cavities  that  were  also  measured  in  PL  (blue  closed  circles). 
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Figure  8.3:  (a)  Scanning  confocal  image  of  photoluminescence  from  DNQDI  doped 
PMMA  float-coated  onto  a  photonic  crystal  membrane.  Pixel  size  is  200  nm  and  scale 
bar  indicates  2  /im.  (b)  Scanning  confocal  image  of  DNQDI  PL  after  electron  beam 
lithography  is  used  to  remove  all  molecules,  except  for  the  ones  coating  the  cavity 
region  at  the  center.  The  same  imaging  laser  power  as  in  (a)  was  used.  Pixel  size  is 
80  nm  and  scale  bar  indicates  2  fim.  (c)  PL  spectrum  from  the  fundamental  mode 
of  photonic  crystal  cavity  after  selective  removal  of  molecules  by  e-beam  lithography, 
(d)  Atomic  force  microscopy  image  showing  localization  of  DNQDI-doped  PMMA  to 
the  cavity  region.  PMMA  thickness  is  12  nm.  Scale  bar  indicates  500  nm. 
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is  flat  to  within  3.5%,  with  slightly  more  emission  from  the  cavity  region,  likely  a  re¬ 
sult  of  enhanced  outcoupling  from  molecules  coupled  to  the  cavity  mode.  Fig.  8.3(b) 
shows  photoluminescence  from  the  same  cavity,  measured  with  the  same  excitation 
power,  after  electron  beam  exposure  and  removal  of  the  resist  surrounding  the  cavity. 
There  is  still  strong  emission,  though  diminished  by  the  e-beam  process,  from  the 
cavity  region,  but  there  is  no  emission  from  the  nearby  areas,  so  the  contrast  is  much 
larger.  Fig.  8.3(c)  shows  a  PL  spectrum  (Q=4500)  measured  on  the  same  cavity  after 
localization  of  the  resist  to  the  cavity,  demonstrating  that  molecules  are  spectrally 
coupled  to  photonic  crystal  cavity.  An  atomic  force  microscope  image  (Fig.  8.3d) 
conhrms  that  DNQDI-doped  PMMA  is  localized  to  the  cavity  and  is  12  nm  in  height. 
The  AFM  image  shows  a  misalignment  of  approximately  300  nm  between  the  cavity 
region  and  the  lithography  dehned  DNQDI/PMMA  region.  With  optimization  of  the 
overlay  process,  it  should  be  possible  to  reduce  this  error  to  less  than  50  nm. 

8.2.4  Low  temperature  spectroscopy  of  DNQDI 
(a)  (b)  (c) 


Figure  8.4:  (a)  DNQDI  photoluminescence  as  a  function  of  temperature.  For  temper¬ 
atures  below  50K,  narrow  lines  appear,  (b)  Individual  spectral  features  qualitatively 
match  the  expected  emission  from  a  zero  phonon  line  and  phonon  wing.  As  much  as 
30%  of  molecule  emission  is  measured  to  come  from  the  narrow  spectral  line. 

We  also  investigated  the  low  temperature  spectrum  of  the  DNQDI  molecules  un¬ 
coupled  to  a  cavity  to  assess  the  feasibility  of  achieving  a  stronger  coupling  between 
cavity  and  molecule.  At  room  temperature,  the  linewidth  of  the  molecule  emission  is 
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increased  due  to  phonons,  resulting  in  a  large  emission  bandwidth.  On  the  other  hand, 
if  the  spectral  linewidth  can  be  reduced  to  the  cavity  linewidth,  the  coupling  can  be 
enhanced  if  the  molecule  and  cavity  resonant  frequencies  are  matched,  enabling  an  in¬ 
crease  in  the  fraction  of  emission  into  a  particular  spectral  transition[130]  and  speedup 
of  the  emitter  spontaneous  emission  rate  [120].  If  the  coupling  is  sufficiently  large, 
the  strong  coupling  regime  can  be  reached [62]  where  emitter  and  cavity  form  mixed 
eigenstates.  Because  of  the  rotational  and  vibrational  structure  of  molecules,  for  co¬ 
herent  interactions,  it  is  also  important  to  access  the  zero  phonon  line  of  the  molecule, 
which  we  hope  to  observe  through  these  temperature-dependent  studies.  The  results 
of  these  measurements  are  shown  in  Fig.  8.4.  At  low  temperatures  (<50K),  we  in¬ 
deed  observe  narrow  spectral  lines  (Fig.  8.4(a)).  It  is  also  possible  in  some  cases  to 
resolve  individual  narrow  lines  that  appear  to  have  the  shape  of  a  zero  phonon  line 
and  phonon  wing  (Fig.  8.4(b)).  In  this  case,  it  appears  as  much  as  30%  of  the  emis¬ 
sion  may  come  from  the  zero  phonon  line.  We  measure  narrowest  linewidths  to  be  of 
the  order  of  100  Ghz  at  lOK,  as  expected  for  a  molecule  embedded  in  a  polymer  him 
[131].  We  also  observe  power  broadening  of  the  spectral  lines  (Fig.  8.4(c)). 

8.2.5  Spin-grown  crystalline  films  containing  DNQDI 

To  attempt  to  further  narrow  the  linewidth[131],  minimize  spectral  diffusion  [132], 
and  align  the  molecule  polarization,  we  also  investigated  embedding  the  DNQDI 
molecules  into  spin-grown  crystals  [133,  134]  tens  microns  in  lateral  extent  and  hun¬ 
dreds  of  nm  in  height.  To  form  these  crystals,  we  added  a  small  amount  of  DNQDI 
in  toluene  to  a  1%  p-terphenyl  (FT)  in  toluene  solution.  We  investigated  a  range 
of  experimental  parameters  including  spin  speed  (500-3000  rpm),  and  substrate  (Si, 
glass,  gold,  GaP).  We  also  studied  crystal  height  by  AFM  (Fig.  8.6(b)),  as  well  as 
the  crystal  luminescence  (Fig.  8.6). 

The  large  variation  in  crystals  formed  (Figs.  8.5,  8.6)  made  this  technique 

difficult  to  combine  with  nanophotonic  structures. 
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Figure  8.5:  Spin-grown  thin-film  crystals  of  DNQDI  molecules  on  different  substrates 
formed  by  using  different  spin  speeds,  as  indicated.  Scale  bars  all  indicate  10  fim. 

8.3  Photoluminescence  from  Ino.5Gao.5As/GaP  quan¬ 
tum  dots  coupled  to  photonic  crystal  cavities 

Semiconductor  quantum  dot  (QD)  emitters  grown  in  gallium  phosphide  are  impor¬ 
tant  for  both  classical  optoelectronic  and  quantum  applications.  The  close  match 
between  the  lattice  constants  of  GaP  and  Si  (0.37%  at  300K  [135])  is  promising  for 
monolithic  integration  with  silicon[135,  136,  137],  and  the  large  electronic  band  gap 
of  GaP  allows  light  emission  at  visible  wavelengths.  Single  quantum  dots  (QDs)  at 


CHAPTER  8.  VISIBLE  QUANTUM  EMITTER/CAVITY  SYSTEMS 


92 


Figure  8.6:  (a)  Optical  microscope  image  of  a  crystal  formed  by  spinning  at  3000  rpm 
on  a  glass  substrate,  (b)  AFM  image  of  crystal  spun  at  2500  rpm  on  glass,  indicating 
crystal  height  of  few  hundred  nm  and  lateral  size  of  pm.  (c)  Photoluminescence 
from  a  spun  crystal  excited  by  a  633  nm  HeNe  laser,  (d)  Confocal  photoluminescence 
scan  of  region  of  crystal  in  (c)  measured  by  a  Si  single  photon  counter,  c-axis  indicates 
photoluminescence  counts. 
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visible  wavelengths  are  beneficial  for  quantum  applications  since  Si  avalanche  pho¬ 
todiodes  (APDs)  have  maximum  quantum  efficiency  in  the  red  part  of  spectrum; 
additionally,  emission  in  this  part  of  the  spectrum  can  be  frequency  downconverted 
to  telecommunications  wavelengths  using  readily  available  lasers[138,  139]. 

Quantum  dots  emitting  in  the  red  have  been  extensively  studied  over  the  past 
decade  in  materials  systems  including  InP/InGaP[140,  141,  142,  143],  InP/GaP[144, 
145],  InP/AlGaInP[146, 147],  GaInP/GaP[148],  InAs/GaP[149],  and  AlGaInP/GaP[150]. 
Of  these  systems,  clear  single  quantum  dots  with  narrow  emission  lines  exhibiting  an¬ 
tibunching  have  been  observed  only  in  the  InP/InGaP  and  InP/AlGalnP  systems. 
GaP-based  materials,  by  contrast,  allow  either  monolithic  integration  with  Si  or 
growth  on  a  non-absorbing  GaP  substrate  (due  to  the  large  indirect  electronic  band 
gap);  additionally,  the  stronger  second  order  optical  nonlinearity  of  GaP  compared 
to  InGaP  is  preferable  for  on-chip  frequency  downconversion  to  telecom  wavelengths. 
Recently[79],  low  temperature  emission  (80K)  was  measured  from  Ino.5Gao.5As  self- 
assembled  QDs  in  GaP  emitting  in  the  red  part  of  the  spectrum.  This  system  provides 
large  wavelength  tunability  as  the  In  fraction  can  be  varied  from  0.07-0.50  without 
introducing  dislocations;  additionally  it  should  provide  deeper  confinement  for  car¬ 
riers  than  InP/GaP.  Subsequently,  room  temperature  emission  was  measured  from 
Ino.3Gao.7As/GaP  QDs[151];  measured  temperature  dependence  of  emission  and  sup¬ 
porting  tight  binding  calculations  indicated  good  confinement  of  carriers  and  type-I 
emission.  Here,  we  further  characterize  this  materials  system,  integrate  it  with  pho¬ 
tonic  nanostructures  that  enhance  the  emission  of  the  QDs,  and  observe  evidence 
indicative  of  emission  from  individual  QDs. 

8.3.1  Temperature-dependent  QD  luminescence 

The  QDs  are  grown  by  solid  source  molecular  beam  epitaxy  in  the  center  of  a  200  nm 
thick  GaP  membrane  grown  on  top  of  a  500  nm  layer  of  Alo.8Gao.2P  on  a  (001)  GaP 
substrate.  Fig.  8.7a  shows  the  measured  QD  photoluminescence  (PL)  as  a  function 
of  temperature  from  25K  to  300K  in  a  continuous  flow  helium  cryostat  using  700 
pW  excitation  power  from  a  405  nm  continuous  wave  (GW)  diode  laser.  (The  power 
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level  was  chosen  to  maintain  a  constant  integration  time  on  the  spectrometer  CCD 
for  all  temperatures.)  The  center  wavelength  of  the  QD  emission  redshifts  by  30  nm 
from  25K  to  300K.  The  large  full  width  half  maximum  of  the  emission  is  expected 
to  result  from  inhomogeneous  broadening  due  to  variation  in  the  physical  size  of  the 
QDs[79].  Fig.  8.7b  shows  that  for  higher  pump  powers,  the  QD  spectrum  broadens  on 
the  high  energy  side  and  the  integrated  intensity  is  nonlinear  as  a  function  of  power, 
indicating  the  presence  of  excited  states.  To  characterize  the  conhnement  of  carriers 
in  the  QDs,  we  study  the  intensity  of  QD  emission  as  a  function  of  temperature. 
Fig.  8.7c  shows  the  emission  intensity  integrated  across  the  low  energy  half  of  QD  PL 
spectrum  (to  minimize  the  contribution  of  excited  states)  from  Fig.  8.7a.  The  emission 
intensity  decreases  by  a  factor  of  4  from  cryogenic  temperatures  to  room  temperature. 
Fig.  8.7d  shows  a  £t  of  this  integrated  intensity  to  an  Arrenhius  model  (assuming  a 
temperature-independent  radiative  lifetime)  with  a  single  activation  energy  with  form: 


HT)  1 

lo  1  +  C  X  exp 


(8.1) 


where  I(T)  is  the  temperature-dependent  intensity,  Jq  is  the  intensity  at  OK,  C  is 
a  constant,  k  is  Boltzmann’s  constant,  and  Ea  is  the  activation  energy  indicating 
carrier  confinement.  We  measure  77^=161  meV;  this  is  slightly  larger  than  the  134 
meV  measured  by  Tranh  et  al[151],  most  likely  due  to  the  larger  indium  content  in 
our  samples,  which  results  in  deeper  confinement. 


8.3.2  QD  photoluminescence  dynamics 

We  investigate  the  dynamics  of  the  ensemble  QD  emission  by  studying  the  time- 
resolved  photoluminescence  on  a  streak  camera  when  the  quantum  dots  are  excited 
at  400  nm  by  a  frequency  doubled  Ti;Sapphire  laser  with  a  repetition  rate  of  80  MHz. 
The  experimental  setup  is  shown  in  Fig.  8.8a;  the  instrument  response  to  the  pump 
(12  ps)  is  shown  in  Fig.  8.8b  inset.  The  time-resolved  emission  of  the  low  energy 
half  of  the  QD  spectrum  (Fig.  8.8b)  for  all  temperatures  shows  a  biexponential  decay 
with  a  short  component  (pj  ~  250  ps,  averaged  across  all  temperatures)  followed  by 
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a  decay  with  similar  time  constant  (tj  ~  230  ps)  after  a  delay  of  ~500  ps.  The  delay 
is  most  likely  indicative  of  phonon-assisted  transfer  of  carriers  from  the  indirect  GaP 
matrix[144].  Fig.  8.8c  shows  the  time-resolved  PL  measured  at  250K,  indicating  re¬ 
gions  used  for  htting  initial  and  hnal  time  constants.  The  extracted  time  constants  for 
each  temperature  are  shown  in  Fig.  8.8d;  error  bars  indicate  error  from  £t.  The  short 
lifetime  is  consistent  with  a  Type  I  system;  the  minimal  temperature  dependence  of 
decay  rates  indicates  the  absence  of  temperature-dependent  non-radiative  processes. 

8.3.3  Enhanced  outcoupling  of  QD  emission  via  photonic 
crystal  cavities 

To  fabricate  photonic  crystal  cavities,  we  used  a  different  sample  with  a  thinner  93  nm- 
thick  GaP  membrane  for  ease  in  fabrication.  The  photonic  crystals  were  fabricated  by 
top-down  fabrication,  including  e-beam  lithography,  dry  etching,  and  wet  etching  to 
remove  the  sacrificial  AlGaP  layer [8].  A  scanning  electron  microscope  (SEM)  image 
of  a  fabricated  sample  is  shown  in  Fig.  8.9a.  Due  to  a  difference  in  MBE  growth 
parameters,  the  emission  wavelength  of  the  thinner  sample  was  slightly  blueshifted,  as 
shown  in  Fig.  8.9b,  and  the  QD  density  was  lower.  Fig.  8.9c  shows  photoluminescence, 
measured  at  12K,  from  the  quantum  dots  coupled  into  the  linear  three-hole  defect 
photonic  crystal  cavity  [33].  The  fundamental  mode  of  the  cavity  with  highest  quality 
factor  (black  circle)  overlaps  with  the  tail  of  the  QD  emission;  the  brighter  higher  order 
cavity  modes[76]  are  more  closely  matched  to  the  QD  emission  spectrum.  Fig.  8.9d 
shows  the  fundamental  cavity  mode  resonance  at  681.7  nm;  a  Lorentzian  fit  indicates 
a  quality  factor  of  2800;  Fig.  8.9e  shows  a  finite  difference  time  domain  (FDTD) 
simulation  of  the  spatial  distribution  of  the  electric  held  intensity  in  the  center  of  the 
membrane  for  this  cavity  mode. 
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(c)  (d) 


Figure  8.7:  (a)  QD  PL  from  an  unprocessed  region  of  sample  as  temperature  is 
varied  between  25K  and  300K.  The  pump  is  a  405  nm  CW  laser  diode  at  700  /iW. 

(b)  Photoluminescence  as  a  function  of  pump  power  at  lOK  showing  broadening  at 
higher  energies  with  increasing  power,  indicative  of  the  presence  of  excited  states. 

(c)  Integrated  counts  of  low  energy  half  of  PL  spectrum  (to  minimize  contribution  of 
excited  states,  which  varies  with  temperature).  Intensity  is  decreased  by  a  factor  of 
4  at  300K.  (d)  Semilog  plot  of  integrated  counts  of  low  energy  half  of  PL  spectrum 
versus  inverse  temperature.  Fit  to  exponential  (red  line)  gives  activation  energy 
Ea=lQl  meV. 
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Figure  8.8:  (a)  Experimental  setup  for  time-resolved  measurements.  3  ps  pulses  at 
800  nm  from  Ti:Sapphire  laser  are  frequency  doubled  in  a  BiBO  crystal.  A  grating 
is  used  to  filter  the  400  nm  light,  which  is  passed  through  a  dichroic  mirror  onto  the 
sample.  The  photoluminescence  emitted  by  the  quantum  dots  is  transmitted  through 
the  dichroic,  passed  through  a  long-pass  hlter  (LPF)  to  remove  any  residual  pump, 
and  sent  to  a  spectrometer  or  streak  camera  (for  time- resolved  measurements),  (b) 
Time-resolved  streak  camera  measurements  showing  lifetime  integrated  across  low 
energy  half  of  PL  spectrum  for  different  temperatures.  Inset:  excitation  pulse  at  400 
nm.  Red  line  indicates  ht  with  decay  time  12  ps,  limited  by  instrument  resolution, 
(c)  Exponential  hts  of  time-resolved  data  at  250K,  showing  initial  and  hnal  decay 
times  Ti  and  tj.  (d)  Measured  lifetimes  as  a  function  of  sample  temperature. 
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Figure  8.9:  (a)  SEM  image  of  photonic  crystal  nanocavity,  (b)  PL  measured  at  12K 
with  405  nm  CW  pump  from  unprocessed  region  of  thinner  sample  used  for  photonic 
crystal  measurements.  QD  wavelength  is  slightly  blueshifted  from  Fig.  1.  (c)  PL 
measurement  indicating  emission  of  quantum  dots  coupled  into  cavity  modes.  Fun¬ 
damental  cavity  mode  is  indicated  by  black  circle,  (d)  Lorenztian  £t  of  fundamental 
cavity  mode  at  681.7  nm  with  Q=2800.  (e)  FDTD-simulated  electric  held  intensity 
for  fundamental  cavity  mode. 
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8.3.4  Spectral  features  indicative  of  single  quantum  dot  emis¬ 
sion 

By  measuring  the  cavity-enhanced  QD  emission,  we  are  able  to  observe  indications  of 
single  quantum  dot  emission  lines  at  temperatures  below  60K.  Fig.  8.10a  shows  the 
low  temperature  PL  spectrum  with  a  CW  pump  from  another  photonic  crystal  cavity 
on  the  chip.  The  black  circle  indicates  the  higher  order  mode  of  interest  which  we 
use  to  enhance  the  quantum  dot  emission  outcoupling.  The  inset  shows  the  spatial 
prohle  of  electric  held  intensity  for  the  mode  of  interest,  calculated  by  the  FDTD 
method.  Fig.  8.10b  shows  a  high  resolution  spectrum  measured  at  lOK  of  the  QD 
and  cavity  mode  at  610  nm  (black  circle  in  Fig.  8.10a).  At  low  power,  two  narrow 
lines  appear  to  saturate  (as  would  be  expected  for  single  QDs)  as  power  is  increased 
above  about  30  fiW.  At  higher  powers,  we  recover  the  Lorentzian  lineshape  of  the 
photonic  crystal  cavity  mode,  as  the  intensity  in  the  cavity  mode  continues  to  grow 
while  intensity  from  the  individual  QD  lines  has  saturated. 

Fig.  8.10c  shows  the  wavelength  shift  of  quantum  dots  and  cavity  as  a  function 
of  temperature,  indicating  a  quadratic  redshift  in  dot  emission  as  temperature  is 
increased,  as  expected  due  to  the  approximately  quadratic  shift  in  material  band  gaps 
in  this  temperature  range[152].  Quantum  dot  lines  are  measured  at  500  nW  (far  below 
QD  saturation),  while  the  cavity  is  measured  at  100  /rW  (above  QD  saturation).  The 
quantum  dot  emission  wavelength  changes  at  a  different  rate  than  the  cavity  emission, 
conhrming  that  the  narrow  lines  are  not  associated  with  a  cavity  mode.  Fig.  8.10d 
shows  the  change  in  full  width  half  max  (FWHM)  of  the  observed  spectral  lines  at 
1  fiW  power  as  a  function  of  temperature  measured.  The  narrow  QD-like  lines  show 
an  increase  in  linewidth  as  the  temperature  is  raised,  as  expected  for  single  quantum 
dots,  while  the  cavity  linewidth  remains  roughly  unchanged  in  the  same  temperature 
range.  Further  conhrmation  of  single  QD  behavior  could  be  obtained  from  photon 
statistics  measurements.  We  did  not  obtain  sufficient  signal-to-noise  from  the  cavity 
to  perform  such  measurements  in  this  case.  An  improvement  in  the  signal-to-noise, 
for  example  by  improving  the  cavity  quality  factor,  would  also  allow  an  investigation 
of  the  time-resolved  dynamics  of  a  single  QD  coupled  to  the  cavity,  where  Purcell 
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Figure  8.10:  (a)  PL  from  a  different  photonic  crystal  cavity  structure  with  higher  order 
cavity  mode  aligned  to  wavelength  of  strongest  QD  PL  measured  with  405  nm  CW 
pump.  Black  circle  indicates  mode  of  interest.  Inset:  FDTD-simulated  electric  field 
intenstiy  for  higher  order  cavity  mode  indicated  by  circle  in  (a),  (b)  High  resolution 
spectra  of  QD  PL  from  cavity  in  (a).  Measurements  are  taken  at  lOK.  The  two  single 
lines  at  low  power  are  indicative  of  single  quantum  dots,  (c)  Change  in  wavelength  of 
quantum  dots  (measured  at  500  nW)  and  cavity  (measured  at  100  /iW)  as  a  function 
of  temperature.  Solid  lines  indicate  quadratic  fits,  (d)  Change  in  QD  linewidth  as  a 
function  of  temperature. 
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enhancement  is  expected [120]  in  this  regime. 

8.3.5  Conclusions 

In  conclusion,  we  have  studied  the  properties  of  two  photoluminescent  quantum  emit¬ 
ters  (DNQDI  molecules  and  InGaAs/GaP  quantum  dots)  at  visible  wavelengths  cou¬ 
pled  to  GaP  photonic  crystal  cavities. 

We  have  demonstrated  the  coupling  of  DNQDI  fluorescent  molecules  and  pho¬ 
tonic  crystal  cavities  with  resonances  in  the  far-red  and  near-infrared  wavelengths 
and  quality  factors  up  to  12,000.  By  exposing  and  developing  the  molecule’s  polymer 
host  using  electron  beam  lithography,  we  localize  the  molecule  to  the  cavity  region. 
Our  results  show  that  molecules  can  be  coupled  to  high  quality  factor  photonic  crys¬ 
tal  cavities  and  easily  localized  to  the  spatial  location  of  the  nanoscale  cavity  using 
standard  lithographic  techniques.  We  also  investigated  the  low-temperature  photo¬ 
luminescence  of  DNQDI  and  performed  preliminary  characterization  of  spin-grown 
thin  him  crystals  with  DNQDI  embedded. 

We  have  also  measured  the  temperature-dependent  photoluminescence  emitted 
by  Ino.5Gao.5As  quantum  dots  embedded  in  a  GaP  membrane,  indicating  good  car¬ 
rier  conhnement  with  only  a  four-fold  decrease  in  emitted  intensity  from  cryogenic 
to  room  temperature.  We  studied  the  temperature-dependent  time-resolved  photo¬ 
luminescence,  showing  a  biexponential  decay  with  time  constants  of  ~  200  ps.  We 
observed  enhanced  emission  into  the  modes  of  a  photonic  crystal  cavity  and  narrow 
lines  consistent  with  single  quantum  dot  emission.  The  materials  system  is  compati¬ 
ble  with  monolithic  integration  on  Si  and  is  also  promising  for  quantum  applications. 
The  quantum  dot  wavelength  is  matched  to  the  high  efficiency  region  of  silicon  APDs 
and  could  be  downconverted  to  telecommunication  wavelengths  through  integration 
with  photonic  nanostructures [9,  11,  10,  103]. 


Chapter  9 
Conclusions 


This  dissertation  has  presented  a  framework  for  three-wave  mixing  devices  in  III-V 
semiconductor  photonic  crystals.  This  enables  miniaturization  of  nonlinear  frequency 
conversion  devices  by  reducing  the  required  optical  power  and  device  size  and  allowing 
integration  with  semiconductor  optoelectronic  materials. 

The  main  results  presented  were: 

•  The  hrst  demonstration  photonic  crystal  cavities  in  gallium  phosphide 

•  The  demonstration  of  430%/W  (10“^  for  11  /iW  coupled  power)  continuous  wave 
second  harmonic  generation  in  gallium  phosphide  photonic  crystal  cavities,  an 
8-orders  of  magnitude  improvement  on  previous  work  in  InP  [7] 

•  The  first  demonstration  of  continuous  wave  sum-frequency  generation  in  pho¬ 
tonic  crystal  cavities,  fabricated  in  GaP 

•  The  first  demonstration  of  second  harmonic  generation  in  photonic  crystal 
waveguides,  fabricated  in  GaP 

•  The  demonstration  of  the  fastest  optically  pumped  single  photon  source,  based 
on  a  quantum  dot  pumped  by  intracavity  second  harmonic  generation  from  a 
telecom-wavelength  laser 
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•  The  design,  fabrication,  and  characterization  of  the  largest  frequency  separa¬ 
tion  multiply  resonant  photonic  crystal  cavities  (more  than  600  nm  resonant 
wavelength  separation  measured  experimentally) 

•  The  first  demonstration  of  lithographic  techniques  to  align  thin  dims  of  molecules 
to  high-Q  photonic  crystal  cavities 

•  Spectral  and  temporal  characterization  of  InGaAs/GaP  quantum  dot  lumines¬ 
cence,  including  the  first  observation  of  narrow  spectral  lines  indicative  of  single 
quantum  dots 

9.1  Outlook 

This  work  could  immediately  be  built  upon  primarily  by  integrating  and  extending 
the  multiply  resonant  design  presented  in  Ghap.  7  to  improve  upon  the  work  presented 
in  earlier  chapters.  The  design  and  fabrication  of  these  structures  could  be  extended 
to  allow  even  larger  frequency  separation  to  enable  triply-resonant  three-wave  mixing. 
By  fabricating  these  multiply  resonant  structures  on  (lll)-oriented  GaAs  or  GaP,  it 
should  be  possible  to  improve  the  efficiency  of  second  harmonic  and  sum-frequency 
generation,  and  also  allow  the  observation  of  difference  frequency  generation,  which  is 
a  prerequisite  to  optical  parametric  oscillation  and  spontaneous  parametric  downcon- 
version  (which  can  generate  entangled  photon  pairs).  All  the  work  presented  in  this 
dissertation  could  be  scaled  to  longer  wavelengths  deep  into  the  infrared  to  provide 
tunable  light  sources  at  frequencies  critical  for  molecular  and  gas  sensing[153]  This 
cavity  design  is  also  promising  for  enhancing  third-order  nonlinear  optical  frequency 
conversion  processes  such  as  four-wave  mixing  and  stimulated  Raman  scattering. 

Furthermore,  employing  the  multiply  resonant  cavity  design  with  resonant  fre¬ 
quencies  at  a  frequency  in  the  G-band  and  at  the  frequency  of  a  single  InAs  quantum 
dot,  could  improve  the  rate  of  triggered  single  photon  generation  beyond  10  GHz 
for  cavities  with  moderate  Purcell  factor,  which  would  surpass  the  fastest  electrically 
pumped  single  photon  source [154,  155]. 

Finally,  the  quantum  emitters  at  visible  frequencies  discussed  in  Ghap.  8  are 
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promising  for  downconversion  of  single  photons  to  telecommunications  wavelengths, 
which  could  also  be  integrated  with  single  photon  pulse  shaping[156,  157]. 


Appendix  A 

Fabrication  of  photonic  crystal 
structures 


Photonic  crystal  structures  were  fabricated  in  an  epitaxially  grown  structure  on  III-V 
wafer  pieces.  The  process  is  as  detailed  below,  and  is  illustrated  schematically  in 
Figure  A.l.  Starting  in  2010,  the  details  of  each  fabrication  step  for  each  sample  were 
recorded  in  a  spreadsheet  accessible  from  the  cleanroom  as  the  process  took  place, 
facilitating  detailed  records  as  well  as  information  transfer  between  collaborators. 
Fabrication  in  GaP  and  GaAs  was  performed  using  similar  process;  the  main  difference 
was  a  longer  etch  time  for  GaP  compared  to  GaAs. 

A.l  Wafer  preparation 

Prior  to  coating  of  wafer  pieces  with  resist,  samples  were  optically  inspected  using  a 
microscope  and  then  sonicated  5  minutes  in  acetone  to  remove  organics. 

A.  2  Lithography 

Photonic  crystals  were  fabricated  using  a  top-down  approach  starting  with  electron 
beam  lithography.  In  some  cases,  an  HDMS  prime  was  applied  hrst  using  a  YES  oven 
to  enhance  adhesion  of  resist  to  the  samples.  The  positive  e-beam  resist  used  was 
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1 ,  Mask  deposition 


2.  E-beam  lithography 


AIGaP 


3.  Dry  etching 


mil 

AIGaP 


GaP 


4.  Resist  removal 


5.  Undercut 


■  IIIIIH 

AIGaP 


GaP 


■  IIIIIH 


GaP 


Figure  A.l:  Top  down  fabrication  process  for  photonic  crystal  structures,  including 
spinning  e-beam  resist,  e-beam  lithography,  dry  etching  into  III-V  semiconductor 
membrane,  and  wet  etch  to  remove  sacrihcial  Al-containing  layer  and  leave  suspended 
high-index  contrast  membrane  structures. 


ZEP520  in  an  anisole  solvent.  The  ZEP  was  removed  from  the  refrigerator  30  minutes 
prior  to  spinning  in  order  to  warm  to  room  temperature.  The  thickness  of  the  layer 
was  varied  between  250  nm  and  375  nm  for  different  experiments,  according  to  thick¬ 
ness  and  composition  of  the  semiconductor  membrane;  this  corresponds  to  a  range  of 
spin  speeds  from  about  4000  rpm  to  8000  rpm.  Prior  to  spinning  cleaved  pieces  of  III- 
V  wafer  pieces  containing  layered  epitaxy,  a  silicon  piece  was  hrst  coated  with  ZEP. 
The  thickness  of  this  silicon  piece  was  then  measured  on  a  nanospec  tool  to  conhrm 
that  the  ZEP  behaved  according  to  the  expected  spin  speed  curve.  The  III-V  wafer 
pieces  were  coated  on  the  backside  with  blue  sticky  tape  prior  to  spinning  to  prevent 
resist  from  accumulating  on  the  backside  of  the  wafer.  Spun  pieces  of  Si  test  pieces 
as  well  as  III-V  wafer  piece  samples  were  baked  for  2  minutes  at  200°C  to  evaporate 
the  remaining  solvent  and  then  visually  inspected  using  an  optical  microscope. 
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Exposure  of  the  ZEP  mask  was  performed  using  e-beam  lithography.  For  most  of 
the  work  in  this  thesis,  a  10  kV  e-beam  with  7.5  pm  aperture  (Raith  150)  was  used; 
for  later  work  (including  parts  of  Chap.  7),  a  100  kV  e-beam  tool  (Jeol  JBX  6300) 
was  used  to  reduce  lateral  proximity  effects.  Each  sample  was  cleaved  to  a  size  of  a 
few  mm  per  side,  allowing  hundreds  of  devices  to  be  packed  onto  a  single  sample. 

A  standard  development  procedure  was  followed:  40  s  in  Xylenes,  followed  by  a 
40  s  rinse  in  a  mixture  of  25%  methyl  iso  butyl  ketone  and  75%  isopropanol,  followed 
by  a  30  s  rinse  in  isopropanol. 

ZEP  was  used  rather  than  the  more  standard  PMMA  due  to  its  enhanced  etch 
resistivity. 

A. 3  Dry  Etching 

Dry  etching  from  exposed  ZEP  mask  into  GaP  or  GaAs  membrane  was  performed 
using  a  Cl-based  etch  in  a  PlasmaQuest  electron  cyclotron  resonance-reactive  ion 
etcher  (ECR-RIE).  Prior  to  etching,  the  chip  was  glued  to  a  Si  carrier  wafer  using  a 
drop  of  PMMA,  then  baked  at  75°C  for  2  minutes  to  remove  remaining  solvent.  A 
standard  recipe  for  etching  is  as  follows: 


Step  1  (setup) 

Step  2  (ignition) 

Step  3 

Ar  (seem) 

15 

15 

15 

BGI3  (seem) 

10 

10 

10 

GI2  (seem) 

3 

3 

3 

Pressure  (mT) 

2 

2 

2 

EGR  Power  (W) 

0 

400 

200 

RE  Power  (W) 

0 

40 

40 

He  pressure  (mT) 

10 

10 

10 

Temperature  (°G) 

13 

13 

13 

Time  (s) 

70 

10 

250-600 

The  etch  rate  of  GaP  in  this  recipe  was  found  to  be  ~2  nm/s;  the  etch  rate  of 
GaAs  was  approximately  1.5  times  larger. 

As  operating  conditions  varied  in  the  PlasmaQuest  varied  over  the  years  due  to 
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various  state  of  cleaning  and  repairs,  as  well  as  different  etches  by  users,  a  few  etch 
parameters  were  adjusted  to  optimize  etching.  In  particular,  CI2  flow  was  varied  from 
1.5-3  seem  and  RF  power  was  varied  in  the  range  of  40-50W.  This  value  of  RF  power 
corresponded  to  a  DC  bias  of  -50  to  -40V.  Starting  from  April  2010  (in  measurements 
done  with  Sonia  Buckley),  prior  to  etching  photonic  crystal  sample,  a  sample  of  ZEP 
spun  on  silicon  was  hrst  etched  for  100  s.  The  thickness  of  the  ZEP  was  measured 
before  and  after  this  etch,  and  the  average  etch  rate  was  calculated.  A  high  ZEP  etch 
rate  would  constitute  a  flag  in  the  process,  indicating  a  problem  with  the  state  of  the 
machine.  The  backside  helium  cooling  for  the  etches  varied  in  the  range  of  3  sccm-10 
seem.  Over  the  course  of  1.5  years  of  measurement,  the  ZEP  etch  rate  was  found  to 
be  0.3±0.2  nm/s  (Fig.  A. 2);  in  general,  a  ZEP  etch  rate  of  >0.5  nm/s  indicated  a 
problem. 


etch  rate  (nm/s) 


Figure  A. 2:  Etch  rate  of  ZEP  on  Si  using  PlasmaQuest  GaP/GaAs  etching  recipe 
over  76  runs  from  April  2010  to  December  2011.  Etch  rate  was  measured  as  0.3±0.2 
nm/s. 


After  etching,  the  remaining  resist  was  removed  by  one  of  three  methods:  (1) 
sonication  in  acetone  (2)  oxygen  plasma  (3)  overnight  soak  in  1165.  We  found  option 
3  (1165)  to  be  the  most  robust  and  cause  the  least  damage  to  structures. 
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A. 4  Wet  etch  for  undercutting 

A  hydrofluoric  acid  wet  etch  was  used  to  remove  the  sacrificial  aluminum-containing 
layer  underneath  the  membrane  containing  the  photonic  crystal.  For  80%  A1  in 
AlGaAs  and  85%  A1  in  AlGaP,  a  7%  HF  etch  was  effective.  For  higher  A1  content, 
1-2%  HF  was  used.  The  HF  wet  etch  was  performed  by  alternating  few  second  dips 
in  HF  and  water,  in  order  to  keep  A1  from  outgassing  too  quickly  and  breaking 
membrane  structures.  Occasionally  a  hnal  1%  HF  dip  or  citric  acid  dip  was  used  to 
remove  remaining  A1  crystals. 

If  a  contaminant-free,  clean  surface  was  desired  (as  in  experiments  with  DNQDI 
molecules  in  thin  dims  in  Ghap.  8),  a  few-minute  oxygen  plasma  clean  using  standard 
recipe  oxygen. rep  was  performed  after  wet  etching. 

Figure  A. 3  shows  a  single  chip  as  it  passes  through  each  step  of  the  process  from 
spinning  e-beam  resist  through  undercutting. 
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Figure  A. 3:  (a)  Optical  microscope  image  of  ZEP  spun  onto  GaP.  Defects  from  the 
original  wafer  propagate  to  the  epitaxially  grown  layers  and  create  perturbations  in 
the  ZEP  him.  (b)  Optical  microscope  image  of  structures  after  e-beam  exposure  and 
development,  (c)-(d)  Optical  microscope  images  of  structures  after  400s  dry  etch. 
(e)-(g)  SEM  images  of  structures  after  HE  undercut. 
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